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Cl.l HYBRID2 


1 . ProRram Description 

HYBRXD2 computes the fuel and energy consumption of a hybrid 
vehicle with a bi-modal control strategy over specified component 
driving cycles. Fuel and energy consumption are computed 
separately for the two modes of operation. The program also 
computes yearly average fuel and energy consumption using a 
composite driving cycle which varies as a function of daily travel. 
The modelling techniques used include the following: 

Heat engine - represented by a map of bsfc as a function of 
bmep and rpm, together with a curve of maximum torque 
versus rpm. The displacement of the engine for which this 
data is supplied is used as input; the program has provisions 
for scaling the data to other displacements. 

Electric motor/controls - electrical input represented as 
a constant load plus shaft power divided by a fixed efficiency 
in both driving and braking modes. Maximum (driving) and 
minimum (braking) torque as functions of rpm are also required 
Batcery - modelled by a fractional depletion technique using 
the power averaged over a specified time interval rather 
than instantaneous power. 

Engine accessory load - represented by a curve of torque 
required vs. system output (torque converter input) rpm. 
Included in these is the transmission front pump, in addition 
to belt-driven accessories. 

Torque converter - represented by curves of speed and torque 


C-2 


ratios (output/input) as functions of an output speed-torque 

parameter equal to output speed/output torque. An Input speed- 

2 

torque factor, (input speed) / (input torque), at stall must also 
be specified. 

Gearbox - represented by a set of gear ratios with different 
efficiencies for each ratio. Spin loss coefficients 
(exclusive of the front pump) may also be specified. 

Differential - same treatment as gearbox. 

Vehicle road load - represented by a combination of an 
aerodynamic load (proportional to speed squared) and tire 
rolling resistance. The rolling resistance coefficients 
can include a constant term and one which is linear with 
vehicle speed. 

The program structure is modular, with the control strategy 
and shift strategies being contained in separate subroutines. 

A breakdown of the program routines and their functions is as 
follows : 

1. 1TYBRID2 (Main Program) 

Input of case data 
Output 

- Numerical integration 

Computation of yearly average fuel and energy 
consumption from individual driving cycle results. 

2. VEHIC 

Computation of road loads, power flow through the 
vehicle system up to the torque converter output. 
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Coraputatlon of derivatives of all variables of 
integration. 

3 . HYP£AD 

Input of fixed, detailed component data. 

4. GRSHFT 

Controls transmission gear ratio in accord with a 
pre-set shift strategy. 

5 . PMOVR 

Controls heat engine/motor power split in accord 
with a pre-set control strategy. 

6. TQCON 

Computes torque converter input speed and torque 
given output conditions, or output torque and 
input speed given input torque and output speed. 

7 . FILTER 

Filters battery output pcwer to provide a smoothed 
battery output power curve. 

8. INTI, INT2 

One and two dimensional interpolation routines, 
respectively. 
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GCLETA FCCTPftK 1,3 • *;F“I-AL’TP 


(Ct-tO-73) 


C-39 


C 

c 

c 

c 

OQOfI03 - - 


000003 

c 

- - c 

00 0003 


000003 

000003 

000003 

000003 

000003 

C00003 

000003 

000003 

000003 

C 

000003 

C 

c 

c 

000005 

000013 

000015 

00002 ? 

0 0 c 0«40 
000056 
00C060 
I? ■'65 
0 . j 6 ? 

0 0 0 0^2 
00 C 07 C 
3 " C 0 7 6 


F'?c^!paM^'v^pTc? 

‘^CCELS fP^CATT^.^ C^' HY6PTO CV'=’9 SP'^'CI^IEC CkIVI 

CYCLES. 

CATA !»■' CC^'^'C> - 

COHMCM CI5 = L,FPK(?0) ,7nYAX(?C) ,eYPP(?0),?SFC(?0,20),K9PY;MyEF, 
lFMUMtEVUG,FINMr,SPYC(?0 ),TMCYAY(?0) ,TMOMTN(20) ,NSF“0, 
2CTC?F0,T5F(20),TQO(2n)»SFR(20),TSP?(?0)» NTSP » 
3NL0CK(5),TPATTC(5) ,CT1(5),CT2(5),F‘"UT(5) , 

^ W3 ,FBMAX,F>"L'FG ,EYUPG2 , 

5CRATIO,C01,rr2,FY.UC,FTTFE,CTIREl,CTIR£2,V-ASE,nLT,CCA, 

eTFC^TN,rRMAy;RPVTOL»'yYASF2tFFOr'IK, 

7PIJFL SG ,SKALE ,TMFXL ,MH,NCDSCH ,CHGF FF , SC ycLE , NC C HP , CS TA V , 
8NUS'ITS,TH5FT(5),UP'£HFT(5),Cr:FHFT(5),ODISCH{20),CYCLES(20)',MC(3 
PKFPTr(3) »CTC(3) ,TFC(3),rSL'P(30),CSC(30),D3AR(30) 
l,RF“BST»pPFF7tPFGP(20),RFHFCP(20),RD“OPM , E^uC'yT ,RATL|o,®^7GK',K'OF 
2 , PCEK'S ( 20 ) , ETENS (20 ) , KCF ^S 

3tVMAy,TFOHN'?,UPSEL»CSEEL,CWA9H,NAX,RP«AX(?0),TAX(20) 

CCHMCN TTHC(3,2C0),SFECC(3,?00),GAHMA(in,3) 

VARIABLES IF CCYMCN- - 

COHMCN' FA,FR,rAC.F*:FT,TCC,TTC,RP»"C0,RP“''Tr',TSC»RO*-3C»RO'^YC0,TTCC 
IRPHTCNt THr»7EC ,TKC2, 7EC 2 , T G C » P SC , PE 0 , <=50 2 , C , P H 0 2 t^GC , FPR^ ,FC , 
2°0 ♦P7, OFW,PL7C»7HRFC ,T7C IN' 

3,Ro*^6C,FP'^EC2,FI yc2,V“FS 

t* »7F^’,'’P^'SC?,FS^?,PPH7C1 ,RFH7 C?»pR “E°»°RHFp? 


= EAL Ync7(’0),Y(7n),TVTi,isT?,V(r:) ,A(3»,Y7'<'.n(^0) 
real FY(fiO) 

'"E AL DTS7 ( 7 1 ,FCrsS ( 3) rss? (3 ) ,FCOMS ( 3) ,ECCvS? ( 

real FrBAF(2D),PFPAC(Tn),TIHF(?C0),SPZFr'(2nn), 

lPCliP?(30) ,FrFAPC»0) ,rf'PAR2(30)»PASGr(ini ,PC'^A'> (30) ,FCYAV(30) 
RE AL ECHF f3)«rC'^YF(3),FHF?AP(30)»ESYSPP(3n) 

”EAL VAVG(3) 

peal V=^A R ( 25 ) t SP = WR ( 25) , S°FNG ( 25) 

REAL =>^CT(7en0),ERG(2?0n),CEPL(.3)tRC'’AR(30) 

INTCGFR IFO(20) 


CALL HYPEAC (NCPfiO) 


9L'»’ CAT A 
READ 60C,^CASE 


original page is 

OF POOR quality 


DO IRA JCAGF=i,kcASF 
oRI?,’T ??5,JCASE 

REAP 5in,7ECHIN,TECHS2»C5'^'AX,PEP“IM,V«AX 

PRINT R30,TecYTN,TFC^'K2,r5HAX,CFnv|I>.’,V''AY 

VMAX = VYAX/3 ,6 

READ 610,C7FL7o 

no IRQ JCYC=1 iNCYCLF 

PRINT 700 

NT i“c = f Tc (jrYn 


HPRS'Tr*'ri'T^ ( jrYC) 


PE LT='^TC ( JC fC ) 


C-40 


GCLETft 


onoo7? 

p'’ f)l 01 

lOE 
ODO 110 

nnoii4 

000117 

C00121 

000122 

000123 


000125 
000127 
P00130 
000121 
000133 
00013U 
000135 
000136 
000137 
CIOO 1^*0 


000141 
142 
I. . j 1 4 3 
OOOlrl 
00015^ 
900163 

onnies 
000166 
000170 
000171 
000 17^* 
000175 
000203 
000213 
000 217 
000222 
000226 


00C231 
000232 
0Q0 2M 
000243 
000247 
0 " 25 ^* 
c 256 
000261 
000263 
000265 


FC‘^7F«K l,r * /'•L'TC ccl • (01-10-7T) HY55T02 


’F='^FC ( JC'^C ) 

^CT=IFiy(TF/CFL7+2,*(CTFLT2/1.0)+l.) 

Cn 1*7 J = 1,K'7IP5 
Tlf^E rJ)=TI^C (JCYC, J) 

SPEEC(J)=‘^FECC(JCYC,J1 
’ t<7 CO^'TI^■^^E 

I — ■ no 4 8 J=1 ,MCT 
■ -- PHOT (J) = 0. 0 
-- i*8 CCNTIMUE 
C ' 

C initialize V6FIAPLES CP INTEGRATION 
C 

no 50 J=I,2C 

Y ( J) =0. 0 
I80( J) =0 

“~ro ccntinl'c 

T = 0. 0 
ITIBP=0 
KrNPPNT 
JGFAR=1 
JGEAR2=1 
NTH = 0 
C 

C COMPUTE VPLCCITy ANR ArrPLPPATICN 

C 

V ( 1 ) = 0 , 0 
'M3) =3^550 ( 1 ) 

Vl'O = INTI (CELT/2. »TIVf’, 3 C55C,NTI‘^5) 
V(5)-T*''Tl(rELT;TiMF,3r;ppr,sTTMP) 

'M6-) =TNT1 (3 . *rELT/2. , TTVF, SPPPC. VTT^E) 

V { 7) ='/ ( 3) + V (■’J-V ( £*) 

GC YC 70 

6 0 V ( 1) = V r 3> 

'/(3)=V(ll) 

V(2)=(V(l)*V(3))/2. 

V (4) ='/ (6 ) 

V(5) =IKTl (T + CELT.TTRr.EPPEC.NTTM'^) 

V(6) =TMT1 (T*3.»r5LT/2. ,'^T^^E.S0 55P.nTIMri 

70 A ( 1) = C V (U) - V (?) ) / (3, e »^ELT) 

Ad) =A,mtni (A (IdR.") 

A(2) =(V{f )-V(3) ) /(3.6^0rLT) 

A (2) =4 MINI (C (2) ,5.5) 

C 

C COMPUTE VEHICLE C°ERflTIK*G COfiO IT lO*'’ S AT “YAJOP TI'-*5 PCTMT 

C 

IFLArrn 

60 CALL VFHTC (V (3) , A ( 1) ,TdCCT , Y, JG'^AP, JG-AP2) 

NT‘^ = MT‘^ + 1 

IF (p»-'C,GE.O, )P:-‘CT (KTM) = P‘'C/E6U^ + PINNL0 
IF(PGO.NF. 0. nPVClT(NTH) = FG0*5ML!HvOTMNLn 
TPdFLAG.FC, nCC TO o? 

CALL GF.SHFT (TFLAH, JGEAR, JGEAPP.f-CEAF) 

TFCIFLAC.'^n.nC-'' TO 80 
93 RPI =YCPT d 6 ) ^ 1 m 0 . 

A3PT=A=S(oPI) 


C-41 


GOLETfl rCPT=fN 1.2 • TZVI-fiUTC ’FL “ (Tl-10-72) 


HYORTC2 


000267 
' 1270 
J271 
000310 
000312 
000316 
000320 
000321 
00032U 
000324 
000326 
000320 


000332 
000334 
000335 
000337 
000340 
000343 
000345 
000346 
0 00350 
000351 


b >i0 451 
00 0^*54 
0 0 0‘*5‘= 
000457 

oont*ei 
noo«'63 
oon«*66 
n oo«»7o 

000i*72 
000502 
000504 
000505 
00C511 
00051t* 
000517 
000521 
000522 
000534 
000537 
000541 
000543 
000546 
000550 
r 553 
L 555 
000557 
000561 
000564 


no cu j = i , ?n ,2 

JTH = J 

IF (AOPI.GT. 2*J-t,4NT,flFOI.LE.2*J)G0 TO 9*= 

94 CCNTINUF 

95 IFfPFI+2.0,GT.0.n)GC TC 96 
T8P(JTH+l)=ieP(JTHfn ♦! 

GO TC 97 ■ 

96 ieP(JTH)=I5 F (JTH)+i 

97 CONTI^JUE 

DO 95 JV=1,T9 
EK ( JV) =YOCT ( JV) 

96 CCNTIMUE 
C 

C Tlf-E TC FPTNT FE9ULT5 - 

C 

IF (T.EO. TF) GC TG mO 

(K .LT. NFFVT) Gr TO 3 10 
K = 0 

lOQ EK IN = V^'ASS2^VGFE**2/?. 

EB = Y (18) 

FB?=Y(19) 

JG=JGFA9 
JG2 = JGFA'^2 

PRINT 710 »T . V (1 1 ) , Y ( 12) , >^RW . FBRK ,Y ( 1) ,Y ( 4) tEKIN , 
lJG,RFMTClt7<2) ,Y(5) tY(7) ,Y(9),Y(13) ,Y(1P) .PSC.PECtOVC.PGC, 
?JG2 ,fp'^TC?,v(3) ,y(F, ) ,V(P) ,Y(10),Y(19) .PSC^.PECZtPf-lCP 
I»^ (T.LT.tF) GC Tr iin 
ri ST (Jf'YC) = Y ( 1?) 

CO 101 J = 1 , ^ TG 

PMCT ( J) = P»*r T ( J) ’ 1 m 0 , 

101 CONTIN'Jf^ 

CALL TE=(^TV,n,r,nT’^LTP.P^OT) 

DC 102 J = 1,KT)’ 

TF (PKOT( J) .LT.O.DGC TO 105 

ERG ( J) =I NTl (F'^CT ( J) / wc , PCENS.ECFNS , MQENS ) 

ERG r J) =EPG ( J) *WP 
GC TO 102 

10' FPG ( J) =WP*FFVAY/F'JUoG 
10 2 CONTTNMJF 

^ DEPL ( JCYC) =2. •D^'OT ( I ) /'pn ( 1 ) 

N' T N! S = N T « - 1 

CC 1 06 J=1 .NTGS 

C£PL(JCYC)=rFFL(JCYr) + (3.f(-l. ) ^» -^ ( J-1 ) ) ♦ P*^CT ( 1 ) /ERG (J + 3) 

306 CONTINUE 

CEFL (JCYC) =CEFl (JCY r)/(3.*36C0.) 
CEPL(JCYC)=CEFL(jrYC)/r)TST(JCYC) 

VA'/G (.JCYr)=ciET(JCY'~}*3 6 00./7P 
FCCNS(JCYD=Y(ie)/Y (12) 
lCGMS( JC YD = re/ ( 3. 6' Y ( 1 2 ) ) 

FCCNS2(JCYC)=Y(17)/V(17) 

ECGNS?(JCYG)=r3?/(3.f>Y(l?)) 

ECpE(JCVD=Y{P) / Y(12) 

EC':yS(JCYC)=(7(')+Y(P))/'^(1?) 

PRTfslT ot-r|,v(i?) ,=TCN'^(JCVC) .'CCN'pCJCYP) ,"CC'’GUCYD.ECC?'?(Jf" 


C-42 


GCLETfl 


000601 

'0615 

0617 

000671 

000623 

000660 
000662 
0 00 66 3 
000665 
000666 


000670 
000672 
000673 
000676 
000677 
000702 
000704 
000706' 
000715 
000720 
000722 
000731 
1732 
■000743 
C00745 
000750 
000752 
000754 
00C761 
000765 
000 766 
000777 

noiool 

CO 1002 
001003 
001006 
C01007 
001010 
COlOll 
001012 


001013 
001014 
f' '015 
016 
I’ll 017 
001020 
m 1 021 


PCRTCf.^ 1,3 • 5E»-T-CL‘TC • (Cl-10-73) 


HY"R TD2 


C 

c 

c 


r 


= 9INT 94 5,V:VC(JCYr),n‘^CL(.JCYC),Y(lt*),v(lc5 
GO 101 J = l,l<? 

IT ieF=ITI=c+iEP ( j) 

loe CONTINUE 

PRINT 020,IPC(161 ,I?C(16) ,I?P(14) ,I9P(12) ,I9°(10) ,I!3°(M»I5F( 
1I0P(4),IBF(2)»(I“P(J)»J=1,19,2),ITIPP 
00 1 09 J=l,ic 
I0P( J) =0 
105 CONTINUE 
ITIEP=C 

IF (T.r,E, TF) GC TC 190 


110 

121 

122 
13 0 

IF 0 
170 

IRO 

IPO 


INTEGPaTE TC NEXT TINE STEF 
TTMF=T 


pmu, I , 
POUR QUauty 


GO 170 JV=l,-3 
JA=M/ixo ( 2, JV) 

IF ( JV. EC. 2) GO TC 121 
TTMP=TTMC+CELT/? . 
ot=ttmp-t 

IF { JV. NE.3) GC TC 122 

a(3)={V(6)-(VTMF+VTHPL)/?.)/C3.e»OELT) 

A(31 =AMINl(fini ,9.6) 

VT.YPL= VTMF 

VTNP=Y (lI)+EK{llf(JV-n*19)’CT 
TFLAGrO 

CALL VFHIC(VTHP,A{Jt),TTf'P,YrCT,v,jGEA=',JGE4C2) 

TP (TFLaG.FO.l)GO TO 150 

CALL G PS H PT (I FL /'•■•;, JGEAP . JG F A p 2 , NGE A P ) 

IF (IFLAG. EC. T ).GC TO 130 

CC 170 JK'=l,lc 

FK { JKf 1C*JV ) =YDCT ( Jt< ) 

CC^'TIMUf 
OC leo JK= 1 , 1 E 

Y(JK)=Y(JV')4CELT»(FY(JK)+2,^(EV(J<flO)4.ri^(jK,.:5o. )),.rX(J<4-57))/6 

ccntismif 

T=Tf CELT 
GO TC 50 
CCNTINUF 
FCAV=0. 0 
FCAV=0. 0 
FHEa\'=C. 0 
fSYSA7=0.0 
CO 300 JC=l,.\Crvc 


C CO“tPUTF FL'PL ant FN^pCY CCN'SUNPTION CN CCVFCSTTE d°iving CYCLF 

C FOR ^OCFSl ASC 2, 

r 


EGPAP( JC)=0.C 
FceAR{jc)=n.c 
FCCAP2 ( JC) = C. 0 
EC9AF( JC) =0.0 
' C5AF2 (JC) =0.0 
FHERARf JO = C.r 
GSYSPP ( JC) =C, n 


GOLETa FCP^P^^ 1.3 EE“T-AL'7r 'TPL 


('’1-10-73) 


HV<3P TO? 


001022 
nni023 
025 
bu 1 0 33 
001040 
001045 
0J1053 
001050 
001066 
001073 
001101 
001103 
001106 

C 
r 
- c 

001112 


c 

c 

c 

0 01 120 "' 

r 

c 

c 

c 

' ' 123 

.125 
001126 
001126 
001130 
001134 
001135 
00 1135 . 
001143 
001145 
001146 
001154 

001170 

001171 

001175 

C 

C 

C 

C 

001175 

001203 

C 

C 

C 

210 

0ul214 

001217 

001223 


vrec:p = o. 0 

CO 260 KC=1,VCYCLE 

VR£CIP = VPFCIP + G t«“A(JC,KD/VflVG('<C) 

5O0AR(JC) = PC?/SP(JC)*Ga"Ka(JC,'<C)*OEFL(KC) 

FC8ap(jc) = FC5aR(JC)+r-aFKa(jc,'<C)*‘^C0KS(‘<c) 

FCBAR2(JC)=FCeaP2(JC)+Ga^“a(JC,xr) *FCONS 2 (KC ) 

ECBAF( JC)=5C6AP( JC)+GflPMA{JC,KC) *5CONS(KC) 

ECBAF2 (JC) =ECPflP? (JC ) ♦GAHva ( JC ,KC) •cC0MS2(XC) 

PHEBAR ( JC) =EHEPAP ( JC) +GA»^«C { JC ,KC) *FCH5( XC) 
ESYSER(JC)=FSYSEP(JC)+GAK^A(JC,KC) EC SYS (XC) 

260 CONTINUE 

V8AR ( JC) =1 . /VFECIP 

SPPWR{JC)=lCan,*ECOAF(JC)*VBAR(JC)/W3 

CCRPECT FUEL CCNSU»^FTICN FCR NCN-7ERO ENERGY CONSL'YPTION ON MC 

FCeAR2(JC)=(FCRAR?(JC)*FCBAR(JC)-FCPAR(JC)*ECBAF?(JC))/{PCEaF( 
1-EC3AR2 ( JC) > 

CC^CUTE RANGE FCR NEW BATTERY DISCHARGE LIYIT 
RANGE(JC)=CF»^AX/POBaR(JC) 

CCNFUTF FCACTICN C'' total CFIVP'G done on NOOE 1 FCR EACH 
CC‘^FCSI'''E CYCLE. 

IF (JC. NF. 1 ) GC TP 
C L C w = T . 

GC TO 

265 CL CW=QSUF ( JC-1 ) 

266 1*^ (CHAom.LF.CSI.'F (JC) )GC ’’C 2F7 
FFRAC( JC) =0 . 

GC TC ?sn 

2F7 IF (RANGE (JC) .GT.DLCW-CWAFM)GC TC 270 
RFRAC(JC)=RANGE(JC)/C=AP(JC) 

GC TC 2fl0 

270 IF(=»ANGF ( JC) .GT.CSUF ( JC) -OWAOMTGC TO 2^5 

KFRAC(JC)=(RANGP’{JC5/C°AC(JC))’('.-(RA>'GE(JC)-OLCW + CwaRH)^ + e/ 

1 ( 2 , »RA NGe ( JC ) ♦ ( CSU= ( JG ) - CLOW ) ) ) 

GO TC 2?a 

27 5 FFRAC(JC)=l.-C'!A=w/pQ;io(jr) 

2f0 CONTINUE 

CCYOUTF MCOE-AVFRAGFD ^UEL AND FNEcr,Y C0^SU^' = TTC^ FOP CChPCSTT 
0RI7ING CYCLF. 

PCMAV(JC)=RFKAC(JC)*FC=AP(JDf (1,-PFP.aC(J''))*PCOAPP(JC) 
ECNAV(jr)=RFFAC(JC)*PGRAF(JC)+{l.-RFPAP(JC))*EC?flR2(JC) 

COMPUTE OVERALL FUEL ANC ENERG,Y rONSUHPT ICN'. 

FCAv=FnAV4-nNC(JC)*='CHA'/(jn) 

EC A7=ECA V + <^NC ( JC ) ( JO) 

PMc^A\/==’prf. V4PNC( JC) 4CPPPAT (JC) 

E3 YS AV=r RV C C V ♦ RNC (JO’F’EYRRCfjr) 




GCLETA PCKTPAK 1,3 * SEHl-CfTC PP’L • 


H Y 3 P 10 ’ 


001226 
' • 231 

. *233 

00123E 


001237 

0012^*1 

001246 

0012E1 

001302 
00131i* 
001324 
001332 
00 1335 


00 1337 
00 1337 
001337 


L u 1 3 3 7 
001337 
001337 
00 1 337 


001 337 


001337 

001337 


001337 


001337 


001337 

001337 


L - 337 
00 1 337 


FE/SV = 1000.»FL'ELEG/FCav 

WPflV=ECAV/C>'GE‘'F 

300 CONTINUE 
C 

HEEFrPHE AV/ESTSfl '/ 

c estikate fattery life. 

C - 

C0AV=05HAX 

2LIFE=CSTA V^TMl ( CCA V , OC ISCH ,C YCL E S , NOOS CH ) 

PRINT ccQ 

PRINT 960t <FCEAP ( J) ,FCRAF2(J),EC"AP( J) ,RANGE( J) ,FCHAV(J) ,ECHAV( 

1VBAR(J),SFcvR(J),J=1,NCFwo) 

CRINT 970,FCAV,ECAV, PLTFE 
PRINT 975 ,FEAV, WOAV 
PRINT 1010,VEEF 

194 CONTINUE Oj;iCI.\'AL PAGE IS 

2P0 STOP OF Purdl QUALITY 

C FCRRAT STATEf'FNTS 

C 

600 FORMAT (I10,6F10, 4) 

610 FORMAT (TFIO.^*) 

700 F0RMAT(lH0,4X,i^HVI‘'F,5X,5M5FEEC,5X,P.HDISTANCE,4Y,9WO.L.FCWEP,3Y 
16HBR.POWFF, 3X»10HR. L.FNEFGY,3X, 9Hr3P, ENERGY, 2X, 10 HK IN. ENERGY/ 
21X,2HGN,2X,3HFn»',lX,10HC.T.ENFRGV,?X,10HSYE.PNcPGY,?X,l0HFNG.EN 
3V,4X, I 

3CTCP SHAFT ENEPG'',^iy,10HPAT.ENEPGY,3X,^HGYE.PCWEF,7y,o«ENn.PCWE 
i^5X , 17HMCTCP Ef-A = T PCVr = ) 

71C PORMAT(lHr,2El?.A/T?,cF.0,10FlE.:./TE,FF.n,aE12.i*,17X,4El2.4) 

P20 FCR*-'AT(1H ,2E1E.6) 

92E FORMAT ( IHl , 27HCCNT90L PARAMETERS POP CaSE,I?) 

930 '=‘09MAT flHn,5X,?nHMTN, ENGINE TCPOUE =,E11,4,9H (‘'nCE 1) 

3 /2SX,Ell.i-,9H (VOCE 2) /6 v , 2 o rd c tT . DTSCH. LT'^tt =,^11. 4 
2/5X, 21HNOMTNAL '^OTOP OCWRR =,Ftl.U 
3/6X, 1<?HTPANSITTCN SPEEC =,E11.4> 

950 FORMAT (IHP ,A7HFUEL AND ENERGY CCNSUMPTTCN CN CONFOS’rTE CYCLES/' 
lt3HFUEL - '^CPE 1,3X,13HFLEL - Mpr-E ? , 2X , 15FENERG Y - '^OOE 1,1X, 
P34HPANr,E - MCGE l,J*X,nHFUEL - A V, , 9 Y , 1.? HENpPGY - AV. ,6X,9MAV. 
3F0 ,FX, 13HSF. Pnwrp - \) 

9F0 PCR“AT(1H ,eElE,9) 

970 FORM AT (1 KO , 29FYE ARLY AV. FUFL CC''’SUMP TI 0 N = , E 15 . c / 1 X , 3 Y F A RL Y 
.1 BATTERv energy OUTPUT = ,E 16 . 5 /I y , 7T|mEXo fp tec RAttfRY LTFF =»E! 
2) 

975 FORM AT { IHO , PPHYG APL Y AV. FUEL cCCNON'Y --.'le.S/lY, 

129HYEA°LY Ay. WALL FLLG CUTFU^ =,E15.5' 
ci;Q fodmaT (IHO, IFhCYCLE CISTANCE = , F t 2 . Al / 1 x , 2 5PFUEL OC^'SUXPTIQN ON 
IF 1 =,F12. t*/?! X , «H*"CCE 2 = ,E 12 . /I X , 3FHR A TT . E'^'ERGY CCNSLMFTICI 
2 MODE 1 =.E32.i*/2PX, ?H“GG£ 2 =,E1?,t^) 

1 03 0 F09*'AT (IHO ,29FHFAT F^lrxvF F^;^9r.Y FPACTinN =,916.*=') 

9^*5 FOFMATdH .lEHAVERAGE SPEED = , E 1 2 . , 1 X , 1 RMF a TTE R Y DEPLETION =,‘ 

16 » IX ,21HHEAT ENGINE ON TImf = , E 1 f . 5 , ( M OOE 1 ) , 1 X , f le . 6 , ? M ( ►"CD i 
2) 


ICOQ FORM AT ( IHO, ip 1 ?. J*) 

1 020 FOP*"AT (1 NO , ♦ °A TTEpv c^wEP D T ST p T "U T I n m • , 
1 /IX, •35-Ti-,lX,*3i-3Q',3X,*30-25’ 


/50X 

, 1X» 
, 3 X , 


, *nf dative 

•26- , IX 

’ 2- 7 * , IX , • 


• t IX , "FC 

, ‘ P ’- l 5 » 


2 - 6 * 


7 Y , 


» 


GCLET4 FCF^Ff^ 1 , • EE‘'I-M;TT! oc ^ « fni-!0-7T) 


UY=IP I 


1337 


^’3 0-34*,iy,»7t-.’«*,iy,*TrTflL*/U20IF<)) 
ENQ 


oi;ii;r;AL page is 

QUALITY 


C-A6 


GCL‘‘TA FCF'^'■^^ 1.1 • SfMT-«UTC ’‘='L * (Cl-n-7?) 


000012 


000012 


'’0 00 12 


70 00 12 


000012 
7 0 0 0 1 e 
? noo 2 o 
700022 
000025 
700027 
70 0031 
'’00033 
700036 
0000 A 2 


700046 

700051 

700053 


70 0 0«=4 
70 0061 


SUS'^CUTINF V^HIC ( V. A ,T, YTCT , V, vnTfi p?) 

C 

C VEHIC rOMpijTES Fowpo PFCL'IRE HE MF AND LC'^SES F° C *• FEAR WHEELS 

C THROUGH EKGTME 

C 

C CflTA IN CCVHCK - 

c - - 

CCHHCN OISPL,PPH(20).TnwAX(2C),B^'EP(20),BSFC(2n,20),N'RPM,N’H!EP. 
lEKUH,EHUG,PINMr,SPHO(201 ,THCM£X(2 0) ,T*^CHIK(2 0) .KSFHO, 

2CTC7P0 ,TSP(2C),T'3R(20),SFR(?C),TSP2(2n).H’TSP» 
3NLCCK(5),TRATI0(5)»CTl(5),CT2(5),EHUT(5)t 
t^wo ,EEHAX.EHijPG,P“UPG2. 

FCR ATIO ,CD1 ,Cr2, E^UO, PTTPE .CTTRFl ,C TIP'"?, VvaSF.TLI.fGA, 
f TEOHIN,rPHax,PPHTDL,'yvASS2 fPEC^If-. , 

Tf^UPL SG .SKALEi T^'S^L ,NTW, NGGSCW .GHC-F FP . VCLE , NC C H© , 7S ’’ A V 
6NUNITS,THSET(51 tUPSHFKE) ,nr)ISCH(?0) , C YC LES ( 0 7 , MG { 

9NPRTC(3) ,CTC (3) ,TFC (3) ,CSU° (30) ,GNC (30) , C5AP(30) 

1 . RPMEST , FPFST , FE 0 P ( 2 C ). RFHECF ( 2 C ’, PPY 0 OM , EYUCVT , PATUP , PA ' fnN,NC 

2 . PHFNS (20) f SCENE (20) .NCENS 

3 . VHAX, TEOHNP ,UPSEL tCNSEL .CWARH ,NAX , Rcuay (jq) ,TAX(20) 

CCHHON TIHC (3,200) ,SPECC f3,200 ), GAMMA (30,3) 

C 

C VARIABLES IH CCHMCN - 

r 

COMMCN FA,FR,FAC,fnET,TCC,7TC,RP“00,PPYTC,TSC,RPMSC,P'’MTC0,TTC 
1PP‘'TCN,Tmc,TFC,T‘'P2, TEC? ,TGC,PSC,PEO,P‘^C? ,PMC,OMC2,PGO,P5PK ,PC 
?‘=0,PT,°pw,ci.tC,TH’5C'',TTC7N 

3, P°*^Fn,PDvirrT^=LTC?,V^'=^ 

4 , TSC 2 ,'-' F ‘' EC 2 , FSC 2 , pp ‘^ TC 1 , RF ''‘ TC 2 ,'^ R‘-'’,^‘^''E = ’ 

c 

REAL YCnT (30) ,TV’T1,T"'T9 

c 

C CO'^PUTE SPRECS AND lC A C - tnCF =E ^ C - N T LOSSES T‘-''50UGH C°IVE TCAIN 

C 

VMPS=AMAX1 ( 0, 0,V/3. 6) 

NLC'<=NLrCK (NGFAP ) 

NLCK'2=HLCCW (HGFAR2) 

<=cM0C=9.54E2Ee7'»VMF'^/PTl0r 

ooMTC-RPMCC*rFATIO 

oDM3C = '5PMTC>fTCATTn (NGEA-R ) 

PP,MEC2 = RPMTC*TPATIC(NGEAP?) 

TLF7 = CC14-CC2"FF’'T0 

TLFT=CT1 (K'GEAP) + CT? (N'GEAC) *ReM.Sr 
TLFT2=CT1(NGCAP?)*CT2(NGFAR2)*RF“S02 
C COMPUTE PCAC LOAD OOW'^P PECUIREM'NTS 

C AEPCCYNAMTC - 

C 

FA=n,6125’CCA*VVPS’^»? 
f^fa'^vnfe/ 1 ncju ■ -- 

r YjoT (1 ) = pa'7Tdoo 
C " 

C POLLING PESIETAHGF - 

c 

‘^R = 9. 807*VM£ES*(rTIPFl+r7IRF?*V) 

PR = FP*V“ PS/ICC i7. 


/CLETA FCPTCfl^ 1,3 • EF^T-AtlTC =*^1 


( n 1 - I (1 - 7 ■» ) 


., -u ir 


OOOOEl 


noooes 

000067 
000073 
000075 
00 0077 


000100 


000 103 
000110 
000112 
000120 


000122 

0134 

•JQ0137 


000143 
000147 
000153 
000155 
000162 
000165 
000166 
000 167 
000172 


000177 
000201 
000202 
000206 
000212 
0 0 n ? 1 ? 
000215 


000221 
C00222 
0 0 0 ? 2 3 




YCCT (1 1= tFA ) /lOOO. 

C 

C eCCELEf^ATTCh - 

C 

FAC = VMASS2'» A . ' 

FNET = FAf FFtFAC + FG * . i j . , ' 

TOC=FNET*RTIFF ' . . * 

PF’HsO. 0 0 010 4719 7^TOC-»RFyCO 
TFIFKET.LT. 0.0)GO TO 200 
C 

G NON-NEG. NET CRTVIKG FCFCE TC THE WHEELS IS NEECEG 

C DIFFERENTIAL - 

C 

TTD=TOC/ (EMLD’DPATTC) ♦ TL^D 
C 

C TRANSMISSICK - 

C 

TSC=TTn/ (EHL'KNGEAPl »TRATIC(NGEA9)) +TLFT 
FSC=TSC»RFMSC*0. 00010471 67 

TS02=TTC/<cJ'UT(NGEAF?)*T = ATIC(NG'AR2)) *TLFT2 
PS02 = TSC2*RF>'Sn2*C. 0 0CIC47197 
C 

C IS CAR AT REST AKC ICLING - 

C 

IF (ABS(VNRS) .LT. 0. 2. AND. AES (A) ,LT, 0, DGO TC 300 
C 

CALL PMCVPd ;jeF><,NLCK,NLCK?) 

11-C GO TC AGO 

c 

C NET CRIVING FCRC? AT WHEELS IS NEGATIVE 

C 

20 0 TTO=TOn*EML'C/CRATTO*TL9r 

TSC = TTG*EHU7(NGEAR)/TFATTC(KGEAC)i-TL'=^T 
PS C=TSn'^RP“SC*0, 00010 47 197 

TS02=TT0’‘FHL'T(NGEAC?)/TRATIC(NGEAR?) +TLFT? 

PSC2=TS0 2''Cfncc 3*0, nnni C47 1<?7 

MGEAP=NGFAC2 

IF (PSO.GE. C.O)GC TO 210 

CALL P‘"CVPr-l»JPP'<,NLC'6,»LCK2) 

IFU8?‘<,En. OIGP TO 400 

c 

c 

PSO=T3C’RFHSC^0.0n010471P7 

'=S02=°SC 

TT0=(TS0-TL‘^T)*ToaTIC(SGEAR)/E^iJ7(NGEAP1 
TCC= {TTC-tlFO *rRATTr/«^HL": 

GC TC ago 

21 C CALL PVCVP (1 ♦J°F'<,NLC*<,» LC<2) 

PLC GO TP AGO 

C 

C CAR IS AT REST ANC IDLING 

C 

300 A =0.0 

VHCS=0 . 0 

CALL P“CWC ( n . J°c< , MCK i 


C-48 


GCLETA 


000225 

0226 


000227 

000234 

000235 

000235 

000240 

000247 

000251 

000256 

000257 

000257 

000262 

000271 

000273 

000276 

000300 

000302 

000320 

000333 

000335 

000337 

"00340 

3343 

U')C344 

000345 

000347 

000351 

000353 

000354 

000355 

000356 

000360 

000361 

000362 

000364 

000356 

000370 

000377 

000406 

000407 


FCRTPAN 1,3 • SfVl-fiUTO RFL • (01-10-7?) 


VEMir 


C 

C 

C 


FSC=o. n 
PSO2 = 0. 0 


Pour quality 


COMPUTE EVG3VE PUEL FATE flS'C CEPTVATIVES CF VAFTAPLES 


400 IF (PEC.NE, 0. one TC 410 
■■ Fc=n.a - ■ • • 

GO TO 4?0 

410 RRf^EP=125.664'»TF0/niS = L 

SFC=INT2 (RF^EC,°R^'EP^R = H ES'='C , r'RPV, KMEP) 

FC=SFC*PEC 

420 IF (PE02.Nt. 0. C) GO TC 470 
FC2=0. 0 
GO TO 4tn 

430 3RVEP2=125. 6e4*TEn2/niSFL 

SFC=I)!T2 (RPt-ECE ,PRKEP2» R F?^,B’^EFf PSP■C,KRP^",^^'EP) 

- PC2=SFC*PEC2 

440 PC=Tro*Rpvnc*c. 0001047197 
PP.RK = P9W-cn 

P7=TTO'»PPPTC*C.C001C4?107 

Y2nT(?)=(A8E(PT-OD)+AR9(FSO-PT)+A3S(PLTC))/iaOO. 

YDCT n)=(A5E(FT-ao)fA5S(FSC2-PT)»A9S(PLTC2)l/10C0. 

YDCT(4)=ApS(PeRVr)/10D0. 

YOCT (5)=PFO/1COC. 

YOOT (6) = PSC2/10Cn. 

YDCT (7) =PEO/100G, 

YHOT (S)=nF03/inro, 

YOCT (9) = P''C/in0P, 

YDCT (1 0) =PYC2/lCan. 

VCOT (11) =A *7.6 
YDCT (1?)=V’^FF/^ non. 

YDOT (13) =FGC/10Cn. 

YCDT(14)=C. 

IF (PEO.EG. 0. )GO TC 44? 

YOOT(14)=i, 

442 YDOT(15)=0. 

IP (PE02.EC, 0. )GC TO 444 
YDCT (IE) =1 . 

44L YDCT (16) =Fc/3eno, 

YDCT (17) =FC?/76C0. 

YCOT (la) = ( = I»'MC + P^C/E^(UY*FGC*FY'„'G*ECU°G ) /lOOO. 
YOCT (19) = (oiKMC+P*^C2/5^L*' + FGD*E'^UC-*FMi)pr-2) / IQOO. 
RE TUPS 
ENO 


C-49 


GCLETft 


000003 


000003 


000003 


000003 

000021 

onco‘?7 

00 0071 
000106 
000110 

0001 1 1 
0 00 1 1 3 
0 0 0 1 1 E 
000117 
000121 
000123 
00012E 
00013E 
0001«<1 
00Qli»U 
000 157 
0 0 0 17 6 
000202 
000211 


000226 

000250 


prcTr»K 1.3 • ^FL • (01-10-77) 


SU!7?r'JTINF (NP-FfiP ) 

C 

C HY9FA0 I)’PUT5 O^TA FCP HVoqiCZ 

C 

C CATA IK - 

CCMMCN DISFL,FF»'(20)iTCyAX(20),P‘'EP(20),eSFC(20,2n),KRP(<»N)'EF, 
lF“'JM,E)"UG,PTNKLr»SOK'0(2C) ,THCMAX(?D) ,TMOMIK(2 0) ,MSFMC, 
2CTCZPO,TSP(2n),TnP(?0),EFR(20),T5P2(20), NTSP, 

3N'LCCK(5) ♦! PATIO (5) , C T 1 ( 5 ) , C T2 ( 5 ) , EHUT (5 ) , 
fcWR,EEMAX,EHL'PG»FK'UPG2» 

5nRATin,C01 ,CC2tEMU0,PTIPE,CTIREl,CTTqf:2, VMASS , DL I , CCA , 
eieCPIS' ,CB»"aX iRPHTTL, V»-AEE2 tPECKIK, 

7FUELSG tSKALEiTMCKL ,KTH,VCCSCH ,CPGE FF , NC YCLE , NC C’"'’ , CHTA V 
EK'UNITS ♦THEFT (5 ) ,'JPSWFT ( 5 ) ,CNSHFT (5 ) , OniSCH (20 ) » C YCLE5 ( 20r , NTC ( 
PNPRT C( 3) ,070 (3) , Tc-C (3) ♦CEL'P (30 ) trsC (30),C5AP(30) 
ltRPHEST,PEEST,PEOP(20),PF“ECF(?0) , PPNOPM , E^UCVT , PA TUP , R A TO N , NO 

2, cnEHS (20) ♦EOENS (20) ♦NDENS 

3, VHAX,TEOHN2tUPEEL,CNSEL,DWAR“,NAX,oPH4X(20)»TAX(20) 

CCHHCN TIKC(3,200),EPErr(3,2CO)t'‘-AMHa(3o»3) 

C 

C VARIABLES IK CCKVCN - 

C OHM CM FA ,FK i FAC , FNET,T0C, TTC, PPYQO , RPf^T CfTSCtRFYSC ,P°^TCOt TTC 
lPPHTCN,THC,TECtTMn2tTEC2tTGC SC ♦ pE( 7 ♦ °EO 2 ,FHC ♦ P NC ? tPOC , PERK,FC 
2FC»PT,FPW,PLTCtTHpor,TTCIN 
3 ♦PPMEO ,RFf'EC2tPlTC2» WFS 

«‘tTSC2,RPMSC2fFSC2,PPNTri,RFKTC?t°RMEPtOPyEP2 

r 

C INPUT ENGINE CATA 

C 

READ 605 tNPFHtNHEP»FUFLEC,CISPLt'^<ALE 

READ &in,(RFV(J),TOKAX(J)»J = l,NPF'1)t (pM£P(J) ♦J=1iN»'EP) , ( (eS^Ct 
It J = l,N»^f P) ,I=ltVRP-) 

PEAD 600,KCFtPP‘'PST,PPFST 

PEAC 610, (FFCP(J),PP»'ECP(J) ,J=l,^'OP) 

cispl=cisfl*skale 
P8EST=PPEFT* exale 
rC 11 J = 1,^’FPV 
TOMAy(J)=EKALE''’0‘^AX(J) 

11 CCNTTNUE 

OC I? J=1,NCP 
PECP(J)=PECP(J)"EXALF 

12 CCNTTNUE 

PRINT 600 ,PLEL^G tPTPPL 
PRINT 805 
CC 1 J=1,NRPN 

PRINT 310,PF»'(J)tTO*'AV(J), (l),8SFC(J,l) 

PRINT 806t(5NFP(X) ♦°EFC(J,K) ,K = 2t‘(MFP> 

1 CONTINUE 

print ')07,P.PVPET,PRF‘‘T 

PRINT giO, (PECP(J),RPf'F.CP(J) ,J=l,NOP) 

C 

C INPUT >^CTCP /CENFP A TCP CATA 

C 

READ 6 !^0 ,NCF «C ♦ F‘"UK' ♦ F^UC , T NS XL ♦ P ! NN L 0 ♦ '’p x I CL ♦ ® ° " CP 
PEAC 6 1 0 , ('^F^'C ( J ) ♦ T*' rvA X ( J ) , 'rvcNTM ; J) , J= 1 , KSPVO ) 


GCLETA 


000267 
0272 
. -0274 
000276 
000300 
000303 
000321 
000325 

000344 


000352 

000354 

000366 

000400 


000402 
000410 
00 0431 
0Q0435 
00 0443 
000447 
'*'’0466 
0 «i71 
000477 


000502 
000516 
00 0547 
00 0566 
CQ0572 
00 0576 
00 J£27 
nn 06 37 
000643 


000662 
000704 
000721 
000736 
0 0 0742 
00 0746 
00 0 764 
1770 
O- 1005 
001011 


FC'JTFAK 1.^ • 5E“T-£L'Tn CFL • (CJ-lC-73) 


C-50 
uYcp •jn 


C 

C 

c 

c 


c 

c 


c 

c 

c 


c 

c 

c 


c 

c 


CO 2 J = 1,^•5F^■': 

TMCHAX(J)-Tf'CrA7(J)*T“5'<L 

THCMIN(J)=T^CPIK(J) 

2 COK’TIMUE ' ; 

PRINT 511 ill 

PRINT 812»PKUP,F*"UG, FTNNLntPPMTOL, RPMOOM 

PRINT »13 

PRINT 514,(SPRC(J),T>'CHAy(J),T(-CNIN(J),J=l,NPc>^C) 

REaO 600,JCVT 

JC'/T = 0 FPF TCSCUE CCNVEPTCT, =1 ^OR HVT 

IF ( JCVT, Er. C) GC TC 4 
READ 6in,F‘a'CV7,qATUP,RAT0N 
PRINT 305,E^UCVT,PATUF,PPTCN 
GC TC 5 

INPUT TCRC'JE CCNVPRTCR CATA 

4 READ 620»ntEP 

READ 610.CTC7'^C»(TSF(J) ,TQR(J) ,SpR(J}»J=1»NT5P) 

PRINT 860 
PRINT 865»CTCZRC 
PPINT 870 

PRINT 814,(TFF(J),T0P(J) , S PR ( J ) » J= I » N T^P ) 

CO 3 J=1,NTSP 

TSP?(,n=T5°(J)'5nRT(Tr''5(j)) . 

3 CON'^INUF 

INPUT TP A ^'S^'I5STf:M CATA 

5 READ 6 05 ,NGEAF ,NT^,UP5FL ,CNSEL 

READ 601,{NLCCi<(J 1»TPATTC(J) .CTKJ) ,CT?(J) ,E''UT(J) »J=l 
READ 610, f THEET ( J ) , UPSHF t ( j ) , D nRufj { j > , j = i ) 

PRINT 840 

PRINT 8 4 I 

PRINT 842, (MCCX(J),TRATTC{J),CT1(.J) ,CT 2 {n,F^*UT(J),J = 
PRINT 844 ,UFSFL , CNSFL 
PRINT 543 

PRINT 814,fTV’EET(J),UPEHPT(J),rK«^MC-T(j),j = i,NTH> 

jmcijt p^ijxrKY DS’A 

READ 605 ,T'OrECH,NDFNE,W?,EE^'AY,E!-URG,F“UPG2,CHGEPP 
REAP 610,(DCI5CW(J) ,C7rLFS<J),J=l,NC0SCH) 

PE AD 610,(PCEKRfJ),FCENE{J),J=l,NOENS) 

PRINT 850 
F P I ► ' T 2 5 5 

PRINT 815 , WE ,EPNAX,E NUPG ,ENUFG2,CHGEPP 
FPI^'T 050 

PRINT qn,(CCISCM(J),rYCLEE(J),J = l,l'ODSCH) 

PRINT 025 

PRINT Din, (FPpNS ( J) , ECE^ E < J) , J = 1 ,NDpn$) 

TNPUT CTFPE'^Fh TTAL ANC TIRE CATfi 


, N CFAC) 


I ,NGEAR) 


C-51 


GCLETA 


' t02€ 
x050 
OOlOEU 


001076 

001104 

001121 

001125 


001142 

001154 

ooliee 


001172 
001177 
001203 
001206 
00 1223 
001225 
24C 


001250 
001252 
001 255 
001257 
001263 
001103 
001307 
001320 
001326 
001332 
001335 
00 1 35U 
001374 
001 ‘•00 
001405 
001407 
001412 
00141? 

001420 


0 0 1 4?o 
Tni4?n 
001420 


FCC7FA4 1.3 • 5E''I-AUT0 °FL • (01-10-73) 




5FaC 610tr5iTICtrf31,CC2,E‘'fC.R7I'^E,CTI»El,CTI5E2 
PRINT 075 

PRINT «85»0RATTC,CCl,CC2,EYUC,RTIR5,CTTPEl,nTi;?f:2 
C 

C INPUT ACC555CRY CATA 

C 

READ SCO.NAY 

REAC 610 , (PFf'AY ( J) ,TAX ( J ) ♦ J = l, NAX) 

PRINT loao . ^ J 

PRINT °10. (RpwAY (J) ,TAY ( J) ,J = 1 ,NCX) 

C 

C INPUT VEHICLE DATA 

C 

READ 61C.VM ASS.CLI.CDA 
PRINT SRO.V^'AES.DLT ,rCA 
VHASS2 = VHA55*'CLI/RTIRE''*2 
C 

C INPUT CPIVI^R CYCLE DATA 

C 

read S20,NCYCLE 

PRINT P15 

DC 45 JC=1,NDYCLE 

read 525, NTC (JC),N PR TC(JC) ,NL'NITE,0TC(JC),TFC(JC) 
NT=NTC (JC) 

pRAD p10,(TIfC(JC,J), 5PFCC(JC,J),J=1,ND 
IF ( Vl'NTTS. V'E . 1) TP 'll 
r 

C CCNV'FCT NFM TC V'TLC^‘pTFC<^ HCt'P 

C 

30 CO 40 J=1,NT 

EPEDC ( JC ,w > ;PErc ( JC , J M 1 , 5 no i 

40 CONTINUE 

41 oRINT 920 

PRINT P10,(TIFC(JC,J),SPFCC(jn,Jl,J=l,NT) 

45 CONTINUE 

PEAD 6 on ,NCC‘'P , PSTAV .CWAPY 

PRINT 990,PETAV 

PRINT 995 

’’0 320 JC^l.NCC'^P 

READ 5 in ,CELP ( JC ) , DhP ( JC ) , (GA*'-A ( JC , J) , J = l,N'CYrAE) 
PRINT ?lr , DE'JP ( Ji") ,CNC ( jr j » ( G A { jr , j) . j=i ,nCYC LE) 
320 CONTINUE 

PRINT 991,C4Apw 
DBAR (1 ) = CSUF (1) /2. 

CO 4£ J=2,KCCN= 

C9fiP(J)-=(CELF(J)+CSUP(J-l))/2. 

45 CONTINUE 
PETtJPN 

r 

c FCRNAT STATF^entS 

c 

600 ' PNAT (110 ,6F10,41 

ei?i »-ur‘";t (iio,t.pio.i<) 
f G5 FOPYAT (211 0 , E*"! C.4) 


ifJH IS 
\[,1TY 




GOLETft FCPTPA> 1,3 * SE“I-fiinG P^-L • (01-10-77) 


H Y P E /) r) 


001(.20 
P "* i «20 
‘*20 

0011*20 

0011*20 

001420 

0011*20 

0011*20 

001U20 

0uli*20 

001420 


001420 

001420 

001420 

001420 

001420 

001420 

001420 

001420 

001420 

420 

001420 

001420 
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001420 


PIO PCRHtT (7P1 0 . 4 ) 

F20 '■QPMAT (7110 ) 

fZ5 FC9MAT (3110 , 2^10.4) 

eOO PQ'^MAT (IHO, 1 IHENGTSE C A T A / 1 X , ? 3HPUEL SPECIFIC GP^VITY =,F12.4/ 
11X,14H0ISFLACEKENT =,E12.4) 

eO: FC9MAT(lH0,7Y,3HpCY,ox,llHMAX. 7nROUE»OX,4He^'EP,12Y»4HBSFC) 

?06 FORMATdH , 32 X , 2E 1 E. E ) 

807 F09MAT(lH0,2nHBEST CPEFA7ING PCIN’T - S°FFO = , E 1 2 . 4 t EX » 7 HPO WEF = 
1E1?.4//6X,5)-FCWEP,5X,10HCPT. SPEED) 
fiOe PORHAT ( IHP , l^HCVT CH A R A C T E R I S T IC S / / 1 X , 1 2 HE F F I C I E N C Y =,E12.4,4X, 
115HSPEED1JF PATIC = ♦ E 1 2. 4 , 4 x , 16 HSL 0 W DC WM '=ATIC =,E12.4) 

810 FCRHATdH ttElE.*!) 

Ell FCRI'AT (IHO, 20HfirTDP/r-EN'EFA TC° DATA) 

ei2 FORMAT (IPn, 13HEFC, (MCTCP) = , El 2 . 4 , 5 X , 1 2H E PF . ( G EN . ) =,E12.4,6X, 
121HNC LOAD INPUT POWER = , E I 2 . 4 , 6 X , I 2H I DL E SPEED = , E 1 2 , 4/ I x , 22HF 
2 OPERATING SREEO =,E12.4) 

813 '=’OR“ AT (IHO » 7X , 3HRoP,<3X , 1 IHMA X. ' ORQUE » X , 1 1 w “ I M . TCRCUF) 

814 FCRHATdH ,3E1E.6) 

81E FCR“AT(1H ,EE1E.G) 

840 FORMAT (IHO , 12HGEARROX DATA) 

841 CCRHAT (IHOt 5X ,7HL0C‘^UP- , lOX ,EHRATindOX, 22HSFIN' LOSS COEPPICIE) 
18X»11HTCPCUE RFC, ) 

042 FCR“AT(1H , 18 ,fiX ,4F1E.E) 

843 FORMAT (IHO, 4X,10HENG. P 0 VE R » E X , 1 IH IJP S HI F T R P H , 4 X , 1 3 HC C W N 2 H I F T P 
84« FORMAT (IHP, 37HSHTRT PCINTS RCP ELECTRIC OPERATION -»2E12.4) 

8E0 FCR‘^ AT dM n , 1 2H^ A TTPC Y DA'^A) 

2E'" FORhaT (IHO, 8X,4HMASS,7X, 14HENERGY OEMS! T Y , 7X , H A\/ , PEGEN. EFF. 

llXdEHMfly F5GE>. EFF. ,px ,1 3H'=ECHARGE E=’P' . ) 

8f n format (IHO , 2 IHTCRCL'F CCNVER'^CR TATA) 

88 ; format (IHO, 17Hni^»2/tT (STALL ) =,E16.E,)- 

F7C Fo:?u^,j(lpp^3^^j^^H^jn/3,7PT(TO) ,8Y,EHT0/TI,11Y,EHN'C/NI) 

875 AT ( 1 HO , 1 rPAVLE an’C ttre ^ATA ) 

»f,E Fc CH at ( 1 HP , 1 ^HC'^'^F. CATIC = , E 1 F . F / 1 Y , 24H f r x m LOSS C CE I C T F M T S 
12E 16 .6/lx* 1 3HTPFOUF fcp, r , e IE . «= /I X , 16HP OLLT NG RAHIUS =,616.6/ 
21X , 33MROLL TNG '=FSIS‘^AKCE C G*^ = F TC IE N T S =,2E16.6) 

FEO 'FORMAT (IHO , 14HVEMICLE ‘*ASS = ,E 16 . 6 . ?X ,19 HDRI VEL I NE INERTIA =,£■ 
1 ,r X, 19H (DF AG CDEF.)*At;EA =,E16.6) 

F'lO FnQMATdH ,2Elf,6) 

C15 FORMAT (1 HP , 14HDF xvxL’G CyGLFS) 

920 PCRMATdHn,6X,4H-'TYE,l?Y,= HSFEFD1 

92E ="0 RMAT dHP , 1 X , 14 McnEG TFT r FC VE R , 2 X , 1 = wsP E r j p t C F^'cqr,Y1 
08 0 format (IHO , IX d'^'HPFnxi* CF C I SC H . , 6 X , 1 OHC YCL F LIFF) 

EGO FComat ( IHO , 24HTF AVFL G J S t p t BUT I C ^ Q A T A /I X , 11 H A V . USAGE =,E12.4 
o°l FcrmaT (IHO , fharmup CISTAN'CF =*,F12.4) 

ccr FORMAT (IHO ,3X ,10 HmaX. C I T . , 4 X , 1 E HFR A CT , OF T C T A L , 1 4X , 2 1 H 0 R I V I 
lYCLE WEIGHTS) 

inOC FCR‘'ATdH0,14HtrCESSGRY LCAD/6X,-HSPEE0,llX,eHTCR0L'E) 

PNC 
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C 

. 0007 


000007- - 

c 

c 

OOOTJ07 
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f’n 0007 


000011 
000012 
00 0013 


000014 

000020 

C00024 
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r 

c 

c 


c 

c 

c 


c 

c 


SUePCIJTlMr C^^^JPTdFLAG.JGEAC.jrra;??, nGE AR) 

CATA IK CCVKCK - 

CO MM CM niSPL.FP'^120) ,TC»"AX(2C) ,P''EPf20) , 8*^ *^0(20.20) ,NRP*",MMEP. 
lEML'M.EML'G.CJh ^L^.SFWC (2 0 ) , TMCHAY (20) , TMOMIK ( 20) .NSF^O, 
?rTC2FO,TSP(20),Tr3P(2n),SFP(20),TSP2(20)d-'T?P, 
3KLCC»<(8),TPATin(5),CTl(5),CT2(5),EMUT(5), 
t*VIB .EeMAX,EMURG.EM»JRG2, 

rCRATIO.CDl ,rC2,FMU0,PTT.FEtCTIREl ,CTIPE2t VMfiSS ,OL I , CCA , 
fTEnMIM,C8MA';t,RF»'ICL, VMAE52 .PEOMIN, 

7FU£LSG,S^ALE;TMSKL.NTH,KrCSCH ,CHGE FF , KC YCLE . NC C ''P , CS T A V , 
eKUNITS,THSET(5)tUPSHFT(E),DNSHFT{5),DOISCH(2a).CYCLES(20),MTC( 
CiNPRTC(3),rTC(3) » T FC ( 3 ) , C FUP ( 3 0 ) , CMC ( 3 0 ) , CRAP (30) 

1, RPMBST,FPE£T.pEOP(201,RPM£CF(20), RP MCP« , pyUC VT . R A YUP , R A TO K , K C 

2, PCENS (20) .ECEN5 (20) .KCEKS 

3, VMAy, TEOMK2,UP'TEl tCK'SEL ,CWAPM.KAX,PP‘^AX (20) .TAX (20) 

COMMON Tire (3.200 ) .SFECC(3,200),GAMMA(30,3) 

VARIABLES CCY“CN - 

CCMMCM FA,FK,FAC.F^;cT.TC^.TTC,PF^'00,RPMTC,TS0.RP^SC,RP'^TC0.TTC 
1 RPMTCN .TMC .TEC .TH.0 2. TE0 2 .TGC .PEC .PEO , PEG 2 . pMC . pmC2 . pG C , PEP K ,FC 
2P0.PT.PRW.PLTC. thRPO .T-CIN 
3 .RPVEC.RPMEC2.PLTC2. VK=S 

(‘♦TS02.°PMEn2,FSr2.PFMTCl . R FMTC 2 ♦ pR*^E F . RP MEP 2 
REAL TMTl 


IF (OFO.GT. C, n)GC TO E 

MOCE t SHTPT CETPRmIN-0 FY mcTCR/GEK. CH AR A C t E® I S t T 

SHPTIJ° = UPEEL 
SHPirN = r/ EEL 
GO 1? 

‘'CCc I SHIFT rrTFPMIKEC EY HEAT ‘^MGIMF C H A R A C TE F I S Y TC S 

E SHFTL'o = JKYl (PFr,’’HSET,L'PSHFT.M'^H> 
SHFTC*!=IKTi(FFr,TPCR7,CKSHFT,NTH) 
ir TP (PE 02 .GT. C. P)GC Tp IE 

MCOE ? shift CPYFpvTfvCP CY I'CTCF/GE^', pH A R A C TEP T ST I CS 


C00031 

000032 

000033 


000034 

000040 

000044 

toeo 

. J072 
CO C072 
000073 
00Q07E 


SHPUF.2 = !.'°SEL 

SHFCK? = (Th;£rL 

GO TC 20 
C 

r “OOE 2 SHI^T rEYERMIMEP OY HEAT fmGINE C h A F A C TE R T S T ICS 

c 

IE SHPUP? = IN'Tl(FEr?,THSET,L!FSHFT,KTH) 

SHFnf:2 = IMTl(FEC2,THSET,rK'SPFT,N*7H) 

20 IF (PPHTCl . GE . SH=-TUP . A^n . JGEAP. LT . NGF AR) GP TO 2r 
TR ( PPHTCI . LE . SFPTCN. A^•C , JG EAP. GT , 1 ) GO TO 30 
TFLAG = n 
G P TP LI) 

JGr A== JFRAD4 1 
TFL Ar.= i 
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000075 

n'’0076 

100 

OUOIOO 

000112 

000124 

000124 

000126 

000126 

000127 

000131 

000131 

000132 


GC TC 40 

30 JGcflFs JG6f 0-1 
IFLflG=l 

40 IF{RF‘^TC2.GF.SHruP?.fiK’T. JG5AF?,LT.MGEflF) GO ^0 45 
IF (RPMTC2.LE.5H=’'?N2.fiNT. JGEiJF^.GT. 1) GO TO 50 
GO TC 100 

i‘5 JGPAP2 = JGFAR2+1 
TFLAG=1 
GO TC 100 

50 JGEAO?=JGrflF2-l 
IFLAGsl 

100 FETUPM '"AL PAGE LS 

ENP ■ ’ '*^1^ QUALHY 


G0LFT4 FCFTPA^ 1,3 * !iFHI-/ilJTC 


» n 1 - in-7 2 ) 
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'0007 


00 0007 
000007 


000007 

000007 

000010 


0011 
■j 0 0 0 1 ? 

C00013 
000013 
000015 
OOOOIS 
000020 
000024. 
000025 
000037 
■100030 
0000 3^ 
•10 0 0 i »0 
00 004.3 
00 004*4 

10 ooi.e 

000050 


000054 


0070 

.0071 

100073 

110173 


SUS^CUTINF C*'CVP (K'PPV».JBR.<,K‘LCK,NLCt<?) 

CO^’^'CN CI£FL»FP>'(20),TO^£X(20)»®‘'EP(301»°SFC(30»20)tK9PK»K“5F, 
1 EMUM ,FMUG, FINK'LC,5c*^C (3 0 ) , TMCMflX ( ?0) , T“0 I ^ ( 2 0 ) ,K'SFmo, 
?CTC7FO,TSF{20),T'10(?o),SF 5(20),TS°?(20),NTSP, 
3KLCCK(5),TC*lTICf5),CTl(5),CT2(5)»EMUT(5)» 
4*MP,FeMAX,F>'l.'FG,‘^>'UPG2» 

FCRflTI0,C01 tCC2,FMUO,P7TPEtCTIREl,CTIRE2, VMA5S,0LT»CDA, 
eTECMIN ,09PAX iFPf^ICL, VPASS2 » PECr'TN, 

7FUELSG ♦SXALEtTMSKL ,KTH,N'C0£CH ,CHG E F F , kjC YCLE t NC CHP » CSTA V , 
PKUK’ITS *THSET (5) tUPSHFT (5 ) tCNEHFT (5) *nOISCH(?0) , CYCLES (2 0r, K'TC( 
9NPRTC(3),CTC(3),TFC(3),CSUP(30) ,C(iC ( 3 0 ) , CS AF ( 30 ) 
l»RF?'P5T,FP£ST,FFCO(2n)*FF>'ECF(201 , P PMCPM , EyuC'/T , FA TUP , R A 70 N , NC 
2,P0E(vS (2 0) tECENS (20) .K’CENS 

3 ,VMAX, TEOf^N2,UP£EL,Cf.'£EL lOWAPM ,)<:CX ,RPMAX (20) ,TAX (20) 

CCM“CN TI^C(3,200),SFEDC(3,200),GAHMA(3n,3) 

CCH'^CN FA,FF,FAC,FKpT»TCC»TTC,RP‘'0 0,RPHTC,T';c,RO)‘SC»PP“TC0,TTC 
1RPmtCN,TPC,7EC,T>>‘C?,TEC2 , TGC ,p SC , PpO , 0£0 3 , =M G , °mc 3 ,PGC , PPPK ,FC 
2P0,PT,PPW,FLTC»THRPn,7TCIN 
3tRPPEG,PFHTC2,PLTG2,VKFF 

4^,TS(12»FPh'SC2»FSC2tRPHTni,RFHTC2,"RMEP,9R 
REAL INTl,ThT2 
■ CC=0. 0001047IE7 

IF (KTRV) 50,80,10 
C 

C SYSTE*^ rcw'='p CUTOUT IS K CN-KFG AT TVE 

C 

in FBPX=0,0 
Tc-n = o. c 

JpRK=n 

TTCn=T£C 

fomtCO=°F“?f 

CALL TOrOT (1,1 ,KLFK) 

TACC = TN'Tl(FF*'TCN',=FH£X,TrY,N';x) 

TTCT^•=TTCI^ ♦TACC 
RP‘<TCl=PP‘'TC^ 

ttcni=ttcin 

THE’*AX=TK'Tl(PPMTr>j^cpv^TCHAX,NRF^') 

■'‘^MAX=I^'Tl(Rpv'rS',^o*'C,T''G''A^,NSPMO) 

CHEMA'< = CC*TFE^'fiX*PP^'TC^• 
sv••Ay=CC•T“)'AX*PP^'TC^■ 
o C '1 =C C • T T C I ♦ P P T C 
o H E ^ ' I S' = C C * T H C N' T N- •* o c (, T r 
P\CM=AMIf'l (FEC‘"TN,D‘'*'AY) 

C 

C IS PCWFR OECUIREC ►'CFF TPAh NC?'TKAL MCTO= RATING 

C 

IP(PCCM.GT,FhC“.rc,yr = S,CT ♦V‘'AX)C-0 TO lE 
C 

C K'O - TRACTTCh f'CTOP Hi^CLES EVFR^TWINr. 

C 

T£o=n. n 
poyc c= 0, 0 

T M G = T T r I . 

GO Tr 30 


r 
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0*1 CVF 


) 

0074 
000076 
) 000103 

000106 
000110 


000110 

000112 


000113 

000115 

000116 

000121 

000125 

000126 

000130 

000131 

000135 

000141 

000143 

0145 

■../0I47 

000151 

000157 

00CI57 

000160 

000162 

000162 

000164 

000167 

000171 


000171 
000172 
0 00177 
00017U 

000175 
000176 
000177 
000200 
00C2Q1 
000203 
^ 0207 
j213 
000215 
0002 ?) 
000230 


r YES - HEAT FKGir.'E GUT5 I) 

C 

15 i?P‘^EC=FPMTC 1 

17 TF(FCOH.GT. PHE^'IN + PyyAy) rc TO 20 
TE0=AMIH1(TTCIK,TFC*^IK) 

T*1C=TTCIM-TEC 
6C TO 30 

r 

C YES - TFA^TIC^ yCTOP AT KAY, HEA^ fk'GIME PICKS UP RE'^AIK'CER 

c 

20 TEO = TTCIN-T»'yAy 
thc=tmhax 
r 

C HO or 2 CALCULATIONS 

C 

30 TTCO=TSO? 

PP YTCO=PPHEC? 

CALL TOCON( l‘, 1,NLC<2) 

TACC=INTl{FFyTCK,PFVAy,TAy,NAy) 

TTCIN=TTCIM+TACC 

RpyTC?=RPMTCN 

TTCN2=TTCIK 

THEM AX=THT1 (RF»*TCN,PPH,TGyAX,NFFy) 

THMAy=rNTl{RFyTCN,SF»'C,T*'CMAY,NSFMn) 

PH£H AX=rC*TPFyAX»opyTCN 
cvYaY=cc*TyyAx*RPMTrN 
FCO“ = CC•TTCI^'♦FF‘'T''^ 

PHEyiH=CC’'T5CNN2 ■'’’PyTr*; 

IF (PCCM.GT.FHENT'nnn Tn 60 
TEC2=0 , 0 
5pv»rc2 = 0 » n 
T*'n?=7Tci»' 

GO TO 100 
40 FP Ycn.? = cPHTC2 

TEr]2 = A‘^IMl{THF*^ay,TTCIK) 

TMC2=TTCIK-TFC2 
GO TC 100 

r 

C SYSTEM FCWEF CUT°UT IS hCGATIVE 

C 

50 , TEC= n, n 
7E 02 = 0 « 0 
TMC=n, 0 
Tr*O?=0. 0 
RP»EC=n. G 
PPM50?=C, 0 
PPyTCO = RP*'SC 
cpMTf==PNTCC 
DC F 0 J = 1,1C 

TMVIK'rlMl (FFNTCtSPMC.Tf/ryiNtNSF'^O) 
7ACC=INTl(PFKTC,RDMAy,Tfly,NCy) 

TTCTy=7*'v*jK-Trrr 
fall Tc^^^ (-1,1 ,NLCK) 

TF ( £ np { ppw-r-F 3 VTC‘') , L- . = , 0 ) FO TC 65 
cf«TC=°F 
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000231 
*• V33 
.237 
000237 
000240 
000242 
000246 
000247 
0002E0 
OOC252 
000253 
000254 
000255 
000260 
000262 
000263 
000264 

000265 

000266 

000267 

000270 

000271 


0H0271 
272 
0ua273 
000274 
000275 
000276 
000277 
000300 
000301 
000303 
000307 
0 00310 
000311 
000312 
Q00314 
000315 
00C316 
000317 
000322 
000325 
0QC327 
C0P33? 
000 334* 
000337 
000344 
r 251 
L .351 
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60 CCNTImuE 

65 IF (TTCn.GT. T50ro TC 70 
TTCC=7SC 

TS02=T5C ^ 

CfiLL TGC06!(ltl,N'LCK) 

TaCC = IMTl(R = f^TCN,PF“ax,T.ex,NCX) 

TTCIM=TTCTN*TaCC 

TGC=TTCIK 

RPMTn = FPKTCN 

RPMTC2=PP>^TC)‘ 

ttcmi=ttcin 

TTCN2STTCTN 
GO TO 100 
70 TSC=TTCO 
TS02=TTC0 
TTCK'1 = TWMIK 
TTCN2 = TF'MI\' 

TGC-=T4MIN 
SPKTC1=RFPTCK 
RPMTC2=PCKTCN 
J9RK=1 
GO TC 100 
C 

C CAR IS AT PE?T AND ICLTNG 

C 

PC PRRK=0,0 
cGO=n,0 
TEC=n. 0 
TEO?=0. 0 
TTCO=n . c 
PPVFC=0, 0 

pp >jrr2 = i ) , 0 
cp‘i7f'rj-r. pv*<;(^ 

CALL TOCC)‘ (1 ,0 , KLCK) 

TACC = T^T1( = PHTCK,cp>'AX,TAX,K'AY) 

TTCT N = '^TrT^f + TACC 

TMC=TTriM 

T»-C.’=TTCIM 

PPMTCl =CPMTCk 

ppMTC2 = =P^'■^C^ 

TTCM1=TTCTV 

TTC''’ = '^TCIK 
JBR‘-<' = 0 

IOC PEC = CC ’TEC* PPVE'" 

P£02=CC*7Fn2*pcf/pf32 

PMC=CC*TMr ♦PPJ'TC 1 
PMC? = Cr'»TPC2-' = p‘'TrP 
'=GC = rr. *7GC*RP“TC1 

PLTC = CC’'AOS(TTC^1''PP^'TC1-TSC•PPVEO) 
PLTC2=CC'‘APE(TTrt.'7*5cvTC?-TSC2^RPMS0 2) 
RPTIJRM 
cN'G 
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GOLETft 


"n 0006 


00 0 006 
000006 


000006 

000006 

000007 

COOOlO 

000012 

000012 


000014 
000015 
00 0 0 16 
oocn?o 
000020 
000021 
00002? 
00 0021 * 
00 0024 
000035 
000040 
000041 
00 0047 
000051 
000054 
000056 
000057 
00 0061 
00 0 062 
00 0063 
000065 
0" 0 071 
073 
000077 

0 00 1 on 

000104 
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SUPOCUTINP TCCC^' ( JE>‘'',J57CP,KLC'<) 

CC'-fMCH nicrL,FPV(2n).Tn:r'«X{20) ,0VPD(20),P5rC(2n,2n) ,^'RP^',^!M5P♦ 
lEMUH.E“UG»FIhKLrT^P‘'Ct?0 ),TVCHflX(20) »T^O‘'IM?0)»NEP»"0* 
ZCTCZRO,TSF(2n),TTP(?0),SP'?(20),T5D2(?0)»NTSP, 

ZriLCCXC?) ,TCATIC(5>tCTl(5),CT2(5),EMUT(5) , 
4WB,Ee‘IAXtE“L'PG,E'”.JRG2, 

5CRATIO.CD1 »CCZ,FHUCf F TIP F,C TIRE 1 ,CTIRP?, VMfiS 5,DL I » CD A , 

ETECHIN tOBPAX^RPhlCL, VyASS2 tPFC*'!^’. 

7PUELSG »SKALE,7H5<L fKTHtKCOSC!- , CHGF FF , NC YCLE , NC r wo , p*^ JA V , 
eNUNITStTH5PT(5) .UPSHFTCPl.CNlHFTfF) ,0DTEFP(2C),CYCLES(?n J,KTC( 
oN'FRTC(3)»GTC(3),TFC(3)»CSUPf^0)t CMC (30),C5AR{3n) 

l,PP^'PST,PPEST,cppp(2m.oFf^c(;F(20),P.Pr'CPM,EHUC\/T.FATlJP,RATPN.K'C 
?t POEMS (20) ,EPENS (2 0) ,NGFNS 

3»VHAX,TF0MK2,l;P 5FL»CKSFL ,CWAfiMtN’AX,PFVAX (20) .TAX (20) 

COMM CM TIMC(3.2C0) ,SPEPC(3,2C0),GAMHA(3:0,3) 

COMMON fa,fp;fac. pNET.TCC, TTC, RFMDO, RPMTC.TSO.RF ms C.RoMTCO »ttc 
IRPMTCN ,TMC .TFC.TV02.TEC2 , T CO . F 5C , 0 . FEO 2 , FM C , P*' C2 , ^GC , F FR K ,FC 

ZPO.PT.PRW.FLTC.TWPPC.TTCIK! 

3, RFME0,PP^'FC2,pLTC2. VMCE 

4. TS02»PFMSC2.FSC2.»FMTC1»RFMTC2.°RMEF.RRHEP2 
PEAL IMT1.TM2 

IF ( JSTOP.NF. 0)GO TO 5 
RPMTCN=PPMICL 
tTCIK=4PPICL'‘*?/CTC7PG 
GO TO 100 

F TP (JFND.EO.-I )GC TC 50 

r 

c MAVP nUTFL!''' SFF5C AMT OUTPUT TCP<"'IF, 

c 

IF O'LCK.EO. 0)GC TO 7 
TTCTM=TTCr 
RP*^TrS' = pP^TCC 
GO TP 100 

7 TP (TTCO. C. 0) GO TC 10 
TT0IN=n. 0 
POWT(->4 = opHTCC 
GO TO 100 

IP TS°aP=PPMTCC/5IGK!(S0PT(AFS(TTCC))»TT00) 

IF(TcpAp,/-r,(.,p)r,c TO PP 

ttcin=ttco 

'OPaT=l,^/T^Tl (CaSfTcpaP) .TCP.SPP.MTSP) 

GO TO 40 

20 IP (TFPAr.rT. r.0)GO TO 
IF (TTCC.GT. 0. 0) GO TO 2’ 
pp‘^TrN: = ppMi CL 
TTCIN=PPMTCL'»*2/CTC7RC 
TTro=Tnp ( 1 ) ^TTOTK 
GO TC 100 

2? TTCIM=TTC0 /TGF (1 ) 

PF^*TQ*4=cQCT(TTCTM*CTC7Fr> 

GC TO 100 

25 TORAT=TMTl (T«^FAC,TSP,T0P.MTSF) 

TTOI*''=TTCr/TOFiT 

70 cpp/\T=T4T1 f T SF A P , Tso , pP- , N'^SP) 

t.p cpwTCNrRfNTCr/SPPAT 


GCLETO 

FCKTPf r 
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Cl. 2 HYBRID SYSTEM STARTUP SIMULATION PROGRAMS (VSYS & VSYS2) 

1* Program Description 

The programs VSYS and VSYS2 simulate the response of a hybrid 
system to a power demand. Initially all of the power requirements 
of the vehicle are being met by the electric motor under steady 
state conditions. The vehicle is travelling at a constant speed, 
the heat engine is disengaged from the drive train and is turned off. 

In response to a demand for increased power output, a clutch 
between the engine and the motor is closed and the engine is 
started. The resultant torques, rpms, vehicle speed, and vehicle 
acceleration, are determined over time scales on the order of a 
second. The simulation terminates when the engine rpm reaches 
motor rpm. 

The heat engine is modelled as a engine output torque curve 
that is a function of crankshaft displacement. The engine clutch 
is modelled as a torque curve as a function of time. Over the 
range of rpm's considered in the startup simulation, the motor 
is modelled as a straight line of motor torque versus rpm. The 
torque converter is modelled as two monotonic functions serving to 
completely determine input torque, output torque, input rpm, 
and output rpm, from any combination of two of these variables. 

The VSYS2 program version has no torque converter. The transmission 
has gearing, spin losses, and efficiencies. The differential has 
a diiferential ratio, spin loss, and an efficiency. The vehicle 
itself has aerodynamic, rolling resistance, and grade forces. 
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PROGRAM VSYS 


O 

V' 


'■\T 


^ pAGR 


THIS PROGRAM MOCELS A HYPRIO VEHTCLE PROPt'LSION SYSTE'^. 

THE ELECTRIC MOTOR IS FRCHL'CING TORQUE UNDER STEADY STATE 
CCNOITTCNS AND THE HEAT ENGINE CLUTCH IS ENGAGED. 


ALL UNITS IN HKS EXCEPT t‘ 
■^HETA, THFTAN 
V, VPL 
A, A PL 

RPMOO, PPKTC. OPYTCC, p°MTCN 
RPHHE, RPHEM 
WPL, UPL 
CM 

CT2, C02, SHFTUR, CTIRE2, TSp. 


( REVOLUTIONS) 

(KH, PER HR.) 

(G) 

CPEV. PEP HIN.) 
(REV. PER HIM.) 
(REV. PER MIN.) 
(MEWTON»HETEPS/RPM) 
TSPAR, CDA 


COHMCN RTTRE,C01,C02.DRATIC.CTl(6),CT2(e),TPATIC(6),C0A,CTIREl 
1CTIRE2,VMASS,EHUD,EHUT(6) .NGEAP.CLI , GRAD.CH.HEI.EMT, 
?VMASS2,FG,V, JGEAR, 

3TSP( 20 ) , TOR (20 ) ,SPR( 20) » NTSP. 

4 NT CL ,TCL (10 ) .TIME ( 10) ,NTHE ,THE (10), THETA (10). 

EVMPS.RPHDC, TLFD.RPMTr.TLFT .PPMTCQ.FA ,FR, FNET O.TCO.TTO.TTCO.TSF 
FTQRAT.TT Cl N, SPRAT. RFMTCN .TEMO, THETAN . 

7TCLN,THEN,TEM,FMET,paC 
REAL Y00T(4) ,Y(4) .INTI 


input E»'GTNE CATA 


READ 6**0,NTHE,HET 

FEAO 61P,(THE(J),THPTA(J) ,J=1.NT«E) 
PRINT 420, HRI 

PRINT «2 0. ( TPF ( J ), THETA (J) ,J--1 ,NT^E) 
INPUT GLUTTH CATA 


READ 62n,)'TCL 

READ 610.(TCL(J),ti*^c(j) .J = i,nTCL) 
PRINT 400 

PRINT 820, (TCL (J) , TT»-E(J) , J=1,Ntpl) 
INPUT HCTCP CATA 


PEAC FIO.CM.EMT 
PRINT U30,CM.E“T 

INPUT TCPPUF CCKVFRTCR DATA 


PEAO 620.NTSP 

READ •=)10,CTC2RO,(TSP(J).TQ = (J).S°R(J),J=1,NTSP) 

PRINT 860 

PRINT 865,CTCZP'' 

PRINT 870 

PRINT 814,(TEF(J) ,TQP{J) ,SpR(J),J=1,NTRP) 

INPUT TPA^■S^'IP^IO^' DATA 
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T223 
. J233 
000254 
000260 
000264 
000305 


000313 

000335 

000341 

000345 


000367 

000405 

000417 

000427 


000432 


000441 


000452 

00C453 

000455 

00C457 

nooup? 

000464 

onoi^E!* 

000470 


000472 

000503 

00050*+ 

000505 


■>505 
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c 

c 
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c 

c 
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c 

c 

c 

c 

r, 
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c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 


READ 640,MGEAF,SHFTUP 

READ 610,(TfiATIC(J),CTl(J),CT2 (J) ,EHUT( J ) » J=1,NGFA R) 
PRINT 040 
PRINT 845 

PRINT 810, (TFATIC( J) tCTl (J) ,CT2(J) ,EMUT( J) ,J = 1, NGEAR) 
PRINT' .)5 0,SHPTUF 

INPUT DIFFEFENTIAL AND TIRE DATA 

READ 6lO,nRATIO,C01,Cn2,EHUD,RTIRE,CTIRE 1,CTIRE2 
PRINT 875 
PRINT 880 

PRINT 885,0RATIC,Cni,C02,EHUD, RTIRE , CTIR E 1 , CT I R E 2 

INPUT' VEHICLE DATA 

READ 610,VMASS,CLI,CDA,GRAC,V 
PRINT SROtVyASSTOLItCDA 
POINT R25,GRACtV 
VMASS2 = VHASS + CLI/RTIPE’^*P 

GRAOIENT- 

FG =R. 807*^ VM ASS ‘'SIN (ATAN( GRAD)) 

INPUT RUN DATA 

READ 630, CELT, TF,NPPNT 

COMPUTE STEADY STATE CrNPTTICNS 

JGEARrl 

VMPS=V/3.e 

RPMno=9, 5 4g2067*'VMPE/RTIPF 
TLFD=C01+rC2*pPMOC • 

RPMVC=RPMCC*'nRATTO 
1!" CONTINUE 

TLFT=CT1 (JGEAR)4CT?( JGEAP) *PPVTC 
PPMTCO = PP''TC‘'TRATI C( JGFAP) 

UPSHIFT IF KECESSAPY 

TF ( ( POMTCC. LE.SHPTUP ) .OP . ( JGEAP.GE.NGEAP ) )GO TO 13 
JGEAR= JGEAR+1 
GO TO 15 
18 CONTINUE 

AERODYNAMIC - 

FA = 0.61?5'»CCA*VHOS**2 

ROLLING R5SISTANCP - 


000510 


FR=P,8n7*VyASS*(CTIPEl+CTIRE?*V) 
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000566 
000571 
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00057U 
000604 
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u - J 64 4 
000647 
00065? 
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c 
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flCCELERATTCh - 

FNETOsPA^FG+FK 

yOO=FNETO*RTtrf 

DIFFERENTIAL - 

TTOsJDO/ (FKUC* jRATIC) ♦TLFC 

TRANSMISSION - 

TTCO=TTO/(EMUT( JGEAR) •TRATIC (JGEAR) ) flLFT 

HAVE OUTPUT SPEED ANC OUTPUT TCROUE. 

TSPAR=RPHTCC/SQPT ( A9S (TTCC) ) 

TQRAT=INT1 ( T S^ AR , TSP t TO P , N TSP ) 

TTCIN=TTCC/TCFAT 

SPRAT=INT1 { TSPAR,TSP,SPR »NTSP) 

PPHTCN=PPMTCC/SFPAT 

TEMO=TTCIN+CM*PPMTCN 

SET INITIAL CCNCITTCNS 


K=0 
T = 0, 0 
Y { 1 ) =o.n 

Y ( 2) =0. 0 

Y ( 3) =. in4710755*PPVTCN 
Y (4) =V/3. 6 
PRINT 700 

en continue 

CALL R^TTAl (Y,YCOT tT ,4) 

IF ((T.GE. TF) ,CP. (K.FC. 0) )GO TO 107 
K = K-1 
GO TO 109 
107 continue 
K--NPRNT 

TIME TO PRIM RESULTS - 

ROMHE=Y ( 1) / . 10£*71975F 
RPMEM=Y( 3) /• 104719755 
V=Y(4) *3.6 
A=YDCT (4 ) /9. 207 

ORINT 886. T.RPMHE.RPMEM , FPMTCO .THCN , TCLN ,TE>" , V. A 
109 continue 

INTEGRATE TC NEYT TI^F STFO 

IF (Y (1) . GF. Y (3) ) GO TO 190 
CALL RYTT A? (V . VPQT , OFLT ,T, £*) 

TF (T.LT. TF> GC TP 80 
190 CCNTTNUF 
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000657 

000657 

000657 

000657 
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000657 

000657 

000657 

000657 

000657 

000657 

000657 

000657 

0657 

.0657 

000657 

000857 

000657 

000657 


000657 


700 STOP 

r 

C FORMAT STATEMENTS 

C 

AOn FORMAT (IHO, ’CLUTCH OAT A ’ //5X . ’CLUTCH TO, ’ ,cX , ’T I HE ’ ) 

A?0 FORMAT (IHO, ’ENGINE C AT A’ //5X ,’ENGI NE INERTIA =’,El6,6// 

15X, ’ENGINE T(3, ’,6X, ’THETA*) ' " - • • 

A30 FORMAT ( IHO, ’MOTOR D A TA ’ / /5 X , ’MOTOR SPEED CROCR =’»E15.6»AY, 
I’MCTCR INERTIA =’,E16,6) 

610 format (7F10.A) ■ , . - _ . 

620 FORMAT(7I10) 

630 FORMAT (2F10, A, 21in) 

640 FORMAT (110, 6F10. 4) ' 

700 format ( IHO ,4X,’TIME’ , 7X, ’ENGINE RPH’ ,3X , ’MOTOR RPM’,4X,*TC0 R 
16X, ’ENGINE TC,’,3X, ’CLUTCH TC,’,3X, 

2’HOTCR TQ,’,4Xt’SPEEO’,8X, ’ACCEL,’) 

810 FORMATdH ,4E16,6) 

814 FORMATdH ,2E16.5) 

820 FORMATdH ,2E16*6) 

etO FORMAT ( IHO, 12HGEARROX OATA) 

845 FORMAT(lH0,4X,10HGEAR P A TI C , 6X , 22HSP I N LOSS CCEFFI CIE NT S , 1 OX 
1 , ’cFFICTEhXY’) - - - - 

85 0 FORM at (/T5 , ’UPSHIFT PPM = *,E16,6) 

860 PORMAT (1HO,21HTCROUE CONVERTOR CATA) 

86= FORMAT (IHO, 4X,17HNI”2/TI(STALL) =,E16,6) 

870 FORMAT (IHO, 4X,11HN0/S0RT (TO) ,5 X ,EH T 0 /TI , 1 IX , 5 WN C /N I ) 

875 format (IHO ,18HAXLE AND TIRE DATA ) 

een =ormat dHa,4X,it*HFiN, or. ratio, 2 x , pzhsptn loss ccef=tctemts, i 
1, ’ EFFICIENCY*, 6X, ’ROLLING PACT US* , ?X , ’ROLL. RESIST, CCEFFICTEN 

885 FORMATdH ,7E16.6) 

886 FORMATdH , E 1 2 . 2 , 9E1 3, 4 ) 

8.°0 FCRMAT(lHn.l4HVEHICLF MASS =,E16,6 
1,SX, ’DRIVE LINE INERTIA = ’,E16.f 
2,5X, 19H(0PAG COEF.)’AREA =,F16.6) 

P25 FORMAT (IHO , 16HCYCLE GFACIENT =,E10.3, 

I’STEADY state SPEED =’,E10,3) 

END 


000657 
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000007 


000007 

000007 

000015 

000017 

000023 


000031 
000034 
0036 
•^ 00040 
00004 ? 
000044 
000047 
000053 
000057 
000064 
000066 
000074 

oooni 

000104 

oooloe 

000110 

000112 

ooniiA 

000115 


SUePOUTINE CFPTV (Y ,YDCT, T,N) 

c 

C Y(1) IS ENGINE SOPEO 

C Y(2) IS ENGINE ANGULAR CISPLACEKENT 

C Y(3) IS MCTCR SPFEo 

C Y(4) IS VEHICLE SPEED 

C ' 

COMMON RTlRE,C01,C02,nRATI0»CTl(6) ,CT2(6)tTRATIC(6),CCA, CTTPE 
1CTIRE2,VMASS,EMUC,EMUT (6 ) , NGE A R , CLI , GRAD ,CN,HEI ,EMT, 
2'/MASS2,FG, V,JGEAR, 

3TSP(20) ,TOR(20) ,SPR(20) ,NTSP» 

**NTCL ,TCL (10) ,TIVE(10) ,NT HE , THE ( 1 0) , T HETA ( 10 ) , 

5 VMPStRPHCC,TLFD,PPHTO,TLFT ,RPMTCO,FA , FR , FN'ET 0 , TCC , T TO , TTC 0 »TS' 
eTORAT,TTCIN»SFRAT»PPHTCN,TEHO,THETAN, 

7TCLN»THEN,TEyTFNET,FAC 
PEAL Y(2) ,YCCT(2) ,TNT1 
C 

TCLN=INT1 (T,TIME tTCLTNTCL) 

THET ANrY (?) /e. 283185? 

THEN=INT1(THETAN,THFTA, THEtNTHE) 

YDOT (1) = (THEN + TCLN) /HEI 
C 

C NO negative ENGINE OCTATICNAL ACCELERATION 

C 

IF(YDOT( 1) ,LT.O. 0) YOOTd ) = 0, 0 
YOOT(?)=Y(l) 

PPMnO=g.54p?oE7¥Y(4) /R^IPE 

PPMTC = RP>^nC*CRATTC 
RPMTrO=PPHTC'TRATIG(JGFAO) 

SPRAT=, 1 0A719755*RcmTCC/ Y( 3) 

TSPAR=INT1(SFRAT,SPR,TSF,NTSP) 

T0RAT=INT1 { SFRAT,SPR,TQR ^NTS®) 

TTCO= (PPMTCC/TSP AR) »*?, 

ttcin=ttcc/torat 

TEM=T£M0-C“*Y(3)/.10A719755 

YDOT(3)=(TEy-TCLN-TTCIN) /EHI 

TTC = TTC0»EHUT ( JGEAR) *'TRATIC(JGFAR) 

TOC=TTn»E>"LC*CRATIC 

fnet=tcc/ptirf 

fac=fnet-fneto 

YOCT (4) =FAC/VHA5S? 

return 

F.NC 
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UU0007 

000010 

000011 


SUS‘?CUTIS’F PKTTA 1 (Y1 ,£*<1 ,T ,K') 
REAL YK20) ,EK1 (20) 

CALL OEPIV(Yl,E'a,T,L’) 

RFTUPN 

ENm. 
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SUBROUTINE RKTTA?(Y1,EK1,CELT,T,N) 

real Yl(20),Y2(20) , YTE'^P (20)»EKK20) ,E<2(20) ,EK3(20),EK4(20) 
00 50J=1*N 

YTEMP( J)=Y1 (J) +CELT*5t<l ( J) /2. 

50 CONTINUE 

TT=T+0ELT/2. 

CALL’ OERIV( YTEHF,EK2,TT,Nr 
CO 60J=1,N 

YTEMP( J)=:Y1 (J) fCELT_*E»<2( J) /2. 

60 CONTINUE 

CALL OERIV( YTEMP,E!<3 ,TT,N) 

00 70J=1,N 

YTENP( J) =Y1 ( J) + CELT*EK3 ( J) 

70 CONTINUE 
T=T+DELT 

CALL 0ERIV(YTEHP;EK4,T,N) ■ ■ 

00 80J=1,K 

Y1(J) = Y1(J) + (EK1(J)+2.*(FK2(J»*E»<:3(J))+EK4(J)) (OELT/6, ) 

80 continue 

RETURN 

END 
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RF.AL FUK'CTTCN INT I ( EX , X , Y , VX ) 

C 

C INTI INTEPOCLATES IN CNE VARlAeLE 

C 

REAL X(10),Y(10) 

IF (EX.GE.X(l) )GC TO 5 _ _ 

JRT = 2 ■ 

GO TO PO 
5 DO 10 J=1,NX 

IF (EX.LT.X ( J) )GC TO 15 
10 CONTINUE 
15 JRT=J 

20 IMTl=Y(JRT-l) + (FX-X C JRT- l) ) ♦ ( Y ( JRT ) - Y ( JR T- J ) ) / ( X ( J RT ) - X { JRT-1) 
OETUFN 
END 
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GO LET A 


000003 


000003 


000003 
00001 3 
000030 
000036 


"00053 
0 00 061 
000076 
000102 


000117 

000127 


000137 

000147 

000170 

^'"174 

200 

000221 


PCKTPAh 1,3 • SEHl-iUTO 'JPL • (01-10-73) 


PROGRAM VSYS2 
C 

C THIS PPOGPAK HCCELS A HYPRIO VEHICLE PROPULSION SYSTPM, 

C THE ELECTRIC KCTGR IS PRCDUCING TOROUE UNDER STEADY STATE 

C CONDITIONS ANC THP HEAT ENGINE CLUTCH IS ENGAGED. 

C VHYSYS2 VERSION HAS NO TORQUE CONVERTER. 

C 

C ALL UNITS IK N»<S EXCEPT » 

C THETA. THFTAN 

C V, VPL 

C A, A PL 

C RPNOO, ROMTC 

C RPMHE. RPMEH 

C WPL. UPL 

C CM 

c CT 2 . cne, SHFTUP, ctipes. tsp. tspar, 

c 

COMMON RTIRF.COI.CDE.DRATIC.CTKE) .CT2(6)»TRATIC(6) .COA.CTIREl , 
1CTIRE2,VMASS,EMU0,E*'UT(6 ) .NGEAR.OLI , GRAD .CM, H El .EMI. 

2VMASS2 ,FG,V, JGEAR, 

3T30.EMI2. 

4NTCL.TCL (10) .TIHE(IO) .NTHE.THE (10) .THETA (10) . 

5VMPS .R PM CC .TLFD.RPMTC, TL FT .PPMTCO.FA .FR, FNET 0 . T CO . TTO , 
f TEMO .THFTAN, 

7tcln,tmen,tev,pket,fac 

REAL YOOT(3) ,Y(7) .INTI 
C 

C INPUT ENGINE DATA 

C 

PEAD 640.NTHE.HFI 

P£AO 610,(Tv-Ff I), T HE TA(J),J=1, NT HP) 

PRINT 420, HEI 

PRINT 820 , (THC (J) , theta ( J) ,J = l ,N'thE) 

C 

C INPUT CLUTCH TATA 

r 

REAP 620.NTCL 

REAP 6 to, (TCL (J) . TIME ( J) ,J = 1 .NTCL) 

PRINT 400 

PRINT 820. (TCL (J) .time (J ) . J= l, NTCL) 

C 

C INPUT ‘'CTCR DATA 

C 

PEAO 610, CM, EMI 
PRINT 430,C*'.E“T 
C 

C INPUT TRANSMISSION DATA 

C 

READ 640,MGEAR,SHPTUP 

PEAD 610.(TRATIC(J),CT1(J),CT2(J),EMUT(J),J=1,NGFAR) 

PRINT S40 
PRINT 8t«5 

PRINT 810 , ( TRATIO (J ) ,CTl ( J ) .PT :> ( j) , EMUT( J) , J=l, NGE AR> 

PRINT *'E0.SHFT!.l = 

C 


( REVOLUTIONS) 

(KM, PER HR.) 

(G) 

(PFV, PER MIN.) 
(REV, PER MIN.) 
(REV. PER MIN.) 

( NFWTnN*METERS/PPM) 
CDA 


GOLETA FCFTCAK 1,3 • SE*^T-AUTO 9FL 


(Cl-10-73) 


VSYS7 




) 

0227 
U00251 
) 000255 
000261 


000303 

000321 

000333 

0003(»3 


0003A6 


000355 


000366 
000367 
000371 
1373 
. .0376 
0 00‘*00 
OOOAOO 
000^*0u 


0 0 0 4 0 c 
onot.17 
0 OCt ?o 
000A21 


000421 


000424 


000430 

000433 


C INPITT OIFFPceKTI AL AMC T IPE HATA 

r 

PE AO 610,CRATIC,C01 ,CC?,EHUC,RTIFE ,CTIPE1 ,CTTPE2 
°RINT 375 
PRINT 880 

PRINT 885,nPATI0,C01,CP2 ,E NUO, FT IRE , CT IP E 1 ,C T IP E 2 

C 

C INPUT VEHICLE CATA 

C 

READ 610,V6ASS,CLI»CDA,GRAC,V 
PRINT 890» VHASE,OLI,CCA • 

PRINT 925,GRAC,V 
VMASS2=VMASS+CLI/RTIRE**2 
C 

C GRAOIENT- 
C 

pG=9,«07*VMAES*STN(ATAH(GRAO)) 

C 

C INPUT PUN DATA 

C 

READ 630, CELT, TF,NPRNT 
C 

C COMPUTE STFACY STATE CCNCITICNS 

C 

JGEAP=1 

VMPS=V/3.f 

PPMnc=9.r492967*VMPS/PTIRE 

TlFC=CD1<-CC2*PP“CO 

P°MTC=RPMCn''CRA'^I'' 

15 CONTTNt.'E 

TLFT = CT1(JGEAR)*CT2(JGEAP)'»PPMTC 
Pdme;m=;Ppmtc»TRATTQ( JGEA<=> 

c 

r UPSHI^’T IF KECFSSADY 

C 

IF((PPMEM . LE.SPFTUo’I ,CF . { JGEAP.GE. vGEAP n GC TO 18 
JGEAP=JGEAotl 
GO TO 15 
IP CONTINUE 
C 

r AEROCYMAMJC - 

r 

'^A = n,6125*CCA»VY«S**2 

C 

C ROLLING RESISTANCE - 

C 

FR=9.8 07*V'‘AS£» ( CTIRF i+C IRE 2^ V) 

C 

C ACCELEPATICK - 

C 

FNETn=FA + FG-»FR 
TOO = FNFT O^RT IFE 
C 

C OIFFEPENTT AL - 

C 
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'/SYS? 


000435 


000441 

00044E 


000451 


C 

C 

C 


C 

C 

C 


C 

C 

C 


TTOs^no/ (evucopatic) ♦Turn 
transmisstck - 

TSC =TTO/ (EMUT ( JGEAR) *Tp AT 10 ( JGE AR) ) «-TLFT 
TEMO=TSO+CM*RFHFM 

COMPUTE EQUIVALENT INERTIA DRIVE T^AIN 

EMI2=EMIf (VHASS2*RTIRE**2. 

1 /(EMUO^EMUT ( JGEA9) »(nRATT0*TRATI0(JGEAP))**2,n 

SET INITIAL CCNPTTTCNS 


000466 



>< = 0 

000467 



1 = 0. 0 

000470 



Y(1)=0,0 

000471 



Y ( 21 =0. 0 

000472 



Y(3)=. 104719755*9PMFM 

000475 



“RINT 700 

000500 


eo 

CONTINUE 

000500 



CALL RKTTAl (Y,YCOT»T,3> 

000503 



IF({T,GE.TF),CR. (^.FO. OllGC TO 107 

000513 



K=<-1 

000514 



GO TC 109 

000514 


107 

CONTINUE 

0514 

r 


<=M0CNT 


c 

p 


TIME TO PRINT results - 

000515 



RP«wF=Y( 1) /,10471975«^ 

000517 



RPMEM=Y( 3)/, 104719755 

000521 



V = V(3) »RTIRE*3. 6/ (rPATIC*TRATIC (JGEAPn 

000526 



A = YDOT(3)*RTIOE/(nRATIC*TRATIP (JGEAR) * 9 . 8 071 

000533 



PRINT aeS.TtRpMHEtRPMEM.THENiTCLNtTEM.V, A 

000556 

p 

109 

CONTINUE 


c 

p 


INTEGRATE TC KEVT TIME STEP 

000556 



IF (Y (1) , GE. Y 13) 1 GO TO 190 

000561 



CALL RKTTA2 (Y,vPOT,DELT,T,3) 

000564 



IF (T.lT, TF) GO TO 80 

000567 


ISO 

CONTINUE 

000567 

p 

200 

STOP 


c 

p 


FOR*"AT statements 

000571 

L 

400 

FORMAT (IHO, ’CLUTCH D ATA * //EX ,* CLUTCH TO,’,6Y,*t 

000571 


420 

FORMAT ( IHC , ’ENGINE DA TA * //5 X , ’ENP INF INERTIA =’ 




15X, ’ENGINE T0.’,6X»’7HETA’5 

000571 


430 

FORMAT (1 HO , ’MCTCR 0 A T A ’ / /t x , ’M CTCR SPEED DROOP 




I’MOTCR INERTIA =’,Fie.6) 

571 


f 1 0 

POR“AT (7Fn,4) 

0un571 


620 

FOPMAT (7 n 0 1 

000571 


630 

PQRMAT (2^10. 4, 2110) 


4Y , 
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000571 

0571 


000571 

000571 

000571 

000571 

000571 

000571 

000571 

000571 

000571 

000571 


000571 


000571 


FtO FOOHfiT fliO.e'^lO.t*) 

700 tlHO. AXt’TIMP* , 7X, ♦ENGINE R?H» . 3X , •yOTC‘? RPN*, 

1^X,*ENGINE 70, *,3t, ♦CLUTCH T0.*,3X, 

2 •MOTOR TO. *,4X, •3PEEC*,»!X,*ACCEL.*» 

810 FORMAT (IH ,4E16.fO 

820 PORMATdH .2E1&.6) 

fli.O' FORMATdHO, IZHGEARBOX CATA) 

845 FORMATdHO,4X .lOHGEAR R A T I C ,6X , 22HS P IN LOSS COEFF IC TENT S , 1 OX 
1,»EFFICIENCY») 

850 FORMAT (/T5 , *UOSHIFT RPM = ♦.EIE.F) 

87® format (1 HO. 18HAXLE AND TIRE CATA ) 

880 format (IHO , 4X,14HFIN, HR. R A TI 0 . ?X , 2 2HSP I N L OSS COrFFIC lENTS , 1 OX 
1 .•EFFICIENCY". 6X, •ROLLING R A Cl US^ . 7X . *ROLL . RESIST. COEFFICIENTS 

885 FORMATdH .7E16.5) 

886 FORMATdH , F12. 2 . 9E 1 3. 4 ) 

"90 format (IHD , 14HVEHICLE MASS =.E16,5 
1,5X,*0®IVF LINE INERTIA = •.Eie.6 
2 .5X. 19H( DRAG CCEF.)^ARFA =,E16.6» 

925 FORMAT (IHO. 16HCYCLE GRACIENT =,Ein.3. 

1*SYEA0Y state speed s»,E10.3) 

END 
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C 

c 

c 

c 

c 

000007 


000007 

c 

000007 

000015 

000017 

000023 

C 

C 

C 

000031 

000035 

000036 

0042 

0047 

000047 


SU3<?CUTINF CEFTV (Y, YOCT, T,Kj) 

Y(l) IS ENGINE foPEG 

Y(2> IS ENGINE ANGULAF CISPLACEY^nT 

Y(3) IS MCTCR SFERD 

COMMOM RTI‘ E,Cni,C0 2,DRATICtCTl(f ) ,CT2(6) , 1RATIC(6) ,CDA,CTIRE1, 
lCTIRE2,VMASS»E“.un,FYUT(6) ,NGEAR,PLI tGRAO ,CY, HEI ,£MI » 

2VMASS2 ,FG,V, JGECRt 
3TS0,EMI2, 

‘*NTCL,TCL(10),7IyE( 10) ,N'7l-EyTHE (10) , THETA (10) , 

5 VH PS ,PPhOC »TLFO,RPHTc,TL FT tRPYTCn, FA , FR, FNETl) ,TCO,TTO, 

6TEM0 tTHETAN’t 
7TCLN,THEN,TEN,FNET, PAC 
PEAL Y (2) ,YCOT (2) ,INTl 

TCLN=INT1 (T,TIMc fTCL.NTCL) 

THETAN=Y (2) /6, 203185? 

T HEN =I NT 1 (THET AN ,THETA,T HP ,NTHE) 

YDOT (1 ) = (THEW+TCLN) /HET 

NO NEGATIVE ENGINE OOTATICNAL ACF.E LE R ATI ON 

YOOT (l)sAHAXl(0.0»YCOT(l)’ 

Y00T(2)=Y(1) 

TEV = TP>!0-CH*Y( 3) / .I 04715755 
YCOT(3)=(TrK-TrLN-T?0)/E“I2 
RE TURN 
END 



GOLETA 


10010 
-00010 
000011 
000017 
000022 
000024 
000D27 
000034 
000041 
000045 
000050 
000055 
000061 
000064 
000065 
000070 
000075 
000107 
COO 111 
000112 
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SUBROUTINE «KTTA7(Y1»EK1 ,DELT, T,N) 

REAL Y1 (20) ,Y2 (?0) , YTEVP (?0 ) .EKl <2 0) »EK2 (?0) ♦ EK3 (20) »EK4 (? 0) 

CO 50J=1,N 

YT£MP( J) =Y1 ( J) +CELT-EK1 ( J) /2. 

50 CONTINUE 

TT=TfOELT/2. 

CALL OERIV( YTEMF,EX2»TT»N) 

00 60J=1,M 

YTEKP( J) =Y1 { J)+DELT*EK2( J) /2. 

60 CONTINUE 

CALL 0EP.IV( YTEMF,EK3,TT,V) 

00 70J=1,N 

YTEMF( J) =YK J) ♦DELT*EKY( J) 

70 CONTINUE 
T=T+OELT 

CALL DESIV( YTEHF,EK4,TtN) 

DO 80J=l,N 

Y1(J) = Y1(J) + (EKl (J)+2,’ (EK2( J) ♦F'<'3( J) )+EF4(J) )* (HELT/F,, ) 

80 CONTINUE 
RETURN 
END 
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000007 

000007 

000011 

000012 

000013 

000015 

000020 

0U0&22 

00002 ^* 

000040 

000041 


c 

c 


REAL FUrCTICM T T 1 1 F X , 7 , V , NX ) 

TNTl TNTERFCLATFS TN CNF V AR I A EL F A' ' ^ ’ 

( i’ 'i'll ^ 

REAL X (10) ♦ Y (10) 

IF (EX.GE.X(l) )GC TO 5 
JRT = 2 
GO TO 20 
5 DO 10 J=1,NX 

IF(EX.LT.X( J) )GC TO 15 
10 CONTINUE 
15 JRT=J 

20 IMTl =Y ( JRT-1 ) ♦ (FX-X ( JOT-1 ) ) * ( Y ( JO T ) - Y ( JR T- 1 ) ) / ( X (JRT ) -X ( jR T- 1 ) ) 
RETURN 
ENC 
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CjF poor QllAUTY 

Cl. 3 CRASH SIMULATION PROGRAM (CRASH) 

1. Program Description 

CRASH is a program which simulates the impact response of a 
general non-linear spring/mass system. The masses are constrained 
to translation along one axis; i.e., planar or 3 dimensional 
translation and notations are not simulated. The springs which 
interconnect the masses can be nonlinear and have different 
characteristics in unloading than they do in loading; i.e., they 
can dissipate energy in heat as well as storing it. The equation 
of motion of the masses are solved using a Runge-Kutta integration 
routine. By appropriately modelling a vehicle (or vehicles) as such 
a spring/mass system, the vehicle response in a barrier or vehicle- 
to-vehicle impact can be simulated. 


A 3.2 Equations 

th 

Equation of Motion for i Mass 


M. 

1 



dFI 


L 1 


ij 


- E 

i/: 




i=l, ■ ■ ‘n, 
j= 0 ,"‘n 


where n is the number of masses in the system and is the force 

generated by a spring coupling mass i to mass j . The sign convention 

on the F. . is such that F. . >0 if it is in a direction to accelerate 
ij 13 

mass i and decelerate mass j . See Figure A 3-1 for an example of 
the nomenclature and sign conventions. 


Generation of Spring Forces 


F, = 


f J . (x.-x.), loading with x.-x. - x 
ij J 1 ® 31 sav 

dF 

MAX ( 0,f (x ) - (x - (x.-x.)) 
’ sav dx sav j 1 


unloading or loading with x^-x^< 


X.^ Changes 
ij 

Sign (- to +) 


Figure A3--1 Generation of Generalized Non-Linear 
Spring Forces 




C-llA 


where x is the last value of x. *• x. at which unloading started 
sav J i 

(i.e., Xj - x^ changed sign from positive to negative) and dx is the 
unloading slope. See Figure A for an illustration of how is 
generated. Note that, if there is a spring element between and 
such that the force can go negative for negative values 
of deflection, it is necessary to represent it as two springs. 

For example, the coupling between and M 2 in Figure A would be 
represented by two functions f ^^2 ^21 
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'PPEPnoETEES 


EounnoN 

PfiOGFPm 

Curl T2 


DESC^JPTJ Oh 




A/C/I5S. 

Nrr)f!5S 

t^dCLLP 

Df^F 

£/r)/95S^/d) 

NCMLCPTSiso^ 


t (6 q, to ) 

\F(5o,td) 

SLC-pt(5:>) 

TP 

D'eMT 
\ h rRH T 


m 

hjm 


6EC 

SEC 


nnjiec^ 


Nurr\6tfl Of C/)5S'j 
fimm G)EA of /tdFI^sSl 
N iLmBE^ or CoaPU NG5 
Dh>NPrnid F) /rtVci F/ICTOF. 

/>?/)53 OF SLZn)BHl5 

CoLULmHO ± Fno cS Cj>N7FiN /r)fi5 5 

NumtSE^S OF OOUPLEO Fr)fit 5 Lt^ TMIFD 

CoLULmN Coh'TFJNS /^um3£A or 

TN f fitJD X F^/ZlCf)tf 5 FO^ oPSClFje.D 

COUlPLLNC 

DEPLEcTion vfllles ro(^ epc.fi (^^uplihc 

PoFCB VflLU£3 FOfC fO/2C2. VSPLdC.l LON 
CFfiCF(LT£JlJ5TrC5 FOP EPCH COUJ^LINC 
Scope of cine PLoriG UPf^lCU UN - 
LOFOlHG occurs FOE SPcF COUPLING 
FIN Pc TxrnE 
Tims JNCEEmENT 
fESQUShiCh^ OP OalVui ~?C1NT1NQ 
CxN ^SLpiiDN 70 7? mt. Step) 

INin Pl VELOCIIN^ 



i 
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T)£5cei 'PT rah' 


JN-b-nRTfiliSous vntue. oF ]/Lloci'^'^ roCefiCH CouPlii 

THST ftTf)NLOU3 value OF Dr FO/^ niATlOt'^ 

Fo/i. Coupling 

ySLoCXTV Of /T)A55 (j,P^ . ,, lo) 

ACCELe CATION OP mfi55 (l,2, . . . lo) 
vALut OF Foccs AT Point loNece 
(UNLOADING BZGTNS 
iNsrANTANcOa^ VALUe. OF FOUCc. 

EXFXVEJD FO^ FJ\CA CdUPLZNG 
FLFG5 OCCllF/^ENCL ^OF UNLOADING 
FoC. EAcA Coupling (unloading -- i ^ 
Loadjlno - ±) 

STORAGE Fan (^AaiAeLEt OF INTtGFFlii 
.NtLrm BeC. OF /)')A55E5 
NumBSfZ. OF TlmZ SIEVE 
VALUt OF BE FaFTAAlI ON AT iOVlCV 

unloading BeGiNa 

DISPLACE rnF NT OF mms (l,0, ... lo) 


\ 


1 

! 
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C/P/9 5 H 

ST A AT ^ 



ti£fl AC 

V/)/ti»6L£.S 


.yCauP 


''i' 1 • ■' "I (.'rALlTV 


He.ADJNS 


© 


NPTS * 

Ndou.PT 8 IS, 3) 


NC.8SS. 


(^(.S,k), FlS,<), 
Ks-t, NPTS 


JdAS&^l 


NC8SZ 


src P 


/V/*78SSNCOaP ! 


SLOPUCs) 


NPTS 


F(S, K ) = 
F(T,k)* /ooO.t CAP 


/ 

DAt 

1 

y~ 

f. 


£./nPssCs), 

j* I, N/nAi/s 


'SLoP^iS) - 
iu>Pdis^^ tooo. 


(I, £/n/JSS (z\, 

1^1, N/nPSs) 


<(NCouPTB(J.K\ 

K*l,Z),S^i,NdouP) 


NLolLPTQ (S,i)^ 

NCouFTB(S,a), 
ScoP£.(7) 


PCS,e.) 
A » / HPTS) 
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GOLETA 


000003 

000003 

000003 

000003 

000003 


000003 

000007 

000017. 


000022 

00002'r 

00033 
J00041 
000054 
000060 
00007<^ 
000100 
000120 
0 0 . 012 ? 
000124 
000145 
000153 
000155 
00016? 
0Q016‘* 
000167 
000205 
000225 
00023t* 
000242 
000250 
000 260 
11027 3 
000277 
000313 
000317 


fCPTPAK 1,3 • TEHI-AUTO PFL * (01-10-73) 



C 

C 

C 


C 

C 

C 


C 

C 

C 


PROGFAM CCASH 


page is 

POi>[; (^CAUTY 


-CRASH SIMULATES THE 3EHAVLCR. OF UP TO 10 CCNNECT.ED MASSES IN AN 
impact SITUATION, CALCULATING THF 01 SPLA CEMENT, VELOCITY, ANC 
ACCELERATION OF EACH MASS AND THF FORCE EXERTED ON THE CCUPLEO 
„ -MASSES, 


DEFINE VARIABLES 

COMMON T^AsiTl?) 7FlTo^lb)TxT50,10) ,FT(50) ,XSAV(50) ,FSAV (50) , 
ISLOPF(SO) ,NCOUPT8(50 ,3) , IFLAG ( 50 ) , OE L X (5 0 ) , DY (2 0 ) , NMASS T , NCOUP 

-..REAL Y1(20).,KJ(20) 

REAL TPL (203) ,XPL (10,200) ,XPL1 (203) ,OXPL (203) , C0XPL(2Q3) 
DIMENSION NPMASS(IO) ,NSYM8 (10) 

. DATA .NSYMB./l,2,6,7,10,12,14tll,4,26/ 


PRINT HEADINGS 
PRINT 500' 

READ 400,NCASE,NPLOT 

-_IF.(NPLOT,EQ,l)CALL PLOTS 

CASE LOOP 

DO 300 JCASE=1,NCASE 
READ 400, NMAS’=‘, N'CCUP 
READ 4?0,OAF 

READ ‘‘?0, (EMASS (J) , J = l, NMASS) 

PRINT 502 

PRINT 503, (I, EMASS(I), I=1,NMAS51 
PRINT 505 

READ 400, ( (NCOt'PTBt J,K) , K = l,3), J = l, NCOUP) 

DO 10 J = l, N'CCUP 

MPTS=NCCUPTB{J,3) 

READ L20, (X(J,K), F(J,K), K=1,NPTS) 

. READ 42C, SLOPE (J) -.1 . . 

DO 15 K=l,NPTS 
F(J,K)=F(J,»0*1000. »DAP 
15 COMTINUE 

SLOPE( J) =SLCPE(J)'»1000 

PRINT 506 , NCOUPTB( J,l) ,NCCUPT8( J,?) ,SLOPF{J) 
10 PRINT 507, (X(J,K), F(J,K), <=1,N0TS) 

PF.AO 420, TF, CELT 
READ 400, NFPNT 
READ 420, VO 
READ 4 00 , NPM, N3P‘" 

READ 400, (NFMASS( J) ,J=1 ,NFM) 

PRINT 503 

PRINT 509, CELT,TF,VO,N'PPNT 
PRINT 510 
PRINT 515 

INITTALIZE VARIABLES 


000321 


T1 = C. 0 
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GOLETA 


000324 
00032& 
..00 0 3.3 3. 
000334 
00033E 
JJ 0 0 J_36 
000337 
000342 
_00.0 3.4‘*. 
00034E 
000350 
000351., 
000353 
000354 
000355 
000356 


000361 

000362 

000363 

000366 

000370 

000371 

000373 

000375 

000376 

000404 

000406 

000411 


000414 

000416 

000424 

000425 

000427 

000440 

000452 

000456 

000462 

000465 

000466 

000477 

000503 

000506 

000510 

000514 

000517 


FCRTPAN 1.3 • SEMT-AUTO =>FL • (01-10-73) C^ASM 


NHAS5T = N'MASS*2 
NTIME3=(TF/DELT+0.5)+l 
. . JPL=0 . . 

JPRNTsO 

00 20 J=1,NMASST 

„ .Y1(J).= 0.0 . __ _ 

20 OY(J)=0.0 

00 25 J=2,NKASST,2 

... jri(J)=_vo_ _ _ 

25 CONTINUE 

00 30 J=1,NCCUP 

.. -IPLA.GC J). = l. . „ _ 

OELX (J) =6.0 
FT(J)=0.0 

XSAV(J)=0.0 _ _ _ _ 

30 FSAV(J)=TJ.d 

C 

C LOOP TO INTEGPATE EOUATICNS, STORE PLOT RECORDS ANC 

C RECCPOS 

C 

JPRN_T=.l _ _ _ 

DO 100 JTiyES = l.NTI)'Ts 

CALL RKTTAl ( Y1 , K1 , T 1 , N^ASST) 

CALL. FLAG<Y1,NCCUF) 

IF('npL0T.EQ. 0)GO TO 106 
JPL=JPL*1 
TPL(JPL)=T1 
00 105 K = 1,ND*' 

XFL (K, JPU =Y1 ( ( NPMflSS (K) ♦NPMftSS (K) ) -1) 

continue 

CXOL(J°L)=vi {N3PR*-N3FM) 

CQXPL ( JOL) = CY (M3PM4-M »P'') 

PRINT ROUTINE 

IF(JPRNT,NE. JTIMFS)GC TO 110 
PRINT 520, T1 
L = 1 

DO 120 J=1,N4ASST,2 

PFINT 525, Y1 (J) »Y1 ( J+ 1) ,OY { J + n 
PRINT 540, NC0UPT8(L,1),NC0UPT-3(L,2) ,FT(L) 
PRINT 550 
L = L + 1 

IP(L.GT.NCOUO)GO to 140 
. DO 130 K=L,NC0UF 

PRINT 540 , NCCUFT? (K.l ) ,NCOUPTO(K ,2) ,FT (<) 
PRINT 550 
CONTIf'UE 

JPRNT = JPPN'T+NPRNT 

CALL RKTTA2 ( Y1,K1 ,CELT ,T1,NMASST) 

CONTINUE 

IF (NPLOT.EO. 0) GO TO 300 
C 

C PLOT 

C 


105 


C 

C 

C 


106 


120 


13 0 
140 
110 
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T00520 

J00523 

000526 

000531 

00054S 


CALL SCALE(TPL,10.0,JPLtl» 

CALL SCALE(CXPL,10.0,JPL,1) 

CALL SCALE(CDXPL*10.-fl.JPL*l) ... - -- 

call AXIS(0.Q»0.n,4HTTME,-A»10.0,0.0,TPL (JPL+1) tTPL (JPL + 2) »0) 
CALL AXIS(O.OtO.O,23HVELCCITY CF MASS ,23,10.0,90.0, 


000561 

000575- 
000601 
000605 
0006UI 

■ — 

000615 


000617 


000620 



000621 


000625 

210 

000627 


000631 

200 

000643 


000645 

_.300 

000650 


000652 

410 

000652- . . 

-.4 0 0. 

000652 

420 

000652 

500 

00652 .. 

502 

000652 

509 

000652 

503 

00065? 

505 

000652 

506 

000652 

507 

000652 

506 

000652 

510 

000652 

c , c 

000652 

520 

000652 

525 

000652 

540 

000652 

550 

000652 



CALL AXIS(-C.5,0.0,27HACCELE9ATICN OF MASS 
100XPL( JPL+1) ,C0XPL(JPL+2) ,0) 

XALU NUKaEP.(.-.28 ,5. 2, 0. 15, N3PM ,90. 0 ,2HI2 ). 
CALL NUMBER(-.78,5.5,0.1S,N3FM,g0. 0,2HI2) 
CALL LINE(TFL,DXPL, JPL,1,5 ,0) 

.CALL. LIN£..(TFL,CCXFL,JFL, 1,5,5) - 

XAXIS=-0,5 
DO 200 JPM=1,NPV 

0.0.2 10 _L.EL=1,JEL. 

XPLl (LPL )=XFL (JFM, LPL) 


XAXIS=XAXISr0.5 - . . 

CALL PLOTl (TPL, JPL,XPL1, NPMASS ( JPM ) , NSYM 9 ( JPM) , X AX IS) 
CALL PLOT(15. 0,0. 0,-3) 


, 27, 10. 0, 90. 0, 


STOP 


PCRMAT (IH- ,1 OX, ♦COUPLING 5 X ,• UNLOADING SL CPE ♦ , 5 X , • Df.PL E CT I ON 
IFORCE* ) 


1TY»,5X,*PPINT INTERVAL*). 


IMASS J,10X,21HPCRCE BETWEEN I AND J) 
FORMAT ( IHO) 


7 0X,I2,9X , 12 ,1AX ,E1G.6) 


EN! 
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.000C07. 

000007 

000010 

.000011 


SUeoCUTINE R<TTA1(Y1,EK1,T,K') 

c 

C_. RKTTAL. SETS EKl. EQUAL TC INITIAL CONDITIONS FOR USE IN 

C INTEGRATION RCUTIN'F ^ 

C 

peal YK5.0) ,EK1 (£0.) 

CALL OERIV(Yl,EKl,T,N> “ “ 

RETURN 

- _ . c.N0 
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GCLETA 


000010 
. 000010 . 
000011 
000017 
000022 
000024 
000027 
000034. 
000041 
000045 
00C050. 
000055 
000061 
000064 
000065 
000070 
00 007.5. 
000107 
000111 
000112 . 
000112 


PCR7FAS 1,3 • SEf'I-AUTO 9PL • (01-10-73) 


SUePnUTINF KKTTA2 (Y1 ,EK1,0ELT, T,N) 

C 

.C 5KTTA2 USES CUMGE-<UTT A. TECHNI CUE. FOR NUMERICAL INTEGRATION. 

C 

REAL Y 1(50) ,72(5 0) ,YTEMP (5 0 ) ,E K1 ( 5 0 ) ,EK2 ( 50) ,EK3 (5 0 ) , E><4 15 0 ) 

no..5 0J=l,.M ■ 

YTE.NP( J) =Y1 (J) ♦CELT*EK1 ( J) /2. 

50 CONTINUE 

. . .TJ.?T.tOELT/2. . . 

CALL 0ERIV(YTEMF,EK2,TT,N) 

DO 60J=1,N 

„ ...YTENR( J).=Y1(J)J-CELT?.EK2( J)/2. _ .— 

€0 CONTINUE 

CALL 0ERIV(YTEHP,EK3,TT, K) 

..0Q..7.0J = 1,M _ - - 

YTEHP( J)=Y1 (J) +DELT*EK3 { J) 

70 CONTINUE 

T=T.+ DELT „ . .. ... 

CALL OERIV( YTEMF,EK4,T,N) 

DO F0J=1,N 

Y1_U.J=Y1.(J) + (E!<1{J)-+2«!L(EK2(J).+E!53(J))+EK4(J) )* (OELT/6. ) 

80 CONTINUE 
RETURN 

650 FORMAT (1H_ , 8E14, 6.) . 

END 


1 
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SUBROUTINE CcR I V ( Y . K ,T , N ) 

C 

-Z OERIV COMPUTES. FORCE VALUES FOR EACH. .COUPLING AND ACCELERATION 

C EACH MASS 


-00 00.01 C0MMCS_EMASS(1.Q) j.F( 5 0,10.) ,X (5.0 .10) ,f.T.(5 0 ) .XSAV.(5 0J ,FSAV (50) , . 

1SLCPE(50) ,NCOUPT9(5 0,3) , IFL AG ( 50 ) ♦ DEL X( 5 0) , 0 Y ( 2 0 ) , NMA SST , NCOUP 
000007 REAL Y(20). K(20) 

000007 ,_00 5..I = 1,N ... .„ .. . 

000010 K(I)=0.0 

000012 5 CONTINUE 

0H0013 . 0C .7„ Isl.N . , . 

000015 K(I)=DY(I) 

000017 7 CONTINUE 

000021 . .00.10 J=l.NCOUP. .. . ... . 

000022 IF (NCCUPTB (J.l) . EO. 0)GC TO 15 

00 00 24 OELX( J) = Y( ( NCO'JPTB ( J , 2 ) ♦ NC OUPT 3 ( J . 2 ) )-D - Y( ( NC CUPT9 ( J , 1 ) ♦ NC 

1 „ TB (J.l) )-l) ... . . ... 

000034 GO TO 20 

000035 15 DELX(J)=Y( (NCOUPT0( J,2)+NCCUOT3 ( J.2) )-l) 

.000043 20 IF..(IFLAG{ J) .EQ..-DGC. TO 25. . . . .. 

000046 FT (J)=ENTP ( J) 

000051 GO TO 10 

000054 _. 25 FT( J) -AMAX1 (0. . FSAV { J )-SLCPE ( J) * (XSAV ( J)-DEL <.(J) ) ) 

000066 10 CONTINUE 

000071 DO 3C J=2.N,2 

00072 FIJ=0.0 

J00073 FJI=0. 0 

000074 . JT=J/2 

000075 DO 40 L=l.NCOUo 

000077 IF(NCCUPTB(L.l) .N£.JT.ANr.NCOU“T9(L,2) .NE. JT)GO TO t^O 

OOOlOe IF( KCCUCTB (L ,2) . GT. JT) GC TO 50 

000112 FJI = FJI*«?T (L ) 

000 114 GC TO t-0 

000114 50 FIJ=FIJfPT(L) 

000117 40 CONTINUE 

000122 DY(J-1)=Y(J) 

000125 CY (J)=(FIJ-PJI) /EMASS(JT) 

000131 30 CONTINUE 

000133 525 «^ORMAT(lH , 1 0 X , 4 (6X , E12 . 4 ) ) 

000133 RETUPN 

000133 ENT 


IS 
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OOOOOE 
00 00 Of. 

000005 

000006 

000010 

000015 

000016 

000023 

000034 

000036 

000040 

000042 

000046 

00004?' 

00.1052 

000052 


SU59CUTINF FLA6(Y,N) 

C 

c flag checks each coupling por unloaoing 
c 

real Y(20) 

- _ . . CO.HMON EM_ASS(10) , F (5 0,_1 0 ) » X (_i0 , 1 0) ,FT(50) , X_S A V( 5.C)., FSAV ( 5 0 ) » . 
iSLCPE(50),NCCUPTg(56 ,3) , IFL A G ( 50 ) ♦ OELX (5 0 ) ,0 Y (2 0 ) , NKASST, NCOUP 
00 10 J=1,N 

IFJNCgUfLT2( J»l) .EO. 0)GO_TO_15_ . . . 

COFLX = Y (K’cbUPTBC Jt2) +NCOUP TG ( J » 2 ) ) -Y ( NCOUPT 0 ( J , 1 ) f NC CUPTR 
1 ) ) 

.. GO_ TP_1 fl _ . . 

15 ODFLX=Y CNXCUOT3 ( J,2) ♦NCOUPTB ( J , 2 ) ) 

le IF (OOELX.GE. 0 ,CR. 1 FL AG ( J ) . EQ . - 1 ) GO TO 20 

I.F.LAGJ J) =-l 

FSAV (J )=’fT ( j) 

XSAV ( J) =CELX ( J) 

20 IF(DELX(J).LE.XSAV(J))GO TO 10 

IFLACf J) =1 
10 CONTINUE 

_ .RETURN. _ _ __ 

ENO 


‘ page is 

POOR QUALITY 
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COLETA 


000003 

000003 

oooooe 

000007 

000010 

000013 

oooom 

000031 

00003<t 

0 0 0043’ 
000044 


FCFT^A^ 1,3 • SEHI-AUTO 9FL • fOl-lO-73) 


peal FUNCTICK ENTP(N) 

C 

_.C _ ENTP.INTERPCLATES FORCE FOR .SOECIFIEO COUOLING . 

C 

COMM0^4 E'<ASS(10)|F(50»10)»X(50,10) .FKGOj^XSAVtEO) ,FSAV(50) , 

1SLOPE<50) ,KCOUP.T.e.(.50.,.3J.,.IFLAG([5.0).,OELX(5 0.)tOY(30) , ^HASST,NCOUP 

IF ( 6 eLX(N) .GE.X fN,!) )GO TO 10 
J = 3 

G.O_JC_30_ 

10 NXsNCOUPTB (N,3) 

00 30 J=3,NX 

.IF {OELX(N) .LT. X,(N,J) )GO TO 30 _ . 

'36 CONTINUE 

30 ENTP=F<N,J-l)*(CELX(N)-X(N,J-l»)*(CF(N,J)-F(N,J-l))/(X(N,J)- 

_1X.(N, J-l).))_ _ .. 

RETURN 

END 



120 120001.00 

11.0 5 <> 010.10 
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Cl. 4 HANDLING PROGRAMS (HYSIM and HYSSG) 

1. Program Description 

HYSIM and HYSSG are programs designed to provide a coarse 
Indication of how changes in mass distribution, vehicle dimensions, 
and tire characteristics affect handling properties. HYSIM 
simulates the transient response of an automobile to a step In 
steer angle, and HYSSG computes the steady-state gain of yaw 
rate to steer angle as a function of forward speed. Both programs 
are based on a mathematical model which represents the automobile 
as a single mass with 3 degrees of freedom: lateral displacement, 

yaw and roll. 

In HYSIM, the dynamics of tire force buildup are neglected, 

which means the front and rear tire force equations become 

algebraic expressions for the tire forces Y. and Y . Thus, if 
. f r 

Y 

Y - i and Z - Y 

the above differential equations can be written in the form 
AY * BY + CZ 
Z - DY + E 
Y - A~^((B + CD) Y + CEO 

This is the form in which HYSIM solves the differential 

equations of motion for the vehicle. 

• ••••• • 

For HYSSG, y^^, ♦ , 0 , Y^, and Y^ are all set equal to zero 

and the steady gain of the vector Y, (and hence of its component 

0 , the yaw rate), to s , is computed as follows: 

1 -1 

BY - -CZ « -C (DY + ES), or « Y - - (B + CD) CE 
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Su.rr\mari^ o-f ^.cjoa'i (or \5 -P^r C/reo-T^ ^YiS'idn’i -5p<2ed /hodcl 


LaJerdl : 






cf,. a 

r»fz. 


0 


= ir^ y, 1- ff, (Ah^ + - h^)<p *■ I Yf - b 


'^oH: 

m. Ah, <y; - d !p * (i^ *m,Ah^)'S = 

PC 

-rn,Ah,x,v - (<L^^ ■>■ C^,) i> ^ {Yn,<^Ahr- K^^-K^,)4‘ 
-(h,-le)Y^ - (h,-xbe)Y, 

Fro n’t Tirc.^ 

%-Y, = -Yf- z. V, 

X, X, 

near TLrcs ■ 

Sit =-Yr- d^rZ, U - C. 'o f 

')/SC£, //dLTidxO O.S : 

Pp^~ 0> (o!Q5 ^ 

z^'-rn^o, -z^ 



/y)^ - Fgh^ 
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melons 

CLNTT5 

Dc’i'. ' Cl 1 7 1 ^ 


'Pf-^Qf f\>r, 


/VG.A51 

— 

A'amfitil «r>F GPzr : 


N5PD 

— 

A'annfiEfe Of SPE-ZV:' 


NSTLPS 

— 

NLLmeZR. OF 1 PlECr f)l 2 OR 5‘iEP5 


NP5TEP 

— 

'PRinr fRLQllbpc^ 


belt 

S£.C. 

nme tngrepient 


5P££.dO&) 

Mrr\ 

SPEEDS 


Try) 

Y(^) 

5B(L 

±^Ve^HL '^VE^OC- 2 1 i,\ y'nul RPrE FdLL BPl e 
" l^ouL fiitiC»L£ 

m 

£/Y) 

KG 


niRSE or \YFH2CLS 

Ir. 


KG-m^ 


rpomET^TS ofjNeRtjh) o'-2 t h Respect 

Ir, 

r„ 

€hRT<i> 

enRTz. 

KG-ro- 

MG-rr\^ 


Jo Y.,7^ AND Z RYES 


Cxz, 

KC-m^ 

TRddu. ti of itiEd Tj /? oJrrH Re5 pec 7 TO rsz 

1. 

ea 

m 

CJHE6tSP5E LENGTH 

t 

pi 

Pi 

length front Rule to dsNTEEof GC/iuiry 


H 

PL 

Roll dEHTEC HEIGHT 

dh^ 

Ch! 

Pi 

(CENiEG 0PGL7'iJ7i‘,' -Roll dFNifr J NEjHij 7 

Q 

THtrP 

— 

FdBC.TIoH cor. aN front U)NE£J- UJJ-5TBTX(L CoNvrTJoN: 

6 

E.PBLN 

Rpd 

inclination of Roll fitzs 


R 

rn 

Ill^.zri fr.Z'J'J: 


Hh 

Pn 

fiERo-B^NFrnlC GENteJL OF EFfoQJ 


CDH 

_ 

Drn.G Orzf ' ? "i': - F Hr n: r: HHe H 

cf 

BELTf) 

■Rpb 

STEEd angle 

C^f 

dVF 

nJ Fr.L 

Koi c Dfin-^Tl A -- * T rr ■' ' 


(LPR 

h.ijLf' z- 

Roll D/lmTjT\G- RffiT 

^(t>f 

£KVF 


Roll 512 if neh 0 - FLoN r 

^<t>r 

' CKVC 

NjrrL 

Rol ! 611 ; / NEL ■■ - RsH R 


ZPr 

p/ Cbl 

lfjeehl force -front 

C. 

Cjik 

nIrb-d 

c/} 7E£ n L For CL - PEF R 

cr 

eiGn'Pi^ 

m 

icE lay H Tj on LE NG 1 R - FCd 

ur 

\ 

\ 

‘zlpr'fip 


A'finr/rucN length - getf 






n^ 


C-140 


Vn^Zfi>6L^5 

UHli'Z 

IDF'-FF'I PI :> r>l\' 

£'?>nrrjoN 

'PPXr.f'/'t rn 

t 

1 

1 

ktCiP T 


L/nEjZFL Fddfi <c^' /n:iF/-, 

t 

f\IhvU,A') 

— 

rP)!tF5e OF nrr. 


Ar)^U,4.) 

— 

CcOE J f J ci£P'T fFni B'd y 

b 

s 

ni 

£L1~£L {6cc !;ti \-i ■n'rh'i 


mD(4,4) 

— 

Bff- m Fn Fi y v ED - n-\Fn l'J y 


’Bn)(4-,4) 

— 


(Ldefficj-f NT piFieiy 


'BrhT(3,3) 

— 


nnm foj CPC USED x/v niPTjBiy 


2nm(3,3) 

— 


[TNi/E£6ION Kdijul /\'£ 


C(4,4) 

— 

FlhV -mFieiY ^ 'h£V- rVP 7 /'■ 0 Y 


dfiFZF 

n/£P)D 

LfiTE^FL FoFt£^ CfiTF^n L UFD 0 PFdH7 JiCE. 


iFi-RZe. 


LF7€j2/)L FOF.CJL ^ LFTFMl COFO TTf>F. 

1 

DD 1:^,4) 

— 

d 0£FF I Cl F-^r /V p-j F'] ) 


SX,(4,a) 

— 

S£~F)iF)i e 1 y ^ 'UD-P-nPi CZX 


£ZL4,2) 

— 

(Ld'c.ffic:j£ni m/ndix 


S.F(4) 

— 

€E-rr\Fi F j F -v f-niFin y 

rrihK 

£mcF 

— 

FF)F5S ^ (C£0T££G£/h/rT9-'RoLL F£JG FT 


^rr)T{3,3) 

■■■— 

FPmCjy L'Dic IF n^Fip y 2 f'i/FF-idf 




nhrFTTDNFC ZNBFlTrF OF CZ hi EELS 


F(p) 

— 

3c FT T 3 CC " Y'l '}■' - ■ ■ ■ ■ 


FF 

— 

Ffi(Ui<B.F=^ STEFUD 

\0.^/3i(L^A 

rPCc^FF- 

— 

D.6/5 ^ DenG CjyFFFlCiFMT ^ PFdHtRl FCfP 

s 


HjFFrD 

GFf j: 7 V 


GC~ {4 ) 

— 

Pi NX'- Fnii C-] / -y r •" ■ r- ~ A v ; 

) 

Ki (id ) 

— 

. VFClAbLS.5 OF JNTSGf/'7J ’' 

[ 

5uBf)(3,r ) 

— 

\ny'ic:czt u-^F.z j-a' /) •■■■-' a/ . 


SuBRNV(3=) 

— 


^ 2 /i' yp ,L '"Z .'■ '. / y. ~ 1 /\ r 

i t 

1 

TRFS(3,3) 

— 



1 

J 

TFJ(%3) 

— 


\ 

y 

! 

XDor 

c. 

57A-rD Fop SPECIFJCISV nPEEC A^rj 7 

1 

ri ( li ] 

— 

y/i / J /) 6ls zj oh IN'ISOL F : j * •; 

Zp j 

! 


LPTE^PL LoPD on PdchT Tz££S 

i 

z(: j 

■ KG 

LPTE£.PL LDPD oh BEftP TieES 


lAL PAGE IS 
X p (JUATiTY 







~r r n 



















sr/» w r 


HVSIr^ (ma^IN F/^oG/^firri) 


C-1A2 


dcftme, 

VAArrt£L£5 


?AELIfnTHAe<f 

ComPuTftTroHs 


nCASE.^ 

AISPJD 


NSTe,7>S^ 

f^Psre.p 


D£LT, 

SPE£JDCr). j»/, 
r^sp6 




[ QomPture, 
ATf\-fnATlU% 
CaiFPt c,i enrs 


Att\- re\A7 e. t X 


\ . m»nnv \ 

/JuTa^sISWhT 
/ fl/nr B/nrr '/ 
f^fnTjiAr,T»(ih,3/ 


SP6£D 

0£P£AD£/vr 

Q.c>n\purftTltHS 


Com7>u.re j 

’bm-rnkTHlx , 

Coi-FPlCitMTi J 
(tpeto oe.vtnos.HT) 


Brn-n\PTA It 


naAst. /ispD, 
NSTEVi, 
hPS7£.P 


de.lt, 

SPEED (Z ). 

T'/, 


T» /, i 

; — 1 7~ 

3 Cfl« 3 £*-l 1 

NCASe 

3 / 

f' 

STOP 


VS-HtCLE 

' CHAgACXeg. IS ■ 
TICS 


yen I CCS 
t HA A ACTce lines 


Ainv-mATAiy 


dampurs. 

'brA-mPTA IT 

Co£FFie.iE/TTS 
(smo XHoeveNtuvi^ 


ComPu.r£, 

EE-mAT/EfY 

EotPF I ae NTS 



JSPDArl 

1 

HSPD 




tom PuTE 

DD-rnATAiX 

totfritiEHTS , 
(sptto VCPENOtHT) 


mATiriPi '\^ 


f£S,DO,eD, 
o-;a. '4 " 


S£oCj,k)= I 
Bm(T,K) T 1 

£0 ( J, K ) I 

\ ^ \ 

MATfnVY \ 

) i 

/ AiNV, BED, C, 

/ 


<irr,j), 3*/, 4, 

x*t ,4 


i'4- 
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C 

_C 

c 

000003 

O0JOO3 

0000(13* 
030003 
000003 
00*00 0 3 

C 

C 

000003 

000013 

000023 

000031 

003044 

000060 

00.0364 

000100 _ 

00 Oil's 
000121 
000133 

C 

C 

C 

J03137 

C 

C 

C 

n 

c 

000141 

00C156 

C 

C 

C 

000206 

C 

C 

000232 
000236 
003254 
003260 
00C303 
000334 
C00320 
00 0324 
000340 
00 0344 

C 

C 

C 


PRCGRAM HYSIP 

REAL AM(4,4) tRH(4,4) ,AINV (4.4).SU0A(3«3) ,SUBANV(3«3) 
REAL EE(4,2) ,00(2*4) «F(2) ,E0(4,4) ,EF (4) t 3EO(4,4) 

REAL RMT (3*'3) *D.*itT (3^3) * E*Nf (*373) .TRI (3,3')7 tRAN"( 3* 3) 
REAL SPEEOdO) 

REAL Vl(lO) ,ja(10)_.Y (^J 

COMMON *C ( 4 ,"4 ) , GG ( 4 )■*■ 

READ INPUT 


READ 410, NCASE,NSPO 
READ 410,NSTEPS,NPSTEP_ 

READ -4 00, DEL T 

RFAO 40C,(SPEED( J) * J=1,NSPD) 

REAR 400,TM, (Y(J),J=1,4) 

PRINT 500‘ ~ 

PRINT 498,NCA3E,NSPD,NSTEPS,NPSTEP 
PRINT 49g,_0ELTj_(^SPEE0JU) *J=1*NSPD) 
PRINT' "SOS" 

PRINT 520, (Y(J) , J=l,4) 

PRIM_T__^0 7 

CASE LCQo 

no 350 JCASE = l,NCASf=’ 

INPUT MASS PROPERTIES 


INPUT geometrical PROPERTIES 
READ 4 0 0,EM,ENRTX, EM.RTY, ENRT2,CXZ 

rear 4 0 0,EL1,EL, H,CH,THETA,EPSLN,R,H A,CnA,OELTA 

ttre properties* " - - ~ - 

RE AO 40 0,CPF,CPR,EKPF, EKPR,CAF,CAR,SIGHAF,SIGMAR 

p’Int vehicle characteristics 

PRINT 51.0 

PRINT 520,EM,ENRTX,FNRTY ,ENRT7,CXZ 
PRINT 530 

PRINT 520, ELI, EL ,H,CH, THETA, EPSLN 
PRINT R43 

PRINT 520, R, HA, CCA, DELTA 
PRINT 553 

PRINT 520,CPF,CPR,Ek:PP,EKPR 
PRINT 560 

PRINT 520,CAF,CAF,SIGMAF ,SIGMAR 
pRILTmIMARY COMPUTATIONS 


000363 


8=EL1-EL 


C^152 


GOLETft 


000362 

00364 

000366 


00 0‘370“ 
000371 
003372 
00 0'375 
000400 
000432 
00 0403' 
000405 
000406 
000407 
000413 
C03413 
000413 
000420 


000436 

000440 

000441 

30446 

J00451 

000461 

000463 

010464 

030474 

000476 

030476 

000503 

000507 

000511 

000514 


000532 

000534 

000535 

000543 

03C543 

000546 

030552 


'00553 

J0555 

030556 

0)0561 


PO«tTRftN 1,3 • SE“.I-®UTO <RFL ♦ fOl-lO-73) HYSIM 


PACdEF=0.6l2*COA 

G=5.507 

EHCH=E1*CH 

c ■ “ 

C A-HATRIX COEFFICIENTSISPEED INDEPENDENT) 

C 

*no 10 0'J=1,4 
DO 100 K=l,4 
AH(J,K)=0.0 
■"100 fONTiNUE 

AMCl, 1)=EN 
AM(l,3)=EMCH 
'AH"(2f'2)=ENRTZ" 

AM (2,3)=-CXZ 
AM(3tl)=_EMCH 

“ 'amis'.zTs-cxz ■ 

AM( 3» 3) =ENRTX^EM*CH**2 
4M (u,4)=i,c 
■■ PRINT 570 

PRINT 5fl0, ( (AM(I,J),J=1,4) ,I=1»4) 

C 

c Inverse of a'-matrix 
c 

DO 110 J=l,3 

DO 110 K=l, 3 
SUFA ( J,K)=AM U,K ) 

110 CONTINUF 

CALL VATINV( SUnA,SU8ANV,BMT,nMTT,EMT,TRI ,TRAN,3) 
DO 120 J=l,4 
DO 120 K=l,4 

IF ( J.LE.3. AN0.K.LE.3)G0 TO 13C 
AINV (J,K)=C ,C 
GO TO 120 

130 ATNV (J,K)iSU8ANV (J,<) 

120 CONTINUF 

AINV(4,U)=1.0 

■ ■ PRINT‘59C - - 

PRINT 5 80 , ( (AINV (I, J) ,J= 1,4) ,1=1 ,4) 

C 

C 3-MflTRix COEFFIC IEMTS(SPEED INDEPENDENT) 

C 

no 140 J=t,4 
DO l*f0 X=l,4 
8M ( J,K) =0. 0 
140 CONTINUE 

RM( 3, 3) =-CPF-CPR 
BM ( 3,4) =E ')CH»G-EKPF-EKPR 
9M (4, 3)=1. 0 
C 

C E-MATRIX COEFFIC JENTSCSPEED INDEPENDENT) 

C 

on 150 J=l,2 
DO 153 K=l,4 
FE (j,«o=c.: 

150 CONTINUF 
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900564 EE(l,l)=i.O 

90566 EE(l,2)=l.O 

000567 ^EE(2,t)=EL 

00 0570 ~ EEI2",2V=-'¥ ' ' ' 

090572 EE (3,1)=-H+EL»EPSLN 

000575 EE (3,2)=-H-9*EPSLN 

0~00577 ‘oRiNT 490, ( (EECI, J) , j=i% 2)‘,I = 1,4V 

C 

C SPEED LOOP 

■ C 

000614 DO 390 JSPO=l,NSPD 

000616 JPLsC 

000617 TlsTH 

000621 DO 165 J=l,4 

000622 Y1(J)=Y(J) 

000624 'l'65 CONTINUE 

000626 XOOT=SPEED( JS=D) 

000533 PRINT 699,XOOT 

‘C 

C SPEED DEPENDENT COMPUTATIONS 

_ C _ 

003635 ■■■ ‘ Fr=FA’COEF*XDCf **2 ’ 

000640 ZF = THETA*£M*G-FA-^HA/EL 1 

000645 7P = EM»G-^F _ 

00 0 647 ■ ■ 'cAF'zFsC’AF^ZF' 

000650 CAP70=CAF*ZR 

10552 ENRTY9=enrTY/R*»2 

-00655 PRINT 619, ZF,ZR 

C 

C R-MATRTX COEFFICIENTSCSPEED DEPENDENT) 

C 

C006 64 9M(l,Z)=-X0nT-*(EMfENRTYR) 

00 0667 RM (2,l) = XDOT*ENRTYR 

00 0673 PH (2,4)=FA* (CHfH-HA) 

000674 0M(3,2)=-XDDT»EMCH 

•000676 PflNT 620 

00 0702 PRINT 580; ( (eMaVi) ,J= 1, 4) ,I = i,4) 

C 

C O-f^ATRlx C3EFFIC lENTS (SPEED DEPENDENT) 

C 

000720 DO 173 J=t,4 

000722 DO 170 K=l,2 

C39723 on tJ,K)=0. 0 

000725 170 CONTINUE 

000731 0n( I, l)=-CAFZF/XOOT 

000733 OD( l,2)=-CAFZF*EL/xdOT 

000735 CO (2, 1 ) = -CA>1ZR/XDOT 

030’’37 DD(2,2)=CARZR''8/XOOT 

C 

C COMPUTE C-MATRIX 

C 

-' 0741 CALL MAT-10Y CFE, OD, ED, 4, 2, 4) 

3745 DO 190 J=l,4 

000747 DO 160 <=l,4 

000750 RE0(J,t<)=3M (J,K) ■►EO(J,K) 


HYSIM 
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00757 180 CONTINUE 

. .00763 CALL HATMPY ( AINV » BED ,C ,4 ,4,4 ) 

00P7_6_7 PPIh^ejO _ _ _ 

0b07"73 “ 'print 5R0i { (C(I, j) ,J=1,4) ,I=1,‘4) 

C 

C F-HATRIX COEFFICIENTS(SPEEO DEPENDENT) 

c ' " ”■ ~ ' ■ 

' 001011 F(1)=CAFZF 

001013 F(2)=0.0 

■ ■ ■ c ■ - -- - - - - ■ 

C CO*<PUTE GG-HATRIX 

C 

001014 ■ CAIl HATMPY(EE,F,'EF,4,2,1) 

00 1020 PRINT 491, (EF(J) ,J=1,4) 

0010 32 CALL HATMPY ( AINV ,EF, GG , 4 ,4, 1 ) 

00i036 ‘‘ ""PRINT 492, (GG( J) 7j"=l,4) ' ' " 

001050 00 160 J=l,4 

00 1052_ _ GG( J) = GG(J) »CELTA 

001054 "leO CONTINUE - - - 

001056 PRINT 492, ( GG(J) , J=l,4) 

C 

C ■■ pr"int "headi ngs" 
c 

001067 _ PRINT 6tG 

c”" -- - ■ - - - - - - - - 

C INTEGRATE ECUATIONS 

' C 

U01073 KosTEP=l 

001074 DO 200 KSTE P=l, NSTEPS 

001076 CALL PKTTA 1 ( Y 1 , K 1 , Tl ,4 ) _ 

001101 ■ IF (KFSTFP..NE.KSTEP)'GO TO' 21C ■ 

00 1103 ACCLAT = XD0T*Yl(2) KKl) 

00 1106 PRINT 650 ,T1, (Yl (J),J=l,4) ,ACCLAT 

001124 ■ KPStEP=KPSTEF+NPSTEP 

00 1126 210 CALL R!<TT A 2 ( Y 1 , K 1 , OELT , T 1 , 4 ) 

001132 200 CONTINUE 

00 1135 ■ 30 0 CONTINUE ' ■ - - -• ---- - - 

031137 350 CONTINUE 

001142 STOP 

00 1144 !*03 FQPMAT (7F10 . 4) 

331144 410 FORHAT(2I10) 

00 1144 498 FORMATdH ,414) 

001144 499 FORMATdH ,10E16.6) 

001144 490 FORMATdH , * EE-M A TRI X»/.8 E 14 . 6) 

001144 491 FORMATdH , *EF-M ATRI X*/ 8 E 14 . e) 

001144 “ 492 FOPMATdH , » GG-'1 ATRI X*/ 8 E 14 . 6 ) 

001144 500 format ( IH l, * VEHICLE CHARACTERISTICS*) 

001144 505 POPHAT dHj , * Y-VECTOR*) 

031144 507 FORMAT ( 1H-, ‘FRAME*) 

001144 510 FORMA T ( IHO , 8X , *HASS* ,9X , * X-INER T I A * , 7X, * Y - INE RT I A * , 7 X , * Z-I NERTI/ 

16X,*y-Z PRODUCT*) 

1144 520 FORMATdH ,3E16.6) 

J1144 5.30 FORMAT ( IHO , *GECMFTRY ‘//6 X, *WHEEL8ASE* ,6X , *FR, AXLE CG*,6X,*RC HI 

1HT*,5X,*(CG - PC) HT.*, 9X, ‘THETA*, 11X,*FPSLN*) 

001144 543 FORMAT (IHO, 5X,*HHL. RAO.*,4X,*CG - AERO. LOAO*,5X,*CO AREA*,9X,» 
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GO LET a 


'0 llLi« 

OOI14V 

00J144 
00li^4 
00 1144 
001144 
00 1144' 
00 1144 
00 1144 
001144' 

00 1144 
001144 


'=’0«?TPAN 1.3 • SFMI-auTO «<FL • >01-10-731 HYSIK 


ILTA*) 

550 FORMATdHO, *SUSPENSION*//10X,*ROLL DAMPI NG* , 2 3 X , *ROLL STIFFNESS 
1X,*FR0NT»,12X,*REAR*,11X .♦FRONTS ,12X,*REAR*> 

56 b' F6'RHAT'( IHO , •'TIS!e''PRbPERT'iES*'//'9X *L A TERA l' "FORC E 19X,*RELA XATIO 
1ENGTH*/8X, CFRONT*, 12X,*R EAR* ,1 IX , ♦FRONT* ,12X, *RE AR* ) 

570 FORMAT <1H0 ,*A-MATRIX*) 

’580 FORMATdH V4Er4.6) ' ” ~ 

590 FORMAT (IHO , *A-INVERSE MATRIX*) 

600 FORMAT (lH-v*SOEEO =*,E14. 6) _ 

'6 1'O FORMAT ('IH , *‘ZF’ =*", E14 . 6, 5y;*ZR =*,E14.6) 

620 FORMATdHO , *e-MATRIX*) 

630 FORMAT (IHO ,*C-MA_TRI_X*) _ 

64 0 F 0 R M A 'T ( 1 H C"/ 3 X , ♦ 1 1 M E ♦ » 6 X\ * Y I'O O't * , '7 X » * P S 1 0 C T * , 6 X 7* P H I 0 0 f * i 7 X , * P H I 
16X,*LAT. ACC.*) 

650 FORMATdH , E 8 . 2 , 1 0 El 2 . 4 ) 

ENO - ■■■ ■ 


PA(}r is 
■ ‘‘j AfJArjTY 
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000013 

00 00 13 

000314 
000016_ 
0'0"0017 
000023 
000324 
00 00 25 ■ 
000033 


000040 

000043 


000051 

000053 

000060 
10061 
-00066 
00 0070 
00 007U 
00 0103 
030102 
030134 
000114 
030120 
030124 

03 0127 
003133 

030143 

000151 

000153 

000163 

000173 

000176 

030201 


'10217 
. 0221 
00 0225 
00 1230 


(jp 
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pA 


subroutine MATINVl a, ATEMF.B.BTEMP.E. TRI, TRANc M) 

C MATINV OPTERf-INES THE INVERSE OF A REAL MATRIX 

C ATEMP IS_THE INJ/ERSE CF THE INPUT MATRIX A. 

REAL AfM.MYi^T'EMPlM.'M) ,3CM,M)V8 TEMP'(M.mT,TRI(M.M) ,TRAN(M,M) 
1E(M,M) 

N=M _ _ 

C ■ SfbRE'THE "iCENtltY MATRIX 
00 20 I=1«N 

DO 20 ^1,N 

IF ( I'LEQ. Jl'GO TO 10 


E(I,J)=C. 
G0_T0_2^ 
lb E(I,J)='l. 
20 CONTINUE 


C 

C iNit'tALIZc TR ANSFOR”mEC AND TRANSFORMATION MATRICES 


C 


C 

C 

C 


C 


c 


c 


c 

c 

c 

c 


CALL TOENKTRI, A,N_) 

C A L L ibEN T'( TR'AN » E , N ) 


TRANSFORM A-MATRIX TO UPPER TRIANGULAR FORM 


IMIN=N -1 
DO 125 I= 1 ,IMIN 

LOCATE nonzero ELEMENT IK I-TH COLUMN 
DO 10 0 J=I,N 

IF(TRI( J,I) .NE.O.)GO to 110 
10 0 COMTIN' e 
PRINT 500 
GO TO 200 

110 NOTZRO=J 

CALL I 0 ENT( 9 ,E,N) 

■?{I,I) = C. 

g (NCTZRO.NOT 2 RO) = 0 . ' 
g(NCT 7 RO,I ) =1 . . 

I .NOTZRO) = 1 . 

GET NONZRO FLEMENT Tn '(I , LO'CA TIO'N) 
CALL MATMPT (B.TRI. ATEMP, N,N,N) 

CALL MATMPY ( B , TR AN , BTEMP ,N,M,N) 

GPT ZEROES IN ( I ♦ 1 , I ) - ( N , I ) LOCATIONS 
CALL nENT(Q,F,N) 

IPLUS=I+1 

DO 120 J=IPLUS,N 

P( J, I)=-ATEMP(J, D/ATEMP <I,I) 

120 CONTINUE 

CALL MATMPY ( 8 , ATEMP, TRI, K,N,N) 

125 CALL MATMPY (B,OTEM?, TRAN, N,N,N) 

transform A-MATRIX to DIAGONAL FORM 


1.3 0 


CHECK DIAGONAL ELEMENTS 
DO 1.30 1 = 1, N 

IF (TRKI.I) ,EQ.3.)GO TO 140 

CONTINUE 

GO TO 150 


OiaGINAL PAGE IS 
' E POOR QUALITY 
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MATINV 


i 


000<33 140 WRITE(6,500) 

100234 GO TO 200 

000240 150 00 170 I=1,IKIN 

C GE'f zeroes"! N UPPER TRIA^GLE' 

) 000242 NMI=M-I 

000244 CALL IOENT(B.E,N) 

" 00024S" ■ 00* 160 J=1*NHI 

') 000253 R(J,Nfl-T)=-TRI( J,N + l-I) /TRI(Nfl-I,N + l-I) 

000274 160 CONTINUE 

.~00‘0277 "“ " ■ CALL MATMPY (B.fRI, ATEMP* N ,N»N)' 

000303 CALL M ATMPY ( B, TR AN , BTEMP , N, N, N ) 

000313 CALL lOENT ( TRI, ATEMP, N) 

00 03 21 i‘7 0 CALL ■lOETjT('TR'ANVBTEMP»N) 

C 

C NORMALIZE THE DIAGONAL ELEMENTS_ 

■ ■ c ■ ■ ■■ 

000335 CALL IOENT(BtE»N) 

000337 DO 190 J=t,N 

000344 '" ■ 3( J, jr='i.7TRI(J, J) 

000354 IBO CONTINUE 

000356 CALL MATNPY ( 0, TRAN, ATEMP ,N,N,N ) 

000363 ■ '2 0 0 RETURN 

000364 500 FORMATdHO, 19HMATRIY IS SINGULAR) 

000364 ENT 
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000006 

000006 

C 

C 

C 

SU3P0UTINE ICENT<A,8,N) 

rDeNT ecUATFS_ MATRICES _ 

REAL A (N,N) ,R(N,N» 

DO 10 1=1, N 

- 

. 

000007 


no 10 J=1,N 



000010 


A (I,J)«R(I, J) 



000017 


10 CONTINUE 



000022 


RETURN 



000022 


EMO 
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SUBROUTINE ► ATMPY ( A.B.C. NRA, NCANRB.N'CQ) 

C 

_ C MAT MPY m ultiplies MATRICES _ __ 

C 

OOOOll REAL A (MRA,NCANRP) ,0 (NCANR0,NC0) ,C (NRA.NCB) 

OOOqU DO 10 Jsl.NRA 

00 0012 OO' 10 K = 1,NCB 

000013 C(J*K)=0. 

00 0^7_ _ DO IG I=1 ,NCANRB 

000020 ■ C‘(J.K)=C{J,K) >A( jriT*B(l'i'0 

000035 10 CONTINUE 

0000A4 RETURN 

c ■ ~ 


000044 


ENO 


t 


GOLETA 


^ 0007 
u0 3007 

00 0 an 
000011 
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SUeROUTINE RKTTAKYl.EKl fTtN) 
REAL VK20) ,EK1( 20) 

CALL DERIV ( Yl,cKl,T,N) _ _ 

RETURN 

ENO 
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GOLETA 


030010 

003010 

00 boil 

000017 
0000Z2 
’boob '24 
000027 
0q_00 34 
00 CO 41 
000045 
000050 
000055“ 
000061 
000064 
00 0065 
000070 
000075 
bo 0107. 
000111 
000112 
030112 
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sue ROUTINE RKTTA2( Yl.EKl .OELT, f,N) 

C 

C_ P<TTA2 USEJ RU_NGE-<UTTA TECHNIQUE_ FOR NUMERIC_AL INTEGRATION 

c - — - - 

real Y1(50),Y2(50),YTEMP(50),EK1(50) ♦ EK2 ( 50 ) . E K3 ( 50 ) .E K4 ( 5 0 ) 

00 50J^l_.ji 

Yt"EMP“( J)=Yl (Jl+b'ELT^EKl ( j) /2. 

50 CONTINUE 

TTgT»OEL T/2 . 

'call 0ERIV( YTEMP,EK2,TT,N) * 

00 60JslfN 

Y T E M P{ J) -Y1 ( J) » 0ELT*EK2( J ) /2. 

■"'60 CONTINUE 

CALL OERTV( YTEHP,EK3,TT, N) 

DO 70J = 1^N _ 

■ YtEMPC J) =Y1( J) ♦b'ELT»EK3( J) 

70 CONTINUE 

T=Tj-nELT _ 

CA"lL 0'ERiV( YtEMP;iE!<4,T,N) " 

DO ftOJ=l,N 

Y1 (J)=Y1( J) ♦ (EKl (J) ♦•2.*(EK2( J) +EK3 ( J ) ) +EK4 ( J) ) •(DELT/8.J_ 

"a o' CONTINUE 

RETURN 

65 0 FORMAT n._H *_3§14i£» 

'en'c ■■ ■ 


"mAL P\CrE IS 
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SUBROUTINE OERIV (Y .K.T.N) 

C 

C OERIV SETS UP K-MATRIX FOR USE IN INTEGRATION ROUTINE 

c ■ ■ " 

000007 CO>*MON r.(4,A),GG(4) 

00 0 007 __£EAL Y(1 0) , K<10) 

000007 oo'io" 

0000 10 |((I)=GG(I) 

oooot^ 

00 00 14 ■ Ka)=K(r)Vc (I^jf^VfJ) 

000022 10 CONTINUE 

000026 RETURN 

00 00 26"' ENO 


VCMICLC CHARACTCRtSTICS 
1 3 sea f 

t>ooooaoE'a2 t.eaaeaoF«at 2.ooaoooE*ei j.aeaoeoEttt 
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<.3jr»a? a. j. r.aso2><aE>aa 

-a.;«.??aoE*a3 -6.^5Jlaa£•03 -s.azezsEE*:^ 
a. i.oaaiacEtja o. 



i 
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EOUlOTJOH 

TRoGflA ns 


PKin^Z 

_ 
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Nspd 


Nutm &EK C F SPEEDS 


S'PEEdIjj) 

^kcc. 

SPEEDS 


£/n 

KG. 

/7)A65 CF S^EMlCLLE 

^r)i, 

£N^TX 

KG- fTt- 

(rromeATs oFZ'NeK.TiA dzih Respect 


£P//?T9 
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\t6 )C^ pkd Z /UES 

^rt. 
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^CL 
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m 
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CdR 

— 
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DELTA 

Rf^D 
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^(pr 

CPF 
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d^r 
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£KPR 
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K^r 
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(LAF 

k/Kfd 
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Clr 

CA/9 

h'/F^n 
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C77 

SlGff\AF 

r'l 
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1 

1 
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: 
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Liwnz 
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— 
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— 
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— 

1 

Yr\fDiEl.aES US?D JN mFiE.iy J Ni/f CSJ ON 


BmTT(4-,4) 

— 

1 

'RclltinE 

d z 

£-p 

OAPZP 

hJ/Zad 



CAAU 




T>d(s,4 ) 

— 

c^£ FFLClfh'i /oim CI y 


euC^.A-) 

— 

£S~nr)Ri C2Y -V- DO- /mni Cj / 


£F(4) 

— 

E£-mAiCiv ^ F-mmczy 


ee (4,s) 

— 

CdEFFiCieriT mpi Civ 


£mT(A,A ) 

— 

nrtF^THiY. U51D 2 niPnCiy J NysEf/idh’ Ct.mint 

rh-Ah, 

Err^Cn 

— 

CnA55 ^ (CchlT<eJl GCA{JU ^ - BjiLLdENje ej//£J5//T 



KC-rn^ 
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— 

C/'S.f y L C'j f.y T / yrn cj } 


Ff\ 

— 

FAYLdSF Y: SVEED 

o.toi^C^ A 

PAtbEF 

— 

O.G/cl^ DCfG dFEFFXCSnErFr ^ FCdHrAi ACea 

s 

- 

h/j/^RD 

,Qcnyyi ^ 

• 


— 
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— 


^Folui ne 

X, 

IDdT 

'^'Isez. 

yeAlCLt SPEED FOC SPfC'j FlFV SPEED L^oi- 


V/o) 

' 
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z. 
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ZE 

KG. 

CRjejERl Loro on ^erc. tires 









FOR':^'^AN CODING FORM 


C-167 







































C-170 


GOLETA FO^T^AN 1,3 • SFMI-AUTO CFL • (01-10-73) 


C 

_ —9 

C 

000003 

000003 

000003 

000003 

0^0 00 0 3 

C 

C 

000003 
0000 13 

C 

C 

C 

OOOO^o 

"C 

c 

c 

000030 
0300 33 

C 

C 

c 

'’000 37 

C 

C 

030055 

C 

C 

C 

000135 

C 

C 

00 01 31 

0001 35 
303153 
000157 
000177 
003203 
003217 
000223 
00 0237 
0 ) 02«*3 

C 

c 

c 

00 0257 
000261 
000263 
0265 

C 

C 


POOGRA?< HYSSG 

OFFINF VARIABLES 

RFAL flH(4,4 » ,9E0INV (4,4) 

^AL E*: (4,2^)_,£gjJ,4J^F|2J jJJ)J4_^4J_,EFm f.3ED(4_,4) 

REAL ■RmT( 4»4) ,BMTT(4.(») , E*MT(4,4) , fRi (4,‘l* ), TRAN (4,4) 
REAL SPEE0(12) 

9 . 9 At .. t ® I » A .1 

READ INPUT 

REA0~4“1 OT~Nl:'/rSETNSPO 

READ 40C, (SFE£D( J) , J=l,NSPO) 

CASE LOOP 

on 350_ JCASEsltNCASE 
PRINT HEADINGS 

PRINT 500 

PRINT 507 

INPUT HASS PROPERTIES 

RFAO 430,EM,FNRTX,ENRTY, ENRT7,CX7 

INPUT GEOMETRICAL FRQOERTIES 

read 4 0C,EL1,EL,F,CH,THE 1 A , FPSL N , R, H A ,CO A , DEL TA 
TIRE- oROPERTIES 

PEAQ 4:o,CPF ,CPR,£KPF,EK PR, CAF, CAR, SIGMA F, 3IG. MAR 

PRINT VEHICLE CH ARACTERISTICS 
PRINT 510 ■ ■ ■ * ■ 

PRINT 520,EH,FNRTX,ENRTY ,ENRTZ ,CXZ 

PRINT 53: 

PRINT 520, ELI, EL ,H,CH, theta, EPSLN " 

PPINT 540 

PRINT 323, R, HA, CDA, DELTA 
PRINT 550 

PRINT 520,CPF,CPR,EKPF,EKPR 
PRINT 560 

PRINT 520,CAF,CAF,SIGMAF .SIGMAR 
PRTLIMINARY CCMPUTATIONS 
R=FL1-FL 

FACCEF=0 .61 2*COA 
G = ‘3..4 0 7 
EMCH = E'^*CH 


B-MAT?ix COcFFICI£NTS(5PEEC INOEPENOENT) 
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GOLETA 


.00267 

000270 

000271 

000274 

000277 

000302 

000306 


000307 

0003il~ 

000312 

000315 

000323 

000322 

030323 

000324 

000326 

000331 

000333 


000353 
3 0352 
000354 


000361 
030364 
000371 
000373 
000374 
000376 
00 0401 


000410 

000413 

000414 

OQ0420 

000422 

000426 


000444 
^ 0446 
.0 C447 
00 0451 
003455 
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C 


00 140 J=l,4 

00 140 K=l,4 

BM(J,K)=0.b ~ 

140 CONTINUE 

_ _ BM(3,3)=-CP_F-CPR 

BM('5V4T=ENCH*G-E'kPF-EKPR 

BM(4,3)=1.0 

C _ _ _ _ 

C"' " E-HATRiX COEFFl'CirNTS CSPEEO INOE PENDENT )‘ 

C 

DO 150 J=l,2 

■ " CO 150 K=i, 4" ■ 

EE (J,K)=0.C 

150 CONTINUE _ 

'EE(i',r)=i'.i) '■ ■ 

EE (1,2)=1.0 
EE(2,1)=EL 

EE(2,2)=-B' ~ 


EE(3,1)=-H«-EL*EPSLN 
EE(3,2)=-H-6*EPSLN_ • 

Print" 49dr( (EE(iTJ)'t j= iTs) ♦'i=it4) 


V. 

c 


c 

c 

c 


c 

c 

0 


c 

c 

c 


170 


SPEED LOOP 

DO 300 JSPD= 1 ,NSPD 
XOOT=SPeEO( JS° 0 ) 

PRINT f,OC,XQOt 

SPEED DEPENDENT COMPUTATIONS 

FA=FAC0EF»XCCT»*2 

2 F=THETA*£M*G-FA*HA/EL 1 

ZR=cM*G-Zr 

CAFZF=CAF*ZF 

CARZP=CAR»ZR _ _ 

EMPTYR'=FNRTY/R *»2 

PRINT 610 , ZF,ZR 

9 - 1 ATRTX COEFPICIENTS(SPEFD DEPENDENT) 


3M(l,2)=-XnCT*(£Mf ENRTYR) 

BM (2,1 )=XQOT*ENRTYR 
9M ( 2,4) =FA* (CH+H-HA) 

EM ( 3,2)=-XD0T»EHCH 
PRINT 623 

PRINT 560, ( *BM(I ,J) ,J=1, :*) ,1 = 1,4) 


OlUGINALP^f;’; 
OF POOR 


0-MATRIX COEFf IC lENTS (SPEED OEPENOENT) 


00 170 J=l,4 
•30 170 K = l,2 

on(j,K)-=o.o 

continue 

30(1,1 )=-CAF7F/XD0T 


C-172 


GOLETA FOPTRAN 1.3 * SFMI-AUTO RFL • (01-10-73) HYSSG 


'00457 
00461 
OJ 0463 


00C465 
000471 
00 04_7 3 
00047'4 
000503 
00 0^7 
000516 
000522 


000540 

000542 


000543 
000547 
000561 
000563 
10565 
000567 
000572 
000575 
030601 
000603 
000606 
03 0607 
00061.3 
000625 
000630 
000632 
C0C634 
000634 
00063<4 
00C634 
00C634 
000634 
000634 

00 36.34 
000634 

000634 

'0634 

000634 


C 

C 

C 


C 

C 

C 


C 

C 

C 


on (l,2)=-CAFZF*EL/XOOT 
DO (2, 1)=-CAR2P/XOOT 
? (1? ? )=CARZR»B/ X p 0 T 

COMPUTE BEOINV MATRIX 


call" MAfMPYCEE,0D,ED.4,2,4) 
no 190 J=l*4 

DO _1^0 )^1 , 4 

BEb( J,K) = BM C'J.K) +Ed( J^K) 

180 CONTINUE 

CALL jIATINV ( BEOt OEDINV tBMT ,BMTT, EMT. TRI_. TR^N, 4) 

PRINT 590 


ORINT 5 80. C (DEDINV( I, J) . J = 1 , 4) , I =1 » 4 ) 


F-MATRIX ■COEFFIClENfStSPErO bEFENOENT) ' 

F(l)=CAFZF 
F(2)=0.0 


COMPUTE Y- MAT RIX 


CALL MATMPY (EE,F,EF,4, 2, 1) 
PRINT 4qi^(EF (J) ,J^^,4) 

no 190 j=l » 4 ' 

EF( J)=EF( J) ♦DELTA 
DO 193 K=l,4 

BEDINVfJ.K) =-"9EbtNV ( J.K) 

190 continue 

CALL MATMPY (RFOINV.EF.Y, 4,4, 1) 
CO 20 0 j=l,4 
VI (J)=Y( J) /CELTA 
200 CONTINUE 
OR I NT 64S 


PRINT 6?C‘ 'Y1 ( J> , J = l,4) 


303 CONTINUs 

350 CnNTINU; 

STOP 

400 FnRMAT{7Fx3,4, 

410 FORMAT(2IlC) 

490 F0RMAT(lHC,»tE-MATRIX»,2(/4E14.6)) 

491 FORMAT (IHO, • FF-M AT RI X * / 2 (4E14.6) 

530 format (IH l , ♦VEHICLE CHARACTERISTICS*) 
507 FORMATdH-, *FRAME*) 

513 FORMAT (IHO, 8 X , *M A 3S * , 9 X, ♦ X- INERT I A* , 7X, 


(h 


KxlNAL PAGE P' 
QUALITY 


l^nOR 


•Y_TMPDTT4*.7y.*7-TWFPTT 


16X,*X-Z p:?0OUCT*) 

520 FORMAT CIH , ,3E16. 6) 

533 FORMATdHO, ♦GEOHETRY*//6X,*WHEEL9ASE*,6X,*FR. AXLE CG*,6X,*RC M 
1HT*,5X,*(CG - PC) HT.*,8X,*THETA*,11X,*EPSLN*') 

540 FORMATdHO, 5X,*WHL. RAD.*,4X,*CG - AERO. LOAO*,5X,*CO AREA*,6X, 
lEER ANGLE*) 

550 FORMAT (1H3 , *SUSPENSlbN*//10X,*ROLL 0 AHPI NG ♦ , 2 3 X , * ROLL STIFFNESS 
IX, * FRONT*, 1 2X,*REAR*, IIX ,*FRONT* , 12X ,*RE AR*) 

560 FORMAT ( iHO , *TIRE P ROPERT I ES* // 9X , *L ATERA L FORC E * , 19 X, * RE L A XA T I C 
1ENGTH*/8X. ♦ FRONT -, 12X,*R EAR* ,1 IX, *FRONT* ,12X,*REAR*) 
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HYSSG 


T00634 

580 

FORNATdH , 

00634 

590 

FORHATdHO, 

000634 

600 

FORHATdH-, 

~ 0006 34 ■* 

610 

FORMAT d'H~ ; 

000634 

620 

FOPMATdHO, 

000634 

640 

FORMAT! IHO, 

‘00 0¥^4 

■ 65 0' 

'FORMAT Hh" ,■ 


t<E14.6) 

♦BED-INVERSE HAT^IX*) 

♦SPEED =*»E14.6) 

♦ZF' =♦; eI4.6, 5X,*ZR'=»',E14.6) 

♦e-MATRIX») 

4X,»Y100T*,7X .♦PSIOOT* .6X,^PHI00T»,7X,»PHI* 
10E12".4) 


00 06 34 E)40 


C-174 


GOLETA 


'odooTf 

OOJOl^ 

000014 

000016 

000017 

000020 

0000J4 

000025 

000033 


000040 

000043 


000051 

000053 

000063 

3C061 

00 0065 
000073 
000074 
000100 
000102 
003104 
000114 
000120 
000124 

0001 27 
003133 

000143 
000151 
000153 
000163 
000173 
003176 
00 0201 


10217 

J0221 

030225 

000230 
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subroutine ♦*ATINV( A.ATEMF .BtBTEMP ,E, TRI, TRAN, M ) 

C MATINV DETERKINES THE INVERSE OF A REAL MATRIX 

C ^4TEMP rS_'Q«E INVERSE_OF TH_E INPUT MATRIX A. 

R>AL“MM,M)'VaTEMP"(M,M)^b'(M,H) , 3TEMP(M,M),TRI(M,M) ,TRAN (M,M) , 
IE (M,M) 

N = M 

C STORE' THE"I0‘ENTITY MATRIX' 

00 20 I=1,N 

DO_ 20 J=1, N 

■ fF(l.EO.J)GG TO 10 

E(I.J)-0. 

GO TO 20 

' iO E('l,J) = i'^ 

20 CONTINUE 

C_ _ ___ 

C' ' INITIALIZE ■ TRANS FO'RMED A NO TRANSFORMATION MATRfCES 
C 

CALL rCENT(TRT,A ,N) 

CALL IOENT(TRTn,E,N) 

C 

C TRANSFORM A-.MATRIX TO UPPER TRIANGULAR FORM 

c' 

IMIN=N-1 

BO 125__I = 1, IMIN _ _ _ _ _ _ 

C ■ ' LOCATE NONZERO ELEMENT I'N I-TH COLUMN 
00 ICO J=I,K 

IF (TRI (J, I) ,NE.O . ) GO TO 110 
IOC CONTINUE 
PRINT ?00 
GO TO 203 
110 NOT 7PO=J ' 

CALL rnENT(!:,E,N) 

8( I,I)=C. 

3 (NCT.7P0,NGT7R0) =0 
3(NCTZR0,I) =1. 

«(I ,N0TZR0)=1. _ _ 

C ■ GET' NONZRO ELEMENT IN’ (1,1)' LOCATION) 

CALL MATMPy (e,TRI,ATEMP,N,N,N) 

CALL MATM«»Y(R,TRAN,RTEMP,N,N,M) 

C GPT ZEROES IN ( I , I ) - ( N , I ) LOCATIONS 

CALL TDENT (3,F,M) 

IoLUS=l4-l 

CO 123 J=IPLUS,N 

B< J, I)=-ATEMF(J, I) /A TEMP (1,1) 

120 CONTINUE 

CALL MATMPY (6, ATEMP, TRI, N,N,N) 

125 call MATMPY (3, BTEMP, TRAN ,N,N,N) 

C 

C TRANSFORM A-MATRIX TO DIAGONAL FORM 

C 

C CHEC>f OIAGONAL ELEMENTS 

00 130 1=1, N 

IF(T = I(T,I) .EQ.O.)GO TO 140 
130 CONTINUE 
GO TO 150 
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0002.13 
; 000234 

00^0 24 0 ^ 

) 000242 

0£0244 
00 024'e 
000253 
000274 
. 000277 
' 000303 
■ 00031 ^ 
0dd3‘21 


000335 
000317 
00 0 3 44' 
000354 
0003 ^ 
00 0 3 63" 
000364 
000364 
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140 W9ITe(6,50C) 

GO TO 200 

_150 DO 170 I=1,IHIN 
'C ' GET ZEROES IN "UPPER TP.IA'^GLE 
NMI=N-I 



00 16 d J=1,NKI 

B( s-TRI CJ,N + 1-I) /TRI(Nfl-r,N>l-n 

16 0 CONT INUE 

call' HATMPY'Cfl.TRI, ATE'MP, N,N,‘N) 

CALL MATMPY t 8 » TR AN , B TEMP , N, N, N » 

CA^ IOFNT(TRI ,ATE MP» Nl 

170 CALL "lOFNT ("TRAN, 0TEMP,N) 

C 

C DIAGONAL ELEMENTS 

c' ■ ■ ■ ■' 

CALL I0ENT(8*F,n) 
no 180 J'sltN 

B(J.J)=1.VTRI(J,JT 

180 CONTINUE 

CALL MATMPY ( BfTRA.NtATEMP ,N,N,N> 

-"'20 0 "RE TURN" 

500 FORMAT (IHO ,18HMATRIX IS SINGULAR) 

END 


MATINV 
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SU-TROOtinE ICENT(a,0,N) 

c 

C lOENT EOUaTES MATRICES 

c ■ 

OOOOQ6 REAL A ( N,N) * B ( N, N) 

OQ0006 00 10 I-ltM 

00 00"b7 00 10 J=i,N 

000010 A(I,J)=B(I« J) 

00^017 JL^ CONTINUE 

000022’ "rETUON 

000022 ENO 




GOLETA 


030011 

000011 

ooooiz' 

000013 
000017 
0 0 00 20 
03 00 35 
03004i* 
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SURPOUTINE yATMPY( A,B,C, NRA,NCANRQ,NC3) 

C 

c mathpy multiplies matrices 
c 

PEAL A (NPA* NCANRfl) . R(NCANRe»NC8) ,C(NRA,NCB) 
00 10 J=1,NRA 
DO IOK=1,NC0 
C(J,K)=0. 

00 10 I=1,NCANRB 
C (“J , K ) iC ( jTk ) ♦ a ( J , I ) * 8 ('I ♦ K V 
10 CONTINUE 

_ _R_ETURN 

CT 


OOOCLP 


ENO 


VEHICLE tH»F»CTERIStICS 
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ii 
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Cl. 5 Life Cycle Costs 

Two changi^s have been made in LYFECC since its summarization 

s 

appeared in the report from Phase I, "Mission Analysis and Per- 
formance Specification Studies Report". First, the method of 
computing the discount factor now matches '(ha one described in 
"H:,brid Vehicle Potential Assessment Interim Progress Report", 
Appendix C, Electric and Hybrid Cost Handbook, R. Left, S. Heller, 
Draft #5030-162, Jet Propulsion Laboratory, Pasadena, California. 

The second change involves the inclusion of replacement 
batteries down payment. This was omitted from operating and 
life cycle costs in the first version of the program. 


or p'xr. 


^ IS 
BBAUry 
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000003 

000003 


000003 

000007 

000015 

000023 


000025 

000027 

00037 


000041 

000043 

000044 

000045 

000046 

000047 


000051 

000071 

000107 

000125 

000140 

000144 

000164 

300170 


000206 

00210 

.00212 

000215 


PROGRAM LYFECC 


C 
C 
C 
C 
• C 
C 
C 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


IFUEL CAN EQUAL 1 THRU 5 MITH THE FOLLOW ING RESUIJ^S 
IF IFUEL=1, PETROLEUM' AND' ELECtRICITY 'ARE 'BOTH NOMINAL 
IF IFUEL=2, PETROLEUM IS >307. 

IF IFUEL=3, PETROLEUM IS 

IF IFUEL=4'r ELECTRICTTV~1S ^307. 

IF IFUEL=5, ELECTRICITY IS -107 


INTEGER KBR(ll)' 

REAL EKPY(11),RKF(11),TRCKF(11) , AOC(5,ll ) , YLCC 15 » 11) f OF( 11 ) t 
10AOC(ll),OYLCC(5tll)fV XT J1 lJ_t^FC KM C 1 1) tPFC KN( 11) « PF CJ<P ( ID « 
1EFCKM(11);EF'CKNT11) ,EFCKP(11),TFCK(5*11) ,TA0C(5) , A0CK(5) « 
1TLCC(5) ,ALCCK(5) ,TOLCC (5 ) ,0LCCK(5 ) 

PRINT HEADING " 

PRINT 500 _ _ _ 

READ 400,NCASE“ 

READ 400,NYEARS 
NYRM1=NYEARS-1 


DISCOUNT FACTOR 

DO 1 J=1,NYEARS 
OF CJ)=1. /(!.♦. 02) ♦♦(J-1) 

1 CONTINUE 

CASE LOOP 

DO 350 JCASE=1,NCASE 
TK = 0. 

BR = 0. 

DC' 5 J=1,NYEARS 
KB«(J)=0 

5 CONTINUE 

INPUT 

READ 41O,CV,HV*W0,PHE,PM ;POP 
READ 410,GGPK»GOPK.EKWHPK,80EM,BRK 
read 410,(\/KT( J) ,J = 1,N YEARS) ,AVKT,EKT 
READ 410, (RKF(J) ,J=1,NYEARS) 

PRINT 510 

PRINT BZOtCV.HVtHB.PHEtPy.POP 

PRINT 530 “ 

PRINT 520, GGP<,GDPK,EKHHPK,BRK,BOEM 

MILEAGE 

FACT=EKT/AVKT 
DO 10 J=1,NYEARS 
EKPY ( J)=FACT *VKT (JS 
10 CONTINUE 


C 
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000217 

000220 

000222 

000224 

000227 

000230 

000231 

000234 

000245 

000247 

000247 

000250 

000253 

000256 

000261 

000265 


000301 

000304 

000307 

000312 

000314 

100316 

J00317 


000324 

000327 

000332 

000335 

000337 

000341 

000345 

000352 

000354 

000357 


000361 
000362 
000363 
000371 
000377 
000405 
000413 
000421 
3 0431 
000432 
000434 


C SET BATTERY REPLACEMENT' FLAG " 

C 

00 15 J=1,NY^ 

BR = BR^EKf»Y ( J) 

TK = TK*-EKPV( J) 

IF(BR.LT.8R»^)_G0 _TO IS 

KBR(J)=1 

BRsBR-BRK 

15 CONTINUE _ _ 

IFCBR+EkPY(NYEARS)'.'LE.8R'i<.OR.PM.EQ.O.)GO TO 18 
EKPY(NYEARS) =BRK-3R 

BRsBRK 

GO TO i9 

18 BR*BR«-EKPY( NYEARS) 

19 TK=TK<-EKPY(NYEARS) 

PRINT' 533' 

print 534 TQ 

^ PRINT 590. <EKPT( J)_.J=W11J ,T!^_ ZX. 

C MAINTENANCE COSTS 

C _ 

HEMC=PHE/;746*.0'0310 6^.1 il818 
CMC=HV».00 00 06 +. 021742 

EMMC=PM/.746*.001 242».03 7273 

IF (PM.EQ.O".")EMMC=0 . 

BMC=W3*.000 248 
FMC=0. 

TMCK=HFMC«-CMC«-EMMC>8MC+F>C 

C 

C REPAIR COST 

C 

HERC=RHE/- 74 6*. 0 0497 + . 17 3939 
C9C=WV*.000 1 24 f .59015 
EMRC=PM/.746*.00 124f.0559l 
IFCPM.EQ.O. )EMRC=0. 

ARC=WW*.000 05 _ _ _ _ 

TRC=(PHE»-PM) /.74'6*i0 00 80'8*.031061 ■ 
TRCK=HERC»-CRC«-EMRC+ARC«-TRC 

00 20 J=1,NYEARS _ _ 

TRCKF(J)=TRCK*RKF( J) 

2.0 CONTINUE 
C 

C FUEL COST 

C 

IFCGGPK.EQ. OIGO TO 30 
00 25 J=1,NYEARS ' 

PFCKN( J) = C 20.90v.7479*( J-1) ) *GGPK 
PFCKP(J) = C27.17»-.9 723*(J-1) ) *GGPK 
PFCKMf J)=( 14 .63+. 52 35* (J-1) ) *GGPK 

EFCKN(J)=(4.23»-. 01U8*( J-l) )*EKHHPK 

EFCKM(J)=C3.81*. GO 97 * ( J- 1 ) ) *EKWHP K 

25 EFCKP(J)=(5.50+.014*(J-1))*EKWHPK 
GO TO 36 

30 00 35 J=1,NYEARS 

PFCKN ( J) =( 19 .36* .6599* ( J-l) ) *GOPK 
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LYFECC 


000<t42 
J004SQ 
000456 
000464 
000472 
000502 
0005Q6 
000526 
000532 
■ 000552 
000556 
000562 
000574 
000606 
000620 
000632 
000644 
000656 
000660 
000664 
000670 
000674 
000700 
000706 
000711 


000724 
000726 
000727 
000730 
000731 
000732 
0007 33 
000734 
000735 
000736 


00 0741 
000743 
000751 
000754 
000756 


000761 

000766 


PFCKP < J) s ( 25. 17». 0579* (J-D) ♦CORK 
PFCKMf J) = (13.55«-.4619*(J-1) )»GOPK 

_ EFCKNCJ) = (4.23«-.0108MJ-_1)_) *Ej<HHPK 

'■ EFCKM(jr=C3.81 + V0097*(J-l) )"*EKMHPK 

35 EFCKPCJI=(5 .50f.014*(J-l) ) *EKWHPK 

36 PRINT 540 

■ PRINT 55O»HEMC7CHC;EMMC,0M(r»FMC,TMCK 
PRINT 560 

PRINT 550,HERC,CRC,EMRC, ARC,TRC, TRCK 
PRINT 570 - 

PRINT 571 

PRINT 572, CPFCKN(J),Jsl,ll) 

"PRINT 573 ; (PFCKP(J) , J=l, 11) 

PRINT 574, (PFCKM(J),J=1, 11) 

PRINT 575, CEFCKN(J) ,J=1,11) 

PRINT '576; (EFCkP ( J)', j= 1,11) ■ 

PRINT 577, (EFCKM ( J) ,J=1, 11) 

00 37 J=l,ll 

TFCK(1,J) = PFCKN(J) fEFCKN (J)*' 
TFCK(2,J)=PFCKP( J) +EFCKN (J) 
TFCK(3,J)=PFCKM( J)+EFCKN( J) 

TFCK(4', J) = PFCKN( J) i-EFCKP (J) 

37 TFCK(5,J)=PFCXN( J) 4-EFCKM (J) 

PRINT 578 _ _ 

PRINT 584, (TFCK(1,'J) ,J = l,li) 

C 

C COST LOOP 

C " ■ ■ ~ 

DO 300 JC0ST=1,2 
TOCC=0. 

OOCK=0. 

00 38 J=l,5 
TAOC(J}=0. 

AOCKC J)=0. 

TLCC(J)=0. . 

ALCCKtJ)=0. 

TOLCCC J)-0. - -- -- 

38 OLCCK(J)=0. 

C 

C PAYMENTS -- - 

C 

PV=CV*2 

IF< JCOST.EQ.l)PV=7646. + 1.25*(CV-3 82.^. ) 

PV = PV«-PV*.05 

AL=PV-PV*POP 

AP=(AL + AL*. 065*4. ) /4. 

C 

C YEARLY DEPENDENT COST 

C _ - - - 

VOCl=33.fPV*.01*>l25. 
YOC2=33.»PV*.006+75. 

C - -- - . 

C BATTERY REPLACEMENT COST 

C 


000772 


0RC=2.*BOEM 
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000774 
00777 
001002 
001003 ' 
00 1005 
001010 
001014 ~ 
001020 


001034 
001036 ' 
001040 
001041 
001042 
001046 
00 1061 
001064 
001070 
001076 
001100 
001102 
00 1104 
001105 
001110 
01114 
.01120 
001122 
00 1131 


00 1135 
001137 
001143 
001147 


001151 

001153 


001157 
001161 
00 1165 
001174 
001200 
001201 
'0 1206 
.1222 
00 1231 


IF (JCOST.EQ. II BRC= 1.25*0 OEM 
BRC=RRC^0RC*.O5 

AOPBRsP.RC*POP 

ALBR=8RC-A0PBR 
AP8R=(ALBR+ALBR*.065*3.) /3, 
PRINT 585 

PRINT 586“ ■ 

PRINT 590,BRC«ALBR, APBRtBR 


C OPERATING'COSTS 

C 

DO 39 IFUEL=1,5 _ _ 

AOC(IFUE(.Vll=VDC'l 

YOCsYDCl 

00 40 J=2,NYEARS _ _ 

IF (J.QE;7) YOC'=YOC2 

AOCC IFUEL, J)=(TMCK+TRCKF (J) +TFCK ( IFUELtJ I > /lO 0 .* EKP Y (J )♦ YO C 
40 CONTINUE _ _ _ 

A OC ( IFUEL, N Y E AR S ) = AOC't IF U EL , NY E A RS) - Y DC2 

IF (NYEARS.LE.6) AOC ( IFUEL , N YE ARS) =AOC ( IFUEL, NYEARS) fY0C2-Y0Cl 

00 50 J=l, NYEARS . ... . 

IF (KBR( j| . EQ.O) GO' T0 "5d 
L=J+2 

00 60 K=J,L __ .. 

if(k.gt;nyears)go' to 50 ' 

AOC (IFUEL,J)=AOC ( IFUEL, J I ^AOPBR 
AOC (IFUEL, <)=AOC (IFUEL, K)+APOR 
60 CONTINUE 

50 T AOC (IFUEL) =T AOC (IFUEL) ♦•AOC (IFUEL, J) 

39 AOCK(IFUEL) =TA0C(IFUEL)/TK 

c ■ 

c oiscounteo operating costs 

C _ __ 

00 80 J=l, NYEARS 
OAOC( J)=AOC (1, J) *OF( J) 

80 TDOC=TDOCi-OAOC (J) _ ^ 

OOCK=TOOC/TK ■ 

C 

C LIFE CYCLE COSTS COMPONENTS 

C 

VSV=.01*PV 

BSV=.5*BRC* (l.-BR/BRK) 

C 

C LIFE CYCLE COSTS 

C 

00 85 IFUEL = 1,5 ' 

YLCC ( IFUEL, l)=AOC( IFUEL, 1) ♦PV* POP 

YLCC( IFUEL, NYEARS) =AOC(IFUEL,NYEARS)-BSV-VSV 

TLCC( IFUEL) =YLCC ( IFUEL, 1) ♦YLCC (I FUEL, NYEARS) 

00 90 J=2,NYRM1 

YLCC (IFUEL, J)=AOC( IFUEL, J) 

IF (J.GE.2. AND. J.LE. 5) YLCC (IFUELTJ )=YLCC (IFUEL, J) i-AP 
9 0 TLCC (IFUEL) =TLCC (IFUEL) ♦• YLCC (IFUEL, J) 
ALCCK(IFUEL)=TLCC(IFUEL) /TK 
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C DISCOUNTED LIFE CYCLE COSTS '* 

C 

001234 00 100 Jsl,NYEARS 

001235 ■ ‘OYLCCflFUEL, J)=VlCC( IFUEL,J)*OFC J) 

001243 100 TDLCC(IFUEL)=TOLCC(IFUEL)+OYLCC(IFUEL,J) 

001250 OLCCKIIFUEL)=TOLCC(IFUEL)/TK _ 

■001252 '■ 85 CONTINUE" ■ 

. C 

C OUTPUT _ 

* c ■■ *■■■ 

) • 001254 PRINT 535,JCOST 

‘ 001262_ PRINT 588 

001266 ■"PRINT '59 i 

001272 PRINT 590, ( AOC C 1 , J ) , J = l , 11 ) ,T AOC ( 1) , AOCKCl) 

001311 PRINT 595 

001315 PRINT "591 

001321 PRINT 590, f OAOC( J) ,J=1,11) ,TOOC,OOCK 

001337 PRINT 600 

001343 PRINT 591 " 

001347 PRINT 590, t YLCC ( 1, J) , J=1 , 11) , TLCC (1 ) , ALCCK ( 1 ) 

001366 00 120 TFUEL=1,5 

001370 PRINT 610"' 

001373 IFCIFUEL.EQ.DGO TO 125 

001375 IFCIFUEL.EQ.2)PRINT 579^_ 

001402 IF{IFI*EL-EQ"."3)PR”INT 58 0 

001410 IF(IFjeL.EQ.4)PRINT 581 

001416 IF(IFUEL.EQ.5)PRINT_582 _ 

001424 125 PRINT 591 

001430 120 PRINT 590, ( OYLCC C IFUEL , J ) , J=1 , 11 ) ♦ TOLCC ( IFIJE L) ,□ LCCK ( IFUEL ) 

001451 300 CONTINUE 

001453 350 CONTINUE 

001456 STOP 

001460 400 FORMAT(7I10) _ 

001460 405 FORMATdHO, 1014) 

001460 410 FORMATS7E10 .4) 

001460 500 FORMAT(lHlj*LIFE CYCLE COST ESTIMATION*) 

001460 510 FORMATdHl, ♦VEHICUE CH AR ACTERIST I CS* /9X , ♦COST'*Vi"l X , *HE IGHT * , 7X 

ITTERY MT.*, 3X,*HEAT ENG. RATING* , IX, *ELEC . MOTOR RA TING* , 2X , *X 
IN*) __ _ _ 

001460 520 FORMATdH ,8E16.6) 

001460 530 FORMATdH ,5X,*GAS LT/KM * , 6X, * 01 ESEL LT/ KM*, 8X , * KWH/ <M *, 4X , *9A 

lY REPLACE/KM*,4X,*eOEM*) 

001460 573 FORMATdHO ,*ANNUAL MILEAGE*) 

001460 534 FORMATCIH ,4X,*YEAR 0*,4X,*YEAR 1*,4X,*YEAR 2*,4X,*YEAR 3*,4X, 

IR 4*,4X,*YEAR 5*,4X,*YEAR 6*,4X^*YEA3 7*,4X,*YEAR 8*,4X,*YEAR ' 
1X,*YEAR 10*,4X,*TOTAL*) 

001460 535 FORMATdHO, *COST CASE=*,I2) 

001460 540 FORMATdHO, *MILEAGE DEPENDENT CO S TS C CEN T S/KM) * /2X , * M AI NT EN ANCE 

15X,*HEAT ENGINE*,7X,*CHASSIS*,7X, *ELEC. MOTOR * , 7 X , *8 ATTE RY *, 9X 
1YHHEEL*,06X ,*TOTAL*) 

001460 550 FORMATdH ,8F15.4) 

001460 560 FORMATdH , IX, *REPA IR*/5 X , *HEAT E NG I NE* , 7 X , *CH AS S IS * , 7 X , *ELEC." 

lOR* ,5X,* ACCESSOR IES*,5X, * TRANS MI SSI ON*, 4 X,*TOTAL*) 

00 1460 570 FORMATdH ,IX,*FUEL*) 

001460 571 FORMATdH*, 17X,*YEAR 0*,4X,*YEAR 1*,4X,*YEAR 2*,4X,*YEAR 3*,4X- 
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0014E0 

001460 

001460 

001460 

001460 ~ 

001460 

001460 

001460 

001460 

001460 

001460 

001460 

001460 

001463 

031460 

001460 

031460 

001460 


031460 

001460 

001460 

01460 


lAR 4*,4X,»YEAR 5*,4X,*YEAR 6*,4X,*YEAR 7*,4X,*YEAR 8*,4X,*YEAR 
13X,*YEAR 10*1 

572 FORMATdH , • RETRO ( NOrll N> * , IX , 1 2F 10 .4 ) 

573 FORMAT! IH ■; •PETRO( f3bX> *,2X, 12F10.4) 

574 FORMATdH , *PETRO (-30Z > ♦ , 2X,12F10.4I 

575 FORMATdH , ♦ELEC T ( NOMI N) ♦ , IX , 1 2F 10 .4) 

576 FORMATdH , •ELECT ! ♦307. )• ,2X,12F10 .41 

577 F0RMAT!1H , •ELECT ! -10%) • , 2X,12F10.4) 

578 F0RMAT!1H ,^TOTAL^) 

579 FORMAT !lHf;27X;^'!PETRdLEUM“ ♦307.) •) 

580 F0RMAT!lHf,27X,* !PETROLEUM -30%) •) 

581 F0RMAT!1H^,27X,^ !ELECTRICITY ♦30%)^) 

582 F0RMAT!1H^,27X,^ !ELECTRlTlfY -107.)*) 

584 F0RMAT!1H^, 13X, 12F10.4) 

585 FORMAT!1HO, •BATTERY REPL ACEHENT^ ) 

586 F0RMAT!1H~, 3X',^BR "COST^, 3X,*BR L OAN* , 2X, *BR P A' YM !♦ d X f IL E AGE- 
IT BATTERY SET*) 

588 F0RMAT!1H ♦•ANNUAL OPERATING COSTS*) 

589 FORMAT !1H , IIFIO .4) 

590 F0RMAT!1H , 1 2F 10 . 2, F 10 . 5 ) 

591 F0RMAT!1H ,4X,*YEAR 0*»4X,*YEAR 1*,4X.*YEAR 2*t4X,*YEAR 3*,4X,* 
IR 4*t4X,*YEAR 5*f4X,'*YEAR 6 • , 4 X , •YEAR' 7 • ♦ 4X , ♦ Y EA'R 3*,4X,*YEAR c 
1X,*YEAR 10*,4X,*70TAL*,5X,*PER KM*) 

595 FORMAT!1HO, •DISCOUNTED ANNUAL OPERATING COSTS*) 

600 FORMAT!lHfl, •LIFE~CYCLE~C13STS*') " 

610 FORMAT!1HO, •DISCOUNTED LIFE CYCLE COSTS*) 

END 
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"ENERGY AND MATERIALS - CURRENT MODEL AUTOMOBILES" 
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ENERGY AND MATERIALS - CURRENT MODEL AUTOMOBILES 
Materials Composition 

Excellent agreement has been obtained among the several recent studies 

n -5) 

which describe the materials composition of modern automobiles. ' ' Percentage 

compositions of ferrous metals, aluminum alloys, copper and copper alloys, 
zinc, glass, rubber, and plastics are changing slowly due to the gradual 
substitution of lighter materials. In any given year, however, the mass fraction 
of each particular material appears to be reasonably constant over a range of 
car weights. 

Table I shows the estimated weight breakdown for automobiles as estimated 
by a 1975 study. The same table gives the author's estimate for 1979 
model cars, in the form of a rounded-off interpolation. * These estimated 
1979 values will be carried forward to the manufacturing energy estimates 

Energy to Manufacture Automobiles 

The investigators Hirst and Herendeen have estimated that it requires 
123 million Btu to manufacture a typical American automobile. This figure 
is similar to a value of 126 million Btu cited by Berry and Pels. McGowan 
and Kirchoff, by contrast, calculated the manufacturing energy per car to be 
only 34 million Btu. 


"k 


The author of this section was also a principal contributor to the materials 
accounting in Reference (2). 
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Table I - Weight Percentage Materials Composition of Automobiles 
(Based on Estimates from Reference 2 for 1975, 1980', 
and 1990 and an updated estimate by the author for 1979) 


Author's 


Material 

Est. for 
1975 

Projections for 
1980 1990 

Est. for 
1979 

Low carbon and alloy steel 

61.2 

56.9 

54.2 

60.0 

Cast and malleable iron 

16.2 

13.6 

7.9 

14.0 

Aluminum alloys 

2.9 

6.3 

11.9 

4.5 

Copper and copper alloys 

1.0 

1.0 

0.6 

1.0 

Zinc 

0.8 

0.4 

0.3 

0.6 

Lead 

0.7 

0.8 

0.7 

0.8 

Other metals, incl. magnesium 

0.3 

0.7 

1.4 

0.5 

Rubber 

4.6 

5.1 

5.0 

5.0 

Glass 

2.4 

2.6 

2.8 

2.6 

Plastic 

3.5 

6.7 

9.2 

5.0 

Other non-metal 

6.5 

5.9 

6.0 

6.0 

Totals 

100.1 

100.0 

100.0 

100.0 
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Energy requirements to produce basic materials and components has not 

been a subject of close agreement. References differ, of course on the 

% 

basis for the estimates; some include only the energy for smelting or processing, 
while others include the activities of mining, ore transportation, rolling 
and component fabrication. One must also distinguish between electrical 
energy and powerplant fuel consumed. 

Table II identifies manufacturing energy estimates from several sources, 
and also sets forth the figures which will be assumed for this present task. 

These last figures represent our best present estimates of total energy typical 
for components made of each material, from mine or petroleum feed stock to 
finished parts. Actual total energies can vary widely with the characteristics 
of the ore, transportation, method of fabrication, and percentages of scrap 
vs. virgin maUrials used. 

A final step is to do a trial accounting of the total energy to produce a 
• complete automobile. We have taken as a 1979 baseline a U.S. -produced 
automobile of 3,620 lb. curb weight, such as the Ford LTD 4-door sedan. This 
is a six-passenger vehicle with a 302 CIO V-8 engine and a wheelbase of 114.4 
inches. The results of the energy accounting are given in Table III. A 
subtotal of 122.6 million Btu is the estimated energy to produce the parts from 
which the vehicle is assembled. To this must be added the energy to transport 
the components to the assembly plant, and then the energy to assemble and 
paint the vehicle. If subassemblies and parts equal to the car's total mass 
were transported an average of 100 miles, the fuel energy for truck transportation 
would be about 430,000 Btu. This is based on a referenced figure of 0.0176 
gallon of diesel fuel per ton-mile of freight hauled. 

A rough estimate of energy expended at the assembly plant can be based 
on an assumption that about 10 kW might be expended over a ten hour period per 
car. This is for conveying "and hoisting, welding, power tools, paint drying, 
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Table II - Estimates of Manufacturing and Processing Energy 
for Materials and Components Used in Automobile 


Production 


' 

% 

Material or Component 


Quoted 
Energy, 
Btu /lb 

Quoted 

from 

Reference 

Alloy Steel 


22,300 

(9) 

Stainless Steel 


34,000 

(9) 

Steel (Material Processing) 


21 ,500 

(3) 

Iron (Material ProcO'.sing) 


15,500 

(3) 

Steel (Material Processing) 


13,250 

(10) 

Steel Auto Hood (Mine to formed part) 


28,000 

(11) 

Aluminum Hood (45% virgin mat'l) 


108,300 

(11) 

Aluminum - Primary (Mat'l Processing) 


110,000 

(3) 

Aluminum - from Scrap (Mat'l. Processing) 

10,000 

(3) 

Aluninum smelting (Present processes) 


22,180 

(12) 

Aluminum smelting (New Alcoa process) 


15,350 

(12) 

Alum (Total energy of metal, from 50% 

bauxite) 

140,900 

(13) 

Copper (Total energy metal, from 1.5% 

ore) 

26,780 

(13) 

Copper (Total energy metal, from 0.6% 

ore) 

49,540 

(13) 

Copper and Alloys 


65,700 

(9) 

Zinc Castings 


45,500 

(9) 

Lead 


14,700 

(14) 

Plastics 


25,000 

(9) 

Plastics (Incl. feed stock) 


78,500 

(10) 

Fi^er Glass-Reinforced Plastic Hood 


40,100 

(11) 

Fuel Value, Resin Raw Mtl. 


20,000 

(11) 

Rubber 


37,000 

(3) 

Glass 


13,000 

(3) 


Value Used, 
This Study, 
Btu/lb 


45.500 

15.000 

78.500 

37.000 

13.000 


(Table II continued on next page) 
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Table II - Continued 


Material or Component 

Other Materials 

Stamped Steel Parts 

Cast Iron Parts, Incl. Machining 

Formed Aluminum Parts (45% Virgin) 

Finished Copper and Cu Alloy Parts 


Quoted Quoted Value Used, 

Energy, from This Study, 

Btu /lb Reference Btu /lb 

37,000 (3) 37,000 

28,000 

20,000 

108,000 

50,000 
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Table III - Estimate of Manufacturing Energy for an 
Automobile of 3,620 lb Curb Weight 
(1979 Model Year 4-Door Sedan) 


liaterial. Component, or Process 

Mass 

Fraction, % 

Mass, 

lb 

Specific 
Energy , 
Btu/lb 

Total 
Energy, 
Thous. Btu 

Steel Components 

60.0 

2,172 

28,000 

60,816 

Cast and Malleable Iron Parts 

14.0 

507 

20,000 

10,140 

Aluminum Alloy Parts 

4.5 

163 

108,000 

17,604 

Copper, Cu Alloy Components 

1.0 

36 

50,000 

1 ,800 

Zinc Castings 

0.6 

22 

45,500 

1 ,001 

Lead 

0.8 

29 

15,000 

435 

Rubber 

5.0 

181 

37,000 

6,697 

Glass 

2.6 

94 

13,000 

1,222 

Plastic 

5.0 

181 

78,500 

14,209 

Other Materials and Parts 

6.5 

235 

37,000 

8,695 

Subtotals 

100.0 

3,620 


122,619 

Transportation, Parts and Subassemblies 



430 

Assembly Plant Energy 




1 ,140 

Est. Total Manufacturing Energy 



124,189 


(Or Approximately 124 million Btu per Car) 
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and general utility supply to the plant. The electrical energy required 
is about 100 kWh or 341,200 Btu this is raised to 1.14 million Btu fuel 
energy input if the conversion of fuel to delivered electricity is 30% 
efficient. 

Our estimate of total manufacturing energy per car is 124 million Btu. 

It is interesting how close this independent estimate is to the 123 million 
Btu calculated by Hirst and Herendeen. 

The energy to produce batteries is of interest from the point of view of 
manufacturers and owners of electric or hybrid vehicles. From estimates 
given by Williams in Ref. (14), we calculate that the energy to manufacture 
an advanced lead-acid battery will be about 15,700 Btu per lb of battery 
weight. Williams has identified such a battery as having a specific energy of 
18 Wh/lb. It should be noted that since 15,700 Btu is the equivalent of the 
4,600 Wh of energy to manufacture a pound of battery, the battery would have 
to undergo the equivalent of over 250 complete discharges before the cumulative 
stored energy would exceed the original manufacturing energy. In other words, 
the manufacturing energy is very significant in relationship with the total 
energy the battery will store over its lifetime. 


THE SCRAPPING AND RECYCLING OF AUTOMOBILES: 
MATERIALS AND ENERGY IMPLICATIONS 
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Materials Recovery 

The materials composition of present and possible future automobiles 
was given in Table I, and is discussed in detail by References (2) and (5). 

Aoproximately ten years after the date of manufacture, cars of these compositions 
will appear in large numbers in the junkyards. 

Ford Motor Co. investigators have stated that the motor vehicle regis- 
trations are terminated on about 9 million cars/yr at this time, and that about 
90% of these are now being recycled for their metal content - primarily the 
ferrous metals. These authors further estimate that of the available 

automotive scrap (in 1974), 85% of the ferrous metal was recovered. Similarly, 
the recycling ratios for copper, aluminum, and zinc were 62%, 53%, and 27% respectively. 
It may be inferred that with currently improving techniques for the separation 
of shredded scrap, the recycling ratios for aluminum and copper will improve. 

Zinc, since much of it is in the form of galvanized panels and trim, will 
continue to be difficult to recover. 

Not much in the way of plastics is now being recycled, although in the 
future some economical methods of doing this may be developed. At least two 
alternatives are to somehow reconstitute the scrap by chemical means, or else 
to use the scrap plastic as fuel. The average heating value of such material 
is probably in the range 15,000 - 20,000 Btu /lb. A present objection to 
burning many plastics is the toxic nature of the combustion products. 

, The lead from storage batteries is eminently recyclable, and each automobile 

.v* 

uses about three or four batteries during its lifetime. Most of the replacement 
batteries are salvaged in the trade-in process, as are those finally ending up 
in wrecking yards. 
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Steps in Resource Recovery 

Usually the steps involved in recycling automobiles are: 

1. The vehicle is transported to a local wrecking yard (dismantler) . 

2. A dismantler salvages parts worth reselling, which is a larger 
fraction in the case of late-model wrecks than for worn-out obsolete 
cars. In the latter instance, the tires, battery, radiator, ard 
other easy- to- remove parts with a high scrap value are removed. 

Sometimes the engine and drive train are also removed at this point. 

The bulk is then flattened or compacted for efficient transport. 

3. The bulk is shipped to a processing center, where it is either baled 
or shredded. Shredding is becoming the preferred technique since it 
results in a higher grade of ferrous scrap and allows segregation of 
the other materials. Many shredding plants can take engine blocks 
and all . 

4. Recovered materials are shipped -to smelters or foundries. Not all of 
the material finds its way back to Detroit; much of the ferrous scrap 
(for example) goes to regional steel mills to be turned into reinforcing 
bar and structural shapes. 

Energy Consumption in Resource Recovery 

Transportation- The average number of ton-miles of transportation is not known 
to us, but one scenario might be the following: 

• The entire car, 3620 Ibis shipped 20 miles to dismantler. 

• After removing 300 lb of material, the remaining 3320 lb is forwarded 
» * 

y. to shredder, 100 miles away. 

• Non-metallics and some lost metal are landfilled near the shedder; 
assume 2,700 lb metal to be reclaimed is shipped 500 miles further to 
regional smelters. 
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• If most of the ferrous material goes to regional steel mills, this 
might still leave about 1000 lb. smelted metal to be sent another 
1000 miles to new users including the auto industry. 

The above shipments total 1377 ton-miles. Again using a payload-specific 

(15) 

fuel consumption of 0.0176 gal diesel fuel/ton-mile, the aggregate trans- 
portation fuel energy would be 3.27 million Btu per recycled automobile. Other 
assumptions or conditions could alter this outcome appreciably. 

Dismantler's Activities - Since dismantling is a labor- intensive activity, 
not much fuel energy is expended during this part of the process. Hoisting, 
flame-cutting, in-yard movements, and flattening of the hulk may consume 
4 kWh per car at 20 % efficiency, or a gross expenditure of 68,240 Btu per 
car. If 3,620 lb scrap is handled, the specific energy would be 19 Btu /lb. 

Shredder - Based on conversations with the operators of shredding plants, one 
automobile can be shredded in 15 seconds to one minute even with engine block 
and drive train still in the vehicle. This is a very energy-intense operation, 
since around 4000 hp is used to drive a hammer mill. Another 2000 hp or so 
is dedicated to segregation of the resulting scrap particles, and this process 
is expected to keep pace with the shredding. Our very provisional estimate, 
based on an assumed 20 seconds for a "car" to pass through each process and then 
adding 2 OX more energy for plant idling and overhead consumption, is 339,000 Btu 
per hulk. 

Total Energy for Resource Recovery - The above energy elements are summed up 
i* Table IV. The total energy for scrapping a 3620-lb passenger car and 
recovering an appreciable fraction of its resources is roughly estimated as 3.7 
million Btu. This does not include smelting or process energy to be charged 
against the next use. 
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Table IV - Provisional Estimate of Energy Required to Recover 
Resources from a Junked Automobile 
Basis ; Original curb weight 3620 lb, and assumptions 
given in text. These assumed conditions are 
subject to wide variations. 


Energy for Activity 

Activity Thous. Btu 

Shipping, aggregate of 1408 ton-miles 3,270 

Dismantling and Flattening 68 

Shredding and Sorting 339 

Total Energy Required 3,677 


(Approximately 3.7 million Btu. Does not include 
smelting or process energy for the next use.) 


Approximate Mass Accountability : 

Removed by dismantler for materials recycling - - - 200 lb 
Removed by dismantler for parts resale ------ lOO 

Metal salvaged by shredder ----------- 2,700 

Lost metal and non-metal to landfill ------ 620 * * 

Original Curb Weight — ------ 3,620 lb 


★ 




* Includes some eventual loss such as battery case materials. 

** Includes 5% loss of iron and steel, 25% loss of aluminum, and 12 lb loss 
of zinc. 
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The largest energy element by far in this total is transportation fuel. 
Even if this transportation component were assumed to be much l§ss (perhaps 
through the economies of rail shipment, rather than by truck), it would still 
dominate the recycling energy requirement. 

h ijcycling energy of almost 4 million Btu per car is to be compared with 
124 jsniion Btu to manufacture that car. While the 124 million Btu already 
includes the energy benefit of some scrap metal recycling, further study 
could reveal something of the sensitivity to ratio of recycled to virgin 
materials. 


1 
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